glume 40 September 1954 Number 9 


‘os 
NCHIGAN 


Proceedings of ifie*** 


National Academy 


of Sciences 
of the United States of America 





The Proceedings of the National Academy of Sciences 


OFFICERS OF THE ACADEMY 


DETLEV W. Bronx 
President 
GEORGE W. CoRNER 
Vice President 
JoHN KirkKwoop 
Foreign Secretary 


ALEXANDER WETMORE 
Home Secretary 


WILLIAM J. RoBBINS 
Treasurer 


EDITORIAL BOARD’ OF THE PROCEEDINGS 


Linus PAULING 
Chairman 


ALEXANDER WETMORE, //ome Secretary JoHn Kirkwoop, Foreign Secretary 
Epwin B. Wiuson, Managing Editor Wivuiam W. Rusey, Chairman of the 

National Research Council 
t. W. BEADLE GREGORY BREIT 
) Doisy L. C. DuNN 


J. W. BEAMS 
R. E. CLELAND 


( : 
] ; 

W. W. RuBey G. C. Evans W. F. Lipsy 
C 


HARLOW SHAPLEY >. H. GRAHAM J. voN NEUMANN 
F. E. TERMAN J. T. PATTERSON Pau, WEISS 

Subscription price is $7.50 per volume (12 issues); the price of a single issue is $.75. 
Make all remittances payable to The University of Chicago Press in U.S. currency or 
its equivalent. Authorized agent for the British Commonwealth, except North Amer- 
ica and Australasia: Cambridge University Press, Bentley House, 200 Euston Road, 
London, N. W. 1, England. 


Editorial correspondence should be addressed to the ProcerpinGs of the NATIONAL 
ACADEMY OF ScreNckEs, Attention Miss Mary D. Alexander, University of Chicago 
Press, 5750 Ellis Avenue, Chicago 37, Illinois. 


Business correspondence should be addressed to the Journals Division, University of 
Chicago Press, 5750 Ellis Avenue, Chicago 37, Illinois. 


Subscribers are requested to notify The University of Chicago Press and their local 
postmaster immediately of change of address. 


Microfilms of complete volumes of this journal are available to regular subscribers 
only and may be obtained at the end of the volume year from University Microfilms, 
313 N. First Street, Ann Arbor, Michigan. 

Entered as second-class matter June 21, 1920 at the Post-office at Easton, Pennsylvania, under the Act of August 


24, 1912. Accepted for mailing at a special rate of postage provided for in Section 1103, Act of October 3, 1917. 
Authorized June 18, 1920. 


PRINTED IN THE U. 8. A. 


THE PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES 
ts published for 
THE NATIONAL ACADEMY OF SCIENCES 
by 
THE UNIVERSITY OF CHICAGO PRESS 








NAL 
cago 


ty of 


local 


ibers 





TA 


re 





Proceedings of the 
NATIONAL ACADEMY OF SCIENCES 





Volume 40 + Number 9 «+ September 15, 1954 


GLUCOSAMINE KINASE OF SCHISTOSOMA MANSONI* 
By Ernest Buepine,t HowArpD RupPPENDER,t AND JOAN MACKINNON 


DEPARTMENT OF PHARMACOLOGY, WESTERN RESERVE UNIVERSITY SCHOOL OF MEDICINE, 
CLEVELAND, OHIO 


Communicated by Harland G. Wood, May 28, 1954 


The parasitic trematode Schistosoma mansoni is dependent for survival on ana- 
erobic glycolysis which occurs at an extremely rapid rate.! As in mammalian tis- 
sues, formation of lactic acid by homogenates and extracts of the worms proceeds 
via the Embden-Meyerhof scheme.2. Opportunities for inhibiting glycolysis of 
schistosomes without injury to the host would be available if the glycolytic enzymes 
catalyzing the same reactions in the parasite and in mammalian tissues were not 
identical. The existence of such differences has been observed recently. For 
example, the functional integrity of phosphohexose isomerase of Schistosoma man- 
soni is inhibited by an antibody which does not affect the activity of the isomerase of 
rabbit muscle.* The lactic dehydrogenases of schistosomes and of rabbit muscle 
can be distinguished from each other not only by their behavior toward specific 
antibodies* but also by their affinity to one of their substrates (pyruvate) and by 
the effect of pH on their activities.’ Furthermore, comparison of the properties of 
the hexokinases of the parasite with those of the host have revealed differences in 
the kinetics as well as in substrate specificities.* ® Sols and Crane’ have demon- 
strated that purified brain hexokinase catalyzes the phosphorylation of glucose, of 
fructose, of mannose, of glucosamine, and of other hexoses by adenosine triphos- 
phate (ATP). The same is true for crystallized yeast hexokinase.* ° On the other 
hand, Schistosoma mansoni contains three distinct hexokinases,* ® each one of which 
reacts specifically with only one hexose, either glucose, fructose, or mannose. None 
of these three kinases reacts with glucosamine.* In the present communication it 
is reported that the worms contain a fourth hexokinase which catalyzes the phos- 
phorylation of glucosamine by ATP. On the basis of these differences in the 
nature of the hexokinases of schistosomes from those of the host, the effect of 
glucosamine on the parasite was tested, and it was observed that this hexose is 
active against Schistosoma mansoni in vitro and in vivo. 

Glucosamine kinase and glucokinase activities were measured by incubating the 
enzyme preparation of the worms with the hexose and with ATP for a period of 10 
minutes at 37°C. The final molar concentrations of the constituents in the 
reaction mixture were as follows: glucosamine (or glucose), 1.5 < 107%; ATP, 
7X 10-*; MgCh, 1 x 10-?; potassium glycylglycine buffer (pH 7.1 for the deter- 
mination of glucosamine kinase and 7.5 for that of glucokinase activity), 5 X 1072. 
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ATP was added to the experimental samples before, and to the control samples 
immediately after, incubation. Proteins and phosphate esters were precipitated 
with barium hydroxide and zine sulfate,'° and the concentration of the nonphos- 
phorylated hexose was determined in the filtrate.'’ '* Hexose utilization was 
ealeulated from the difference between the free hexose concentration (reaction 
mixture incubated without ATP) and the experimental sample. Glucosamine-6- 
phosphate was prepared according to Brown.’ Glucosamine kinase was purified as 
follows: adult male schistosomes were homogenized in 0.01 M glycylglycine buffer 
(pH 7.5; 1 ml. of buffer was used per 140 worms), and the protein concentration of 
the homogenate was adjusted to 5.0 mg/ml by dilution with the buffer. 0.04 ml. 
of alumina gel C’, (23 mg. of dry weight per milliliter of gel) were added per milli- 
liter of homogenate. The mixture was stirred and then centrifuged. Addition of 
0.04 ml. of the same gel per milliliter of supernatant and subsequent elution of the 
residue with glycylglycine buffer (0.06 47; pH 7.5) yielded a fraction containing 
glucokinase activity. All operations were carried out at temperatures ranging 
between 0° and 4° C. 

Incubation of schistosome homogenates with glucosamine, ATP, and MgCl, 
resulted in the disappearance of the free amino sugar. This suggested that the 


TABLE 1 


GLUCOSAMINE KINASE AND GLUCOKINASE ACTIVITIES OF A HOMOGENATE OF Schistosoma mansoni 


GLUCOSAMINE KINASE GLUCOKINASE 
Per Cent Per Cent 
ADDITION Activity* Inhib. Activity* Inhib. 
None 0.21 Sas 0.65 
0.004 M glucosamine-6-phos- 
phate 0:12 43 0.28 57 
0.0016 M glucosamine-6-phos- 
phate 0.18 14 0.54 17 
0.01 M N-acetylglucosamine 0 05 76 Be ia 
0.003 M N-acetylglucosamine 0.13 38 


* Micromoles of free hexose phosphorylated per milligram of protein. 


enzyme preparation of the worms catalyzed the phosphorylation of glucosamine. 
Evidence for the formation of a phosphate ester of glucosamine under these 
conditions was obtained by the isolation of a water-soluble, alcohol-insoluble barium 
salt from the deproteinized reaction mixture. After the removal of the barium, 
analysis of this material for acid-stable organic phosphate and for glucosamine” 
as well as fractionation by ion-exchange chromatography" revealed that the com- 
pound had characteristics identical with those of glucosamine-6-phosphate.® 
These data indicated that the compound was not identical with glucosamine-l- 
phosphate. Because of the limited amount of available organisms, it was not 
possible to obtain this compound in sufficient quantities to determine whether phos- 
phorylation had occurred in positions 6, 3, 4, or 5. 

Glucose did not inhibit the phosphorylation of glucosamine, nor did glucosamine 
affect the activity of glucokinase. Therefore, neither of these two hexoses inter- 
fered with the phosphorylation of the other by schistosome homogenates. It has 
been reported that N-acetylglucosamine inhibits the phosphorylation of glucosamine 
by brain hexokinase.'* Similarly, N-acetylglucosamine reduced the activity of the 
worm enzyme which catalyzes the phosphorylation of glucosamine (Table 1). 
Crane and Sols have demonstrated that glucose-6-phosphate inhibits the activity 
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of brain hexokinase.” ° While this ester did not affect the phosphorylation of fruc- 
tose, it did inhibit the phosphorylation of glucose and of mannose by homogenates 
and by extracts of schistosomes.* Glucosamine-6-phosphate inhibited glucosamine 
kinase and glucokinase activity of schistosomes (Table 1); however, higher concen- 
trations were required to produce this effect than in the case of glucose-6-phosphate. 

On fractionation of homogenates of male schistosomes by the use of alumina gel 
(',, an enzyme extract was obtained which contains glucosamine kinase, free of 
glucokinase activity (Table 2). It has been found previously that purified gluco- 
kinase, fructokinase, and mannokinase of schistosomes do not react with glucos- 
amine.? Therefore, in Schistosoma mansoni phosphorylations of glucosamine and of 
glucose are catalyzed by two different enzymes. 


TABLE 2 


GLUCOSAMINE KINASE ACTIVITY OF A FRACTION OF A MALE SCHISTOSOME HOMOGENATE PURIFIED 
BY Use or ALUMINA Cy 


Experiment Glucosamine Kinase Glucokinase 
No. Activity* Activity* 
1 7.6 0 
2 5.4 0 
3 9.2 0 


* Micromoles of free hexose phosphorylated per milligram of protein. 


The occurrence in the worms of a glucosamine kinase whose activity is not 
affected by glucose raised the question whether a phosphate ester of glucosamine or 
a metabolic product thereof interferes with chemical reactions essential for the 
survival of the worms. Such a possibility was suggested by the inhibitory effect of 
glucosamine-6-phosphate on the glucokinase of the parasite as well as by the fact 
that glucosamine reduced the rate of glycolysis of intact schistosomes.* This 
problem was studied by determining the effect of glucosamine on the parasites in 
vitro and in vivo. In a chemically defined medium the schistosomes survive for a 
period of 5-6 days.'® Addition of glucosamine to this medium resulted in a signifi- 
cant reduction of their survival (Table 3). The action of glucosamine was not 
affected by the concentration of glucose in the medium. Accordingly, there was no 
evidence for a competition between glucose and glucosamine. Reduction of the 
survival of the parasite by glucosamine was observed also when horse serum was 
used as the medium (Table 3). 


TABLE 3 


KFFECT OF GLUCOSAMINE ON SURVIVAL OF Schistosoma mansoni IN VITRO 


Glucosamine Average Survival Reduction of Survival 
Medium (Mg/M1l) (Days) (Per Cent) 

ieee 6.2 oy 
Synthetic) 3.0 a 82 
1.2 2.3 63 
| 0.4 3.9 37 
( : 23.0 - 
Serum } 3.0 6.2 73 
1.2 9.8 58 
| 0.4 15.9 31 


Oral administration of glucosamine to mice infected with Schistosoma mansoni 
resulted in a change in the distribution of the worms in the host. The adult 
schistosomes are located in the mesenteric veins, in the portal vein, and in the liver 
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sinuses. When the worms are damaged by chemotherapeutic agents, such as organic 
antimonials, they lose their attachment in the mesenteric veins and migrate toward 
the liver. Such an effect was observed after oral or intraperitoneal administration 
of glucosamine; under these conditions a significant reduction in the number of 
worms in the mesenteric veins occurred which was associated with a corresponding 
increase of worms in the liver (Table 4). The differences in the distribution of the 


TABLE 4 


EFFEecT OF ORAL ADMINISTRATION OF GLUCOSAMINE TO MICE INFECTED WITH Schistosoma mansoni 
on DISTRIBUTION OF WoRMS IN PORTAL-HEPATO-VENOUS SYSTEM OF Host* 


Tora No. —PERCENTAGE OF Worms 1n——~ 
EXPERIMENT No. oF or Worms Me »senteric Portal 
No. Group Mice (Av.) Veins Vein Liver 
1 {Control 7 28.7 54 40 6 
| Experimental 13 28.7 30 41 29 
‘ {Control 11 25.9 47 41 12 
rs \ Experimental 10 24.7 32 38 30 


* Glucosamine was administered twice daily (10 per cent solution; pH 6.0) at a single dose of 50 mg. for a period 
of 14 days. This schedule resulted in no demonstrable acute toxic ‘effects on the host. 


worms in the mesenteric and portal veins and in the liver between the control and 
experimental groups were analyzed statistically by the use of the x-test. In both 
experiments the p-values were less than 0.001, indicating a high degree of signifi- 
cance. Quastel and Cantero” have reported that the intraperitoneal administra- 
tion of glucosamine resulted in a retardation of growth of sarcoma 37 in mice. It 
has not been determined whether a specific glucosamine kinase is present in these 
tumors and whether the activity of glucosamine against these tumors and against 
Schistosoma mansoni is brought about by the same mechanism. 

It should be pointed out that the selective toxicity of glucosamine for Schistosoma 
mansoni was detected as a result of observations revealing differences in the nature 
and in the substrate specificity of enzymes in the parasite and its host. It is 
reasonable to assume that further information about, and a better understanding of, 
the comparative biochemistry of parasites will result in the rational development of 
other and more effective chemotherapeutic agents. 


* This investigation was carried out with the support of a grant from the Division of Research 
Grants, National Institutes of Health, United States Public Health Service, and under a contract 
with the Surgeon General’s Office, United States Army. 

+ Present address: Department of Pharmacology, School of Medicine, Louisiana State 
University, New Orleans 12, Louisiana. 

1 E. Bueding, J. Gen. Physiol., 33, 475, 1950; E. Bueding and L. Peters, J. Pharmacol. Expl 
Therap., 101, 210, 1951; E. Bueding, L. Peters, S. Koletsky, and D. V. Moore, Brit. J. Pharmacdl., 
8, 15, 1953. 

2 EK. Bueding and H. Most, Ann. Rev. Microbiol., 7, 295, 1953. 

3 —. Bueding and J. MacKinnon, unpublished observations. 

4T. E. Mansour, E. Bueding, and A. Stavitsky, Brit. J. Pharmacol., 9, 184, 1954. 

5 T. E. Mansour and E. Bueding, Brit. J. Pharmacol., 8, 431, 1953. 

6 —. Bueding and J. MacKinnon, Federation Proc., 12, 184, 1953. 

A. E. Sols and R. K. Crane, Federation Proc., 12, 271, 1953. 

8 M. Kunitz and M. MacDonald, J. Gen. Physiol., 29, 143, 1946; L. Berger, M. W. Slein, 
S. P. Colowick, and C. F. Cori, J. Gen. Physiol., 29, 379, 1946. 

9D. H. Brown, Biochim. et biophys. acta, 7, 487, 1951. 

10 M. Somogyi, J. Biol. Chem., 160, 69, 1945; S. P. Colowick, G. T. Cori, and M. W. Slein, 
J. Biol. Chem., 160, 583, 1947. 

1M. Somogyi, J. Biol. Chem., 160, 61, 1945. 











DI 


NEW 


P 
whi 
this 
and 
sho 
mor 
side 
hem 
glob 
new 
was 
hem 


mole 
carb 
in tl 
posi: 
tive 
tive 
nega 
resid 
In s¢ 
the 1 
ize 1 


Su 
globi 
as, f 
high 
all e 
hem« 
boxy 





ALS, 


ganic 
ward 
ation 
er of 
nding 
of the 


ansoni 


ver 
6 
9 
2 
0 


. period 


1 and 
both 
gnifi- 
istra- 
oR 
these 
rainst 


soma 
ature 

It is 
ng of, 
ant of 


search 
ntract 
State 


Expil. 
macol., 


Slein, 


Slein, 











~I 
~I 


Vou. 40, 1954 BIOCHEMISTRY: SCHEINBERG ET AL. 


12 Z, Dische and E. Borenfreund, J. Biol. Chem., 184, 517, 1950. 
13 —), H. Brown, J. Biol. Chem., 204, 877, 1953. 

4. P. Harpur and J. H. Quastel, Nature, 164, 693, 1949. 

6 R. K. Crane and A. Sols., J. Biol. Chem., 203, 273, 1953. 

1 H. Ruppender and E. Bueding, unpublished observations. 

1 J, H. Quastel and A. Cantero, Nature, 171, 252, 1953. 


DIFFERENTIAL TITRATION BY MEANS OF PAPER ELECTROPHORESIS 
AND THE STRUCTURE OF HUMAN HEMOGLOBINS* 


By I. HERBERT SCHEINBERG, RutH 8S. Harris, AND JOAN L. SPITZER 


NEW YORK STATE PSYCHIATRIC INSTITUTE, AND DEPARTMENT OF MEDICINE, COLLEGE OF PHYSICIANS 
AND SURGEONS, COLUMBIA UNIVERSITY 


Communicated by J. T. Edsall, July 19, 1954 


Patients suffering from sickle-cell anemia possess hemoglobin (hemoglobin b) 
which differs from normal adult hemoglobin (hemoglobin a). A molecular basis for 
this difference was first demonstrated experimentally by Pauling, Itano, Singer, 
and Wells in 1949 by means of free electrophoresis in the Tiselius apparatus.! They 
showed that the hemoglobin b molecule possesses a net charge about three units 
more positive than hemoglobin a, within a pH range of about 1.5 units on either 
side of neutrality. Since that time, three additional abnormal hemoglobins, 
hemoglobin c,? with an even greater net positive charge than hemoglobin b, hemo- 
globin d*, and hemoglobin e* have been described. The present paper describes a 
new experimental method of differential titration by means of which information 
was obtained concerning the chemical basis of the difference in charge between 
hemoglobins a, b, and ¢. 

As the pH of a protein solution is increased, various groups within the protein 
molecule lose protons. Neutral terminal carboxyl groups and the neutral free 
carboxyl groups of aspartic and glutamic acid residues become negatively charged 
in the pH range from 2 to 5.5; imidazole groups of histidine residues lose their 
positive charges in the range from 5.5 to 8.5; terminal amino groups lose their posi- 
tive charges from 6.5 to 9.0; epsilon amino groups of lysine residues lose their posi- 
tive charges in the range from 9 to 11.5; neutral tyrosine hydroxy! groups become 
negatively charged in the range from 9 to 12; and guanidinium groups of arginine 
residues lose their positive charges above—often considerably above—pH 11.° 
In some instances groups may be tied off so that they are completely unreactive in 
the native protein,® 7 and the characteristic pH ranges within which the groups ion- 
izemay be modified by neighboring substituents.° 

Suppose that hemoglobin b possesses three more net positive charges than hemo- 
globin a by virtue of possessing three more positively charged nitrogen atoms such 
as, for example, those in lysine residues. In solutions with a pH of about 11.5 or 
higher the difference in charge between the two hemoglobins should disappear, since 
all epsilon amino groups are uncharged. On the other hand, if we assume that 
hemoglobin b possesses its greater net positivity by possessing three less free car- 
boxyl groups than hemoglobin a, then the difference in charge between the two 
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hemoglobins should disappear in a solution with a pH of about 3.5 or less, since 
most carboxyl groups are uncharged around this pH. Therefore, comparisons of 
the charges on hemoglobins a, b, and ¢ at selected pH values may be expected to 
yield information as to the nature of the groups responsible for the differences in 
charge when the comparisons are considered in the light of the pK values of the 
various groups. 

Paper electrophoresis is a simple technique for comparing relative mobilities* and 
therefore for comparing relative charges of proteins which are as similar in size and 
shape as hemoglobins a, b, and, presumably, c.''° If paper electrophoretic experi- 
ments are made on several paper strips, each with the three hemoglobins but each 
in a buffer of different pH, the differences in the distances which the hemoglobins 
migrate on a given strip should bear a close relation to the differences in charge be- 
tween the hemoglobins at that pH. If electro-osmotic flow can be measured*® 
and if, as seems reasonable, the hemoglobins do not bind different numbers of any 
small ions other than protons,! the ratio of the distances which the hemoglobins 
migrate should equal the ratio of the numbers of protons bound by each hemoglobin 
at that pH. Such experiments can be considered to be differential titrations of the 
three forms of hemoglobin and are similar to experiments performed by Durrum” 
in which histidine could be distinguished from arginine and lysine as well as from the 
neutral amino acids. 


MATERIALS AND METHODS 


Hemoglobins.—Blood was drawn into acid-citrate-dextrose solution from sub- 
jects who were homozygous for hemoglobins a, b, or c.!!~ Hemoglobin solutions 
were prepared according to the method of Drabkin,!? and were converted to car- 
bonmonoxyhemoglobin. 

Hemoglobin solutions for use in paper electrophoretic experiments were generally 
made isoionic by passage over columns of ion-exchange resins, according to the 
method of Dintzis.'* For determination of the isoionic points, 5~—7 per cent solu- 
tions of carbonmonoxyhemoglobin flowing from the columns were collected under 
mineral oil, and the pH was measured without added salt. All pH measurements 
were made in a Beckman Model G pH-meter at room temperature. 

Paper Electrophoresis.—Electrophoretic experiments on paper were carried out 
essentially as described by Kunkel and Tiselius.. Two spots of hemoglobin a were 
generally used, and these, together with two spots of a 10 per cent dextran solution, 
provided evidence regarding the uniformity of the electric field and an estimate of 
the electro-osmotic flow. Runs were made, at least in duplicate, at room temper- 
ature, for 1-4 hours, with a potential of 350 volts. A current of 7-10 milliamperes 
owed in most experiments. Papers were dried at 120° C. and stained with alco- 
holic bromphenol blue.’ 

The pH of the buffers was measured at room temperature (22—27° C.) before, and 
after, most runs. The change in pH of the buffer during a run was seldom more 
than 0.03 pH units. The ionic strength of all buffers was 0.06. 


RESULTS AND DISCUSSION 


The isoionic points!‘ of the three hemoglobins, measured, as described above, 1 
the absence of salt, were as follows: hemoglobin a, 7.11; hemoglobin b, 7.21; 
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hemoglobin ¢, 7.42. If the hemoglobins bind no other small ions present in the 
solution except hydrogen ions, each isoionic point should be equal to the isoelectric 
point as measured electrophoretically, except for the electrostatic change of iso- 
jonie point due to the ionic strength of the buffers used in electrophoretic experi- 
ments.'* In a solution somewhat more acid than the lowest pH of the three iso- 
ionic points, all three hemoglobins will move as cations in an electrophoretic experi- 
ment. Hemoglobin c, with the greatest positive charge, will move faster than 
hemoglobin b, while hemoglobin a, with the smallest positive charge, will move as 
the slowest of the three. In a solution somewhat more alkaline than the highest 
pH of the three isoionic points, the relative mobilities of the hemoglobins, now all 
anions, will be reversed, and hemoglobin c, with the greatest positive charge, or 
smallest negative charge, will move as the slowest. 

Table 1 gives the relative mobilities of the three hemoglobins in buffers of dif- 
fering pH. The pH values given for the isoelectric points (e.g., pH 7.0 for a) are 
probably no more accurate than 0.1 pH unit, since electro-osmotic flow was only 
estimated with dextran. Photographs of some of the filter papers of these runs are 
shown in Figures 1-6. 


TABLE 1 


RELATIVE MoBILITIES OF HEMOGLOBINS a, b, AND ¢ AS MEASURED BY PAPER ELECTROPHORESIS 
in Various Burrers (Ionic STRENGTH = 0.06) 


Buffer pH Relative Mobilities* 
Glycine-glycine HC] 3.0 +a = +b = +e 
Glycine-glycine HC] 3.4 +a = +b = +e 
Sodium acetate—acetic acid 4.1 +a = +b = +e 
Sodium acetate—acetic acid 4.8 t+a<+b<+e 
Sodium acetate—acetic acid 5.8 +a<+b<+e 
NasHPO,; -~-NaH2PO, 7.0 a =O + b < +¢ 
NasxHPO.-NaHePO, 4.2 —-ab=0+ce 
NasHPO,-NaH2PO, f 4 —a> —b c=0 
Sodium diethy] barbiturate—diethy] 

barbituric acid 8.6 —a>—b> —c« 
Sodium glycinate-glycine 9.8 —-a>-—b> —¢ 
Sodium glycinate-glycine 10.8 —a>-—-b>-e 
NaoHPO,-Na;3PO, M7 —a>-—b> —<« 
NaoHPOw-Na;PO, 12.0 —-a>—b>—-« 

* + denotes mobility as cation toward cathode; — denotes mobility as anion toward anode; > and < denote 


“mobility greater than’’ and ‘‘mobility less than,’’ respectively. 


It can be seen from Table 1 and Figures 1 and 2 that at all alkaline pH values in- 
vestigated hemoglobin a moves faster than hemoglobin b, and hemoglobin b moves 
faster than hemoglobin c. This shows that the net charge of hemoglobin ¢ is more 
positive than that of hemoglobin b, and the net charge of the latter is more positive 
than that of hemoglobin a, even at a pH of 12.0. Since most lysine and some argi- 
nine groups in proteins have lost their positive charges at a pH of 12.0,5 these ex- 
periments indicate that the differences in charge between the three hemoglobins are 
not due to differences in their content of lysine or arginine residues. 

The paper electrophoretic experiments performed at acid pH values show that 
the three hemoglobins continue to differ in net positive charge down to a pH of 4.8. 
At a pH of 4.1, 3.4, or 3.0, however, the mobilities of the three hemoglobins are 
essentially equal (Table 1 and Figs. 3 and 4). These experiments provide evidence 
that the differences in charge between hemoglobins a, b, and ¢ are due to differences 
In their content of free carboxyl groups. For, in a pH range where only these 
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groups are losing their charges, the differences in charge between the three proteins 
disappear. 

It is necessary to prove that the equality of the mobilities of the three hemo- 
globins at a pH of about 4.1 and lower is not due to irreversible denaturation of the 
proteins. Solutions of the three hemoglobins, brought to a pH of 4.1 with hydro- 
chloric acid, were dialyzed at room temperature against an acetate buffer of pH 4.1 
for 197 minutes. The dialysis sacs were then placed in a veronal buffer of pH 8.6 
and dialyzed overnight. This procedure causes a portion of each hemoglobin to 
precipitate, and the precipitates were removed by centrifugation. A paper electro- 
phoretic experiment was performed with the supernatant solutions in veronal buffer 
of pH 8.6. The result is pictured in Figure 5 and shows that the same differences 
in mobility were observed as had been seen at this pH in experiments in which the 
three hemoglobins had not been exposed to acid (Fig. 1). In a similar experiment 
the precipitated hemoglobins were dissolved in a small amount of sodium hydroxide 
and electrophoresis was performed at a pH of 10.8. The result, seen in Figure 6, shows 
that some of the denatured hemoglobin remained at the starting point but that the 
portions that migrated did so with the usual differences in mobility between hemo- 
globins a, b, and c, which are observed at alkaline pH. The exposure to a pH of 
4.1, therefore, has not irreversibly denatured the hemoglobins with respect to the 
differences in structure which are responsible for the differences in charge at higher 
pH values. Exposure of the hemoglobins for 2-3 hours to a buffer of pH 3.4, 3.0, 
or 11.8, with subsequent dialysis against a veronal buffer of pH 8.6, also yielded 
supernatant hemoglobins still possessing the characteristic differences in mobility 
at pH 8.6. 

Havinga and Itano have cited experiments in which electrophoretic analyses in 
acidic buffers failed to disclose any significant differences in mobility between de- 
natured globins prepared from hemoglobins a and b. They also report experi- 
ments in which the usual differences in mobility characteristic of the two hemo- 
globins were observed in the native globins of these proteins at a pH of 7.0 and in 
a 0.01 molar solution of NagHPO,. They do not report experiments on the de- 
natured globins in alkaline buffers. 

On the basis of these results it would appear unlikely that ionizable groups 
other than carboxy! groups could account for the differences in charge between the 
proteins. Since the imidazole nitrogen of histidine ionizes in a pH range of around 
5.5-8.5,° such a nitrogen atom would be expected to be charged at pH 4.0 and un- 
charged at pH 11. Therefore, a difference in the histidine content of the hemo- 
globins would result in differences in mobility at acid pH values and in no differ- 
ences in mobility in alkaline buffers. This is the opposite of what was observed. 
Differences in tyrosine content should not contribute to differences in charge be- 
tween the hemoglobins, except at pH values near and above the pK value of the 
tyrosine hydroxyl group, which is about 10.8.5 The phosphorus content of hemo- 
globin a is not large enough for differences in the content of phosphate groups be- 
tween hemoglobins a, b, and c to account for their differences in charge.'® Differ- 
ences in specific binding of buffer ions by the hemoglobins cannot account for the 
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Figs. 1-6.—Paper electrophoresis of hemoglobins a, b, and ¢ and dextran, z. The line is the lo- 
cus of the starting points for hemoglobins and dextran. The outlined circles are locations of dex- 
tran spots which stained too faintly to be photographed. The + and — signs indicate on which 
side of the starting line the anode and cathode, respectively, are located. 

Fig. 1, veronal buffer, pH = 8.6, '/2 = 0.06; run for 240 minutes at 400 volts. Fig. 2, phos- 


phate buffer, pH = 11.7, 1/2 = 0.06; run for 126 minutes at 350 volts. Fig. 3, acetate buffer, 
pH = 4.08, ['/2 = 0.06; run for 57 minutes at 400 volts and 140 minutes at 200 volts. Fig. 4, 
glycine buffer, pH = 3.35, '/2 = 0.06; run for 75 minutes at 350 volts. Fig. 5, veronal buffer, 


pH = 8.65, '/2 = 0.06; run for 225 minutes at 350 volts. Hemoglobin solutions, adjusted to 
pH 4.1, were dialyzed against pH 4.08 acetate buffer for 197 minutes, then against pH 8.6 veronal 
pony overnight. Supernatant hemoglobin solutions from veronal dialysis were run in this experi- 
nent. 

Fig. 6, phosphate buffer, pH = 10.8, '/2 = 0.06; run for 95 minutes at 350 volts. Hemoglobin 
solutions, adjusted to pH 4.6, were dialyzed against pH 4.01 acetate buffer for 88 minutes, with 
4 final pH of 4.05 of the hemoglobin solutions. Solutions were then dialyzed against veronal buffer 
of pH 8.6 overnight. The precipitated hemoglobins were centrifuged, and the precipitates, dis- 
solved in a minimum amount of 0.1 N NaOH, were run in this experiment. 
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differences in their charges and mobilities, for two reasons. First, and more im- 
portant, the isoionic points of the three hemoglobins are different in the absence of 
any salt. Second, mobility differences have been observed in several different alka- 
line buffers, and the equal mobilities at low pH have been observed in glycine and 
acetate buffers. Finally, it may be mentioned that ordinary titration studies of 
these hemoglobins!’ have revealed that the total charge per molecule of hemoglobin 
a at various pH values is as follows: pH 3.0, +88; pH 4.0, +47; pH 12, —67. 
These figures indicate that the disappearance of differences in mobility in acid 
buffers is not due simply to. a much greater total charge on the hemoglobin molecule 
in acid than is present in alkali. 

Pauling and his co-workers calculated, on the basis of the difference in isoelectric 
points between hemoglobins a and b, that hemoglobin b possessed two to four more 
net positive charges per molecule than did hemoglobin a.' Using this figure, our 
experiments, therefore, provide evidence that hemoglobin b possesses two to four 
fewer free carboxyl groups than does hemoglobin a. Since the difference between 
the isoionic points of hemoglobins b and c is of the same order of magnitude as, 
though somewhat larger than, the difference between the isoionic points of hemo- 
globins a and b, it seems probable that hemoglobin ¢ contains about five to eight 
fewer free carboxy! groups than hemoglobin a. In order to calculate quantitatively 
the differences in content of free carboxyl groups in hemoglobins a, b, and ¢, we 
have attempted to combine paper electrophoretic experiments, in which electro- 
osmotic flow is measured, with acid-base titration experiments. These results will 
be reported in a subsequent paper. 

The structural basis for the apparent difference in the content of free carboxy! 
groups of these hemoglobins remains obscure. The available analytical data® 
do not indicate a smaller content of aspartic or glutamic acid residues, or a larger 
number of carboxyl groups which exist as amides, in hemoglobin b than in hemo- 
globin a. No amino acid analyses of hemoglobin ¢ have been published. 

It has been suggested that hemoglobins a and b have identical polypeptide chains 
but that the folding of the chains is different in the two molecules, presumably re- 
sulting in some charged groups being unreactive.» 8 An actual masking of acid- 
binding groups, perhaps due to folding of the molecule, has recently been experi- 
mentally demonstrated in native horse hemoglobin by Steinhardt and Zaiser.® ' 
These workers have shown that thirty-six acid-binding groups are not available for 
titration until the protein has been exposed briefly to a pH of about 4.0 or less. 
When these thirty-six groups are broken out, the hemoglobin is simultaneously 
denatured, as judged by spectroscopic changes and loss of solubility at the isoelectric 
pH. Renaturation, and presumably remasking of the thirty-six groups, is effected 
by returning the protein to neutrality. Steinhardt and Zaiser consider these 
thirty-six groups to be epsilon amino or guanidinium groups, since their pK values 
are above 5, presumably too high for carboxyl groups, and since all the histidine 
residues of horse hemoglobin are titrated in the native protein.’ It is tempting to 
consider whether the excess positive charge on hemoglobins b and ¢ could be due 
to an analogous masking of differing numbers of charged groups in the three hemo- 
globins. The fact that the variations in charge persist at pH 12.0 indicates that 
these differences are not due to epsilon amino groups and are not very likely to be 
due to guanidinium groups. There is the possibility that carboxyl groups may 
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also be masked in hemoglobins with more of them blocked in the abnormal hemo- 
globins. Other more complicated masking possibilities may be imagined to ac- 
count for the differences in charge, but the evidence provided by the experiments 
reported above makes it very probable that the three hemoglobins differ in their 
content of free carboxy! groups. 

The method of differential titration described in this paper should be applicable 
to the study of other sets of closely related proteins. 


* This work was supported in part by a grant (A-572) from the National Institute of Arthritis 
and Metabolic Diseases of the National Institutes of Health, United States Public Health Service. 

It is a pleasure to thank Professors John T. Edsall and George Scatchard for their very helpful 
critical discussions. We are also indebted to Dr. Helen Ranney, Miss Grace Vanderhoff, and 
Mrs. Catherine Holavko for their assistance in obtaining the hemoglobins studied and in con- 
structing the paper electrophoresis apparatus. 
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ON THE REPLICATION OF DESOXYRIBONUCLEIC ACID (DNA) 
By M. DetsritcKk 
KERCKHOFF LABORATORIES OF BIOLOGY, CALIFORNIA INSTITUTE OF TECHNOLOGY 
Communicated May 18, 1954 
The discoveries of Hershey and Chase! concerning the role of DNA in transmitting 
genetic information in phage and of Watson and Crick? concerning the structure of 
DNA have brought the problem of the replication of DNA into focus. The struc- 


ture proposed by Watson and Crick consists of two polynucleotide chains wound 
helically around a common axis, tied together by hydrogen bonds between the 
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purine and pyrimidine side chains. These side chains of the two chains are ar- 
ranged so that adenine is always matched with thymine and guanine with cytosine, 
The sequence of bases along either chain is not subject to any restrictions, but once 
the sequence along one chain is given, the sequence along the other chain is com- 
pletely determined. This sequence, then, constitutes the genetic information, a 
linear message written in a four-symbol code. The duplex of the two chains con- 
tains the information in a twofold redundance. Each chain has a directional 
polarity because of the nonequivalence of the 3- and 5-positions through which each 
pentose is linked to the preceding and the following phosphate group in the chain. 
This polarity runs in opposite directions in the two chains of the duplex. 

Watson and Crick* have proposed that this structure is replicated by a process in 
which the chains of a duplex are separated and each catalyzes the synthesis of a 
complementary chain. This is schematically represented in Figure 1. One 
imagines that the synthesis of the complementary chains occurs in a zipper-like 
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fashion along both the parental chains, so that, when synthesis has progressed to 
point n, the niche formed by link n of the new chain and link n + 1 of the old chain 
provides a suitable enzymatic surface for the insertion of the correct nucleotide as 
link n + 1 of the new chain and for the rejection of the three incorrect nucleotides. 

The principal difficulty of this mechanism lies in the fact that the two chains are 
wound around each other in a large number of turns and that, therefore, the daughter- 
duplexes generated by the process just outlined are wound around each other with 
an equally large number of turns. There are three ways of separating the daughter- 
duplexes: (a) by slipping them past each other longitudinally; (6) by unwinding 
the two duplexes from each other; (c) by breaks and reunions. 

We reject the first two possibilities as too inelegant to be efficient and propose to 
analyze the third possibility. 

If one tries to separate the two chains of a duplex by moving the two chains 
laterally in opposite directions, an interlock occurs for each turn of the helix, i.e., at 
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each tenth link. Such an interlock can be resolved in two ways: (a) by breaking 
one of the chains, slipping the other chain through the gap, and rejoining the 
broken ends (Fig. 2, a); (6) by breaking both chains at each half-turn and rejoining 
them crisscross (Fig. 2, b). 

We reject both these mechanisms—the first one because it introduces an asym- 
metry between the two chains (only one of them being broken) which is contrary to 
the symmetry of the structure and the second one because it rejoins chains with 
opposite polarity, which is chemically not permissible. We conclude that it is not 
feasible to separate by breaks and rejoins the two chains of a single duplex. The 
situation is quite different, however, when one considers a duplex during the proc- 
ess of replication. Let us consider a duplex in which replication has proceeded 
synchronously along the two chains up to link n. We will call this point the 
“growth point.” If we now break both the old chains between links n and n + 1, 
we may join the lower terminals of the breaks in a crisscross fashion, not to the 
upper terminals of the breaks but to the open ends of the new chains of equal 
polarity. The upper terminals of the breaks now become the open ends for the 
continuation of the replication process. This process is illustrated in Figures 3 
and 4. In Figure 3 we illustrate the breaks and rejoins in lateral view. Figure 4 
illustrates a series of successive cross-sections through the replicating duplex, viewed 
from above, one gross-section for each successive link, in the vicinity of the growth 
point. The lower cross-sections are below the growth point and illustrate the rota- 
tion of the two chains around the common axis by an angle of 36° for each link, in a 
counterclockwise direction as we proceed upward (right-handed helix). The upper 
cross-sections are above the growth point and illustrate the two daughter-duplexes, 
and in each of them they indicate again the rotation of the chains of the daughter- 
duplexes around the axis of each duplex. As the growth point moves down one 
notch, the two members of the next pair of complementary links of the parental 
duplex move laterally, to the left and right, respectively, without change in ortenta- 
tion, and to each of them a new complementary link is added as the next link of the 
new chains of the daughter-duplexes. As this process continues from link to link, 
obstructions arise at every half-turn of the helix, i.e., at every fifth link. These 
obstructions can be resolved by making at every half-turn breaks and rejoins indi- 
cated in Figure 3 in lateral view and in Figure 5 in cross-sectional view. 

Chemically these breaks and rejoins amount to exchanges of one nucleotide bond 
for another, i.e., to transnucleotidations. No energy is consumed or liberated in the 
process. The activation energy needed is probably similar to the activation energy 
involved in the addition of a new link, since the pool of DNA precursors probably 
exists not as free nucleotides but in such a form that the functional groups to be 
linked are tied up by chemical groups which are exchanged for the new partner. 

The daughter-duplexes generated by this process are correctly coiled from the 
start; in fact, there does not occur at any time any twisting of the parental chains. 

We consider it an attractive feature of this proposed mechanism that the comple- 
ment formation proceeds synchronously down the two chains. This means that at 
any one time only one short section of a duplex is the “working section.”’ Only in 
this section will the structure be different from that of a stable duplex. This is 
attractive for two reasons: (1) the energy required at any one time for the separa- 
tion of the parental duplex is a minimum; (2) since the total length of a duplex, 
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when fully extended, is very large compared to the intracellular space available, the 
actual length of the duplex must be enormously reduced by some process of super- 


coiling (or, less likely, folding). The working section, 
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and reunions the distribution of label would occur only at the first replication, 
At each subsequent replication one daughter-duplex would receive all of the label, 
the other none (Fig. 6, b). At present it does not seem possible to discuss the bear- 
ing of this implication on the experiments of Stent* concerning the mortality, due 
to the decay of incorporated P* of phage infective centers at various stages of the 
reproductive cycle. These moftality experiments are complicated by the phenom- 
enon of multiplicity reactivation, i.e., an interaction between different duplexes, 
the nature of which is still uncertain. 

Another implication of the proposed mechanism is the fact that after each odd- 
numbered exchange the open ends of the chains are those of parental material. If 
the addition of new links is a reversible process, then these open ends of parental 
nucleotides offer the possibility of exchange of parental labeled nucleotides with 
unlabeled nucleotides from the pool. One may estimate that such exchange reae- 
tions could lead to a maximum loss of label of 10 per cent per replication cycle from 
the aggregate of vegetative phage particles. Since there are 5-10 replication cycles 
during one intrabacterial growth cycle, such losses could easily account for the ob- 
served low values (80-50 per cent) of transmission of parental P** to the progeny,’ 

Another implication concerns hypothetical ring duplexes, analogous to ring 
chromosomes. In the first place, such a ring duplex contains necessarily a whole 
number of turns, since a half-integer number of turns would necessitate joining 
chains of unequal polarity. In the second place, the daughter-duplexes are separate 
rings if there are two pairs of breaks and rejoins for exactly every turn. How likely 
such an exact coincidence would be for large rings depends on details of strain rela- 
tionships, about which we know too little at present. In any event, exact equality 
should be more probable than any specified inequality, and its absolute probability 
should be greater the smaller the ring. 

Summary.—Watson and Crick have proposed a mechanism for the replication of 
DNA. This mechanism involves the synthesis of complementary polynucleotide 
chains on each of the complementary chains of the parental structure. In this 
paper a mechanism is proposed for resolving the interlocks which prevent the 
separation of the daughter-duplexes. 

It is assumed that (1) the complementary synthesis proceeds synchronously along 
the two chains, and (2), as synthesis proceeds, the chains break at the growth point 
at every half-turn of the helix (every fifth link); the lower terminals of the breaks 
are immediately rejoined to the open ends of equal polarity of the new chains. 
Some implications of this mechanism are pointed out. 


1A, D. Hershey, and M. Chase, J. Gen. Physiol., 36, 39, 1952. 
2 J. D. Watson and F. H. C. Crick, Nature, 171, 737, 1953. 
3.J. D. Watson and F. H. C. Crick, Nature, 171, 964, 1953; Cold Spring Harbor Symposia Quant. 

Biol., 18, 123, 1953. 

4G. 8. Stent, Cold Spring Harbor Symposia Quant. Biol., 18, 255, 1953. j 

5T,. M. Kozloff, Cold Spring Harbor Symposia Quant. Biol., 18, 209, 1953. 
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THE VAPOR PRESSURE OF SMALL DROPS* 
By W. H. Roprsusut 


CLOUD PHYSICS LABORATORY, UNIVERSITY OF CHICAGO 


Communicated June 23, 1954 


Gibbs did not consider the Kelvin equation for the vapor pressure of small drops 
for the obvious reason that it does not constitute chemical equilibrium. The 
Kelvin equation represents true equilibrium between vapor and a flat surface of 
liquid water, where the liquid is under a pressure greater than the vapor by any 
assigned value AP. For example, in a closed vessel in the presence of an inert, 
slightly soluble gas such as air a definite equilibrium exists which will be re-estab- 
lished if displaced by the transfer of water from the liquid phase to the vapor or 
vice versa. 

For equilibrium between vapor and small drops as described by the Kelvin 
equation this is no longer the case. For a given partial pressure of vapor a drop 
must be of a particular radius. A drop which is slightly larger will grow by con- 
densation, and a drop which is slightly smaller will evaporate. The situation 
does not even correspond to unstable equilibrium, since it does not satisfy the 
condition AF = 0 so far as free surface energy is concerned. Gibbs! pointed out 
that the production of a drop of the required size from the vapor in the Kelvin 
equilibrium required an increase in free energy equal to one-third the surface free 
energy. The vapor in quasi-equilibrium with the drop has an additional free 
energy per mole 2 yV°/r, where y is the surface tension and r is the drop radius. 
But the free surface energy of a mole of water in the form of drops is 3 yV°/r, as 
will be shown later. Incidentally, free surface energy can only be of significance 
when the equilibrium involves the creation or disappearance of surface. This 
only happens when the liquid phase is subdivided into small drops. 

The equilibrium between a vapor and small drops is not a simple affair. Pro- 
vided that the drops do not exceed a certain size, such an equilibrium is possible; 
but it necessarily involves a whole gamut of drop sizes and concentrations. For- 
tunately, we need consider only one particular drop size at a time in discussing the 
equilibrium. We shall proceed to discuss the range of conditions under which 
such equilibrium can exist. Before actually formulating the equilibria involved, 
it is necessary to consider the special thermodynamic properties which are possessed 
by small drops. 

Surface Energy.—The first of these is the free surface energy, which, expressed 
in ergs per square centimeter, is the surface tension. This is ordinarily assumed 
to be independent of drop radius, but this cannot be true down to the point where 
the radius becomes that of a single molecule and the drop vanishes. This problem 
has been considered first by Gibbs and later by Tolman? and others, but no definite 
conclusions can be reached. What we are concerned with is not the validity 
of surface tension as a physical concept but rather the decrease in the free surface 
energy of the drop as the radius is decreased. From thermodynamics, 
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where o is the surface area in square centimeters. The first term on the right is 
larger than 7(0S/0c)7, p and will presumably always be larger, so that we can 
confine our attention to it. Stefan, following Laplace’s theory of surface tension, 
concluded that, if the center of an atom or molecule was in the surface, the particle 
had lost one-half its neighbors and, therefore, that the enthalpy change in bringing 
a particle to the surface would be one-half the enthalpy of vaporization. This 
seemed to be a naive point of view, but the values for mercury are in fair agreement 
with Stefan’s theory. In the case of associated liquids, the surface enthalpy is 
very much smaller than one-half the enthalpy of vaporization, which leads to the 
conclusion that the molecules are able to maintain hydrogen bonding even in a 
curved surface, and there seems to be no reason to assume an appreciable increase 
in surface enthalpy so long as the drop contains one hundred or more molecules. 
Since the co-ordination number in the liquid is probably five at a maximum, 
this will not be seriously constrained until the radius becomes very small. 

While the second term of the right-hand side of equation (1) must always be 
exceeded by the enthalpy term, there is another contribution to the entropy of 
small drops to be taken into account, viz., the entropy of the drop considered as a 
gas molecule. There is no corresponding enthalpy term to be considered, and this 
term must become the overriding factor leading to stabilization as the drops ap- 
proach molecular dimensions. 

Effect of Drop Concentration.—Equilibrium in gaseous systems depends upon 
the partial pressures or concentration of the various species present. This de- 
pendence upon concentration comes in through the dependence, mentioned above, 
of the entropy of the drop, considered as a gas molecule, upon the concentration. 
Liquid water may be considered as having the classical number of degrees of free- 
dom per molecule which contribute to its enthalpy and entropy. In the liquid 
drop, however, considered as a gas molecule, three of these degrees of freedom 
must be assigned to translation and three to rotation of the drop as a whole, and 
these contribute more to the entropy than does the average degree of freedom. 
The contribution of the free surface energy to the free energy per mole of the 
liquid drop varies inversely as the radius of the drop and becomes negligible in 
1 u. The contribution of the entropy of the drop qua gas molecule diminishes as 
the inverse cube of the radius, and, while this contribution depends upon the 
concentration, all these terms become negligible at 1 uw radius. For practical 
purposes, therefore, drops of 1 » radius may be considered as a continuous liquid 
phase with a normal vapor pressure. In addition to the contribution to the 
absolute entropy, the translational degrees of freedom involve the concentration 
of the drops through the term —R In P,,, where P,, is the partial pressure of the 
drops considered as gas molecules. This comes about because of the Bose-Einstein 
statistics which must be applied. 

Spontaneous Nucleation.—As a result of these considerations, we may postulate 
a given concentration for any of the possible drop sizes in equilibrium with the vapor 
at any temperature and pressure. So long as the vapor is unsaturated, the con- 
centration of the larger drop sizes is vanishingly small; but it is enormously sensi- 
tive to changes in vapor pressure, and at some point of supersaturation the drop 
concentration becomes finite at all radii. Another way of expressing this phenom- 
enon is as a rate process. At vanishingly low concentrations the time required 
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for growth to finite (i.e., observable) size is infinite; a slight increase in super- 
saturation will shorten this time to a matter of seconds. 

Kinetics of the Condensation Process.—Since all kinetic processes in gases at 
moderate pressures must be assumed to be bimolecular, we may consider the growth 
of the drop or cluster as a molecule by molecule condensation: 


Fr, +1 
. ee - _ eee a Kp. 
P, 
For purposes of calculation, however, it is convenient to consider the equilibrium 
. . p" ‘ 
We ne. P = ie, (2) 
n 
If we write the equilibrium 
+1 
, > p" , sit 
Xa = (n + 1)X, 5 = Kan, (3) 
Patt 


then we see that these three types of equilibria are consistent and derivable from 
each other. Subtracting equations (3) from equations (2), we have 


b + X = ». ae ci Kp. (4) 


Equilibrium between Cluster and Vapor.—Using the term “cluster” to designate 
the drop of diminishing size, let us write the equation for the equilibrium between 
the cluster of n drops and vapor, 

| p 


Xn — xX, 
n |: nis 


= K. (5) 


This equation, as written, is on the basis of one mole (18 gm.) of water, and the 
corresponding free energies will be so calculated. Using the prime for the cluster 
and the double prime for the vapor, we have for the chemical potentials 


p’ = F,° + RT in P,, 
uw’ = F,° + RT In p. 


For the gaseous equilibrium 


and 
3 yV° I 
BP? a Be + +— Tg". (6) 
r n 
Here V° is 18 cm.*, and 3 yV°/r represents the free surface energy (ergs) of 18 cm.* 
of liquid water in the form of drops of radius 7, S’ ° is the excess entropy of 1 mole 
of clusters considered as gas over that of liquid water, and n is the number of 
molecules per cluster. 
If we remember now that In po”, the saturation pressure of water, is given by the 
equation 
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l —aFo Fy? — F,° 
n 1 = Ese = — oe ’ ” 

R1 RI ) 

where F’,° is the free energy of water vapor at the temperature 7' and 1 atmosphere, 

assuming ideal-gas behavior, and F,° is the free energy per mole of liquid water 

likewise under standard conditions, we may write, from equations (5), 
" ‘ ae) 

p . p : : 3 I \ Pee 

—RT In K = —RT In = Fy° — F,° — + -—TS” 

j t 


” 


Po 7 7 


or 


p" 3yV° S—RInP, : 
. po”  rRT nk (8) 
Comparison with Experiment.—According to data in the current literature, air 
saturated with water vapor at 25° C. will, on cooling to —10° C. by adiabatic 
expansion, show evidence of spontaneous nucleation. We shall accordingly cal- 
culate the terms in equation (8) for these conditions. For the surface-energy term 
we shall assume 7 = 263°, y = 75 ergs/cm, and the volume of the water molecule 
30 X 10-** cm. 
The surface-energy term becomes 


For the entropy term we have the following expressions from Paul :* 
Si, = 26.0 + (3/2)R In M, 
where M is the molecular weight of the drop. (= 18n). 


Srot = 284.6 + (3/2)R In if 
I = (2/5)mr? (m = mass of drop). 


Assuming a symmetric number of 1, radius of 10 A, and 133 molecules per drop, 
one calculates 


Sr = 49.0, Sroc = 36.0. 


It remains to calculate the value of R In P,,, which must be expressed in atmospheres. 
A concentration of 107 drops per cubic centimeter would correspond to P,, = 10~” 
atmospheres; hence 


Rin 10-” = 55.9 E.U. 


107 nuclei should certainly produce a visible fog. A smaller concentration might 
suffice. 

These assumptions give, for the total entropy term, 140 E.U. These values, 
therefore, lead to the equation 
Pp _ 96 _ 70. pe 
Po Vn on 
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This equation gives a maximum at n = 102, and tables of calculated values are 
given herewith (Tables 1 and 2). The data are also plotted in Figure 1. The 
value In (p/po) = 1.37 at the maximum corresponds approximately to the fourfold 
degree of supersaturation that has been frequently reported in the literature. 


TABLE 1 TABLE 2 
VALUES OF P, FOR 
CALCULATIONS OF SUPERSATURATION FROM HQUATION In (p/po) = 1.37 
sy 3s Pn 
n Vn 9.6/V n 70/n In (p/ po) n (Molecules/Cm ) 
64 4 2.4 ie 1.30 64 > 10° 
102 4.68 2.06 0.69 1.37 102 107 
125 5 1.92 0.56 1.36 125 > 108 
216 6 1.60 0.32 1.28 


It is believed that the surface-energy term cannot be in serious error so long as 
the critical drop size contains more than one hundred molecules. The entropy 
term may be in error, but it is almost certainly too small. While an increase of 
several units would increase the critical drop size calculated, the effect is not very 
great because n is a divisor in the term. The maximum as calculated in Table | 
is very flat, but this is somewhat misleading. If one solves equation (8) for In P,, 
using the value for the supersaturation In (p/po) = 1.37, one obtains a sharp 
minimum at n = 102. Thus = 102 is the critical drop size as defined for spon- 
taneous nucleation. Drops of this size exist at smaller concentration at equilibrium 
than drops of larger or smaller size. The equilibrium concentration is enormously 
sensitive to the degree of supersaturation. 








4 8 12 16 20 24 
rA 
Fic. 1.—Logarithm of supersaturation plotted against radius in angstroms. 


No attempt has been made in this paper to consider the rate of growth of drops 
or to give a critical review of data as reported in the literature. These questions 
will be considered in another paper which will give an account of some experimental 
attempts to produce spontaneous nucleation. 

It should be obvious that the thermodynamic method used here can be applied 
to spontaneous nucleation in any sort of phase transition—for example, crystal- 
lization from a supersaturated solution. The entropy term is always manageable, 
but one must seek some sort of experimental approach to the estimation of inter- 
facial energy. 
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COMMENTS ON A “MESO-CARBON ATOM” 
By M. L. Woirrom 
DEPARTMENT OF CHEMISTRY, THE OHIO STATE UNIVERSITY, COLUMBUS, OHIO 
Communicated July 12, 1954 


Schwartz and Carter! report a disproportionate yield (60-40 per cent) of two 
products when (/evo)-a-methylbenzylamine was reacted with 3-phenylglutaric an- 
hydride, and claim that this demonstrates the stereochemical nonequivalence 
of the two groups attached to the CsH;(H)C< system of the 3-phenylglutaric 
anhydride. They designate such a substituted carbon a “meso-carbon atom.” 
Actually, Schwartz and Carter have subjected a compound containing an optically 
active center to an equilibrium reaction in which a new asymmetric center is formed. 
The products are then diastereoisomers, possessing different chemical energies, 
and the reverse reaction will set up an equilibrium involving a disproportionate 
ratio of products, as shown below (the absolute configurations are assumed). 


CeHs CO—CH, CeHs 


ithe T ‘ . iN, Y ‘ ‘ 
Hi—-SNH, + O HC—C,;H; = Ha-—-—NH—CO—CH, 
SZ te = Se 
CH; CO—CH, CH; ri% : 
HM——>C,H; 


CH.—COH 


4h ae ; ; 
HM@—S—NH—CO—CH, 

NG 
ve CH SH 
CH.-—C0.H 


The above type of reaction is well known and is no different in principle from 
the Kiliani reaction, whereby p-arabinose is converted into disproportionate 
amounts of p-glucose and p-mannose. The cause of this disproportionation resides 
in the optically active nature of the initial compound and not in any “stereocheml- 
cal nonequivalence” of two groups in the optically inactive substance, real as that 
may be. In the early days of the development of the Van’t Hoff—Le Bel theory 
of the tetrahedral carbon atom, many experimental attempts were made to resolve 
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meso compounds. A sufficient number of negative instances were established to 
make such an objective unfruitful toward further pursuit. The term meso had 
best still be employed as a prefix to designate an optically inactive and unresolvable 
stereoisomer of one or more optically active substances. 

Schwartz and Carter! also mention an “‘axis of symmetry.”’ ‘This is an ambiguous 
term, since there are two types recognized by the crystallographer: an alternating 
axis of symmetry and a rotating axis of symmetry. The absence of the former 
in the molecular model is a requisite for optical activity, but the absence of the 
latter is not, and indeed such long-known substances as the optically active forms 
of tartaric acid possess such symmetry. 


! Pearl Schwartz and H. E. Carter, these PRocEEDINGs, 40, 499, 1954. 


RADIATION-INDUCED EXCHANGES IN DROSOPHILA FEMALES 
By D. R. PARKER* 
BIOLOGY DIVISION, OAK RIDGE NATIONAL LABORATORY, OAK RIDGE, TENNESSEE 


Communicated by J. T. Patterson, July 14, 1954 


The purpose of the present experiments is to determine the nature of the ex- 
changes involved in detachment of attached-X chromosomes, with particular em- 
phasis on the possible involvement of the Y chromosome. Rapoport! reported 
that the rate of X-ray-induced breakdown of attached-X chromosomes in Dro- 
sophila melanogaster is not affected by the presence or absence of a Y chromosome; 
he concluded from this observation that the process of induced detachment did not 
involve exchange with the Y, differing in this respect from spontaneous detach- 
ment.” * Rapoport’s failure to find evidence of long, 2-break deficiencies or of 
translocations with autosomes led him to believe that a process of simple fragmen- 
tation, with “restitution of the telomere,’ is the mechanism of induced detach- 
ment. However, he failed to obtain any evidence of such fragmentation of chro- 
mosome 3 in 2,411 cultures derived from females given 3,000 r of X-rays. This 
leaves open the question of the nature of the induced separation of attached X’s. 

It is necessary to know the stage of odgenesis in which detachment may most 
readily be induced, in order to determine how the samples are to be taken in experi- 
ments dealing with the frequency of detachment. The writer®> has shown by use of 
egg counts that detachments are readily induced in late odcytes. Attached-X 
females of the constitution y v f car, y v f, or br ec were treated (1,000 r of X-rays, 
50 kvp., 10 ma., 100 r/min) 0-12 hours after emergence and each one mated imme- 
diately to two Muller-5 (sc>! B InS w* sc’) males. Each day for a total of 7 days 
the adults were transferred to new food, and egg counts were made. Over this 
period the females laid an average of 140 eggs. The females were allowed to ovi- 
posit for 3 additional days, but no egg counts were made. Taking the minimum 
number of eggs laid as twice the number of adult offspring recovered during this 
3-day period, the average number of eggs laid per treated female for the entire 10- 
day period was at least 210. The results are shown in Figure 1; each line represents 
a sample of eggs in which a detachment was recovered. Of the twenty-eight re- 











796 GENETICS: D. Rk. PARKER Proc. N, A. §, 


covered in this experiment, from twenty-three to twenty-five were found among 
the first ninety eggs recovered from each of the various females whose progeny 
included cases of detachment. Thus most detachments must be induced in late 
oocytes. 

The viability of individuals with detached-X chromosomes likewise might be ex- 
pected to affect the observed frequency of detachment. It was found that with the 
attached-X y v f car stock the flies recovered with a detached-X chromosome were 
invariably females, suggesting either a high rate of induction of lethals or the pre- 
existence in the attached X of spontaneous recessive lethals. 

The latter hypothesis was tested by using various attached-X stocks. Table | 
gives the sex of flies in which detached-X chromosomes were recovered, together 
with the frequency of lethals in detached X’s recovered in females. The wide vari- 
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ACCUMULATED TOTAL NO. EGGS LAID BY INDIVIDUAL TREATED FEMALES 


Fic. 1.—Recovery of detached-X chromosomes. Each line indicates sample of eggs in which a 
detachment was recovered. 


ation in these values among different stocks is not consistent with an interpreta- 
tion that the treatment induces lethals and sublethals; rather, it supports the hy- 
pothesis that spontaneous recessive lethals might accumulate in some attached-X 
stocks. 

These spontaneous lethals probably accumulate near the centromere, where 
homozygosis frequency is low and selection relatively ineffective. It might be 
possible occasionally to obtain a stable system of balanced lethals through the origin 
of closely linked lethals in each of the homologues. Perhaps the Austin y v f car 
stock is an example of such a balanced lethal attached X. 

The extent to which the Y chromosome is involved in the exchanges leading to 
these detachments was determined by comparing the frequency of detachment and 
the type of exchange involved in the absence of a Y chromosome and with various 
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modified Y’s. The treatment employed a G.E. Maxitron-250 X-ray machine oper- 
ated at 80 kvp., 30 ma., with a 3-mm. aluminum filter. The dose rate was approxi- 
mately 100 r/min as measured with a Victoreen meter. Irradiated females of the 
attached-X stock y? su-w* w* bb were placed individually in vials with two “‘Muller- 
5” (sc?! B InS w* sc’) males each and allowed to oviposit for 4 days, after which 
they were removed. Since the viability of Muller-5 males varies with the type of 
Y chromosomes present, X-0 males being almost completely inviable, and since the 
F,; attached-X females would in all cases have a normal Y derived from their Muller- 
5 fathers, the rates of detachment are based on the numbers of surviving females of 
y’ su-w’ w° phenotype. Frequencies, as given in Table 2, are based on a total num- 
ber of zygotes equal to four times the number of such females, plus the number of 
individuals with detachments. 


TABLE 1 


Sex or Fires 1N Wuicu DetTAcHED-X CHROMOSOME Was RECOVERED AND FREQUENCY OF LE- 
THALS IN DETACHED X’s RECOVERED IN FEMALES 


ATTACHED-X NUMBER OF DETACHMENTS RECOVERED IN FEMALES 
Stock/Y CHROMOSOME DETACHMENTS RECOVERED No. Tested No. Per Cent 
(1p ANY) Females Males for Lethals Lethals Lethals 
yo f car (Austin) /Y 16 0 13 13 100 
yvcar/Y 12 2 10 8 80 
yuf car (Oak Ridge)/sc?- Y 25 7 8 6 75 
y? su-w* w* bb (no Y) 94 87 90 6 6.7 
y? su-w? w* bb/sc8» Y oa 35 36 l 2.8 
y? su-w? w* bb/scY!+ YS 149 100 126 5 4.0 
y? su-w* w* bb/se+ Y" 50 50 41 3 to 
y? su-w* wt bb/ YL" 10 10 10 5 50 
y/Y 12 11 10 1 10 
br ec/Y 17 17 15 l 6.7 


TABLE 2 


FREQUENCY OF X-RAY-INDUCED DETACHMENT AND FREQUENCY OF X-4 LINKAGE AMONG THE DE- 
TACHMENTS, Usinc Various Y CHROMOSOMES IN COMBINATION WITH y/? su-w* w? bb ATTACHED X 








———Doss: 1,000r - Dose: 2,000r ; 
Adjusted Adjusted 
No. Per- No. Per- PERCENTAGE OF 
Cases of No. Fi, centage Cases of No. Fi centage DETACHMENTS 
: Detach- yy? su-wa@ wa Detach- Detach- y? su-wa wa Detach- GIVING 
Y CuromosomE ment Females ment ment Females ment X-4 LinKAGE* 
None 44 5,454 0.20 66 3314 0.50 59.8 (61/102) 
sc. } 8 813 0.25 30 870 0.85 13.2 (9/68)T 
see YE ie OSs nee 97 2839 0.85 13.0(10/77) 
scl! YS site er ais 247 3927 1.55 8.9 (8/90) 
Ni se ve i 20 944 0.53 80.0(12/15) 


* This includes only one of the two classes of products of translocation between attached X’s and chromosome 
4; the complementary class cannot be detected by linkage of X and 4. ) . 
Bt... dicate stocks tested include some derived from mass matings, which are not included in other columns 
of this table. 


Attached-X females lacking a Y chromosome were obtained by using males with 
the attached X-Y chromosome with no free Y. A comparison of detachment rate 
in such females with that in females with the sc*- Y chromosome shows the rate of 
detachment to be somewhat increased by the presence of a Y chromosome (Table 2). 

Detachments obtained from females lacking a Y were tested for linkage with chro- 
mosome 4; 61 of the 102 cases tested were found to have chromosome 4 linked to 
X, with the break in the fourth being proximal to the cz locus in all but one case. 
The complementary product of exchange between X and 4, which should occur 
with equal frequency, is an X with centromere of the fourth but cannot be identi- 
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fied as such by linkage of X and 4. This raises the question whether detachments 
associated with all possible segregation products are equally viable. Thirty-four 
of the detachment stocks showing linkage between X and fourth chromosomes 
were found at the time of testing to carry an extra fourth chromosome; the com- 
plementary class from this type of segregation would be haplo-4 individuals with 
X attached to centromere of the fourth. The failure to find any such haplo-4 
individuals among the fertile exceptional flies indicates that not all detachments 
are being recovered. On this basis, about 14 of an estimated minimum of 136 
detachments must be the result of some type of exchange other than between X and 
t. The lower frequency of detachment in females lacking a Y as compared to those 
with a sc’- Y may in part be due to loss of certain exchange products through cer- 
tain types of segregation, including occasional nondisjunction of the separated 
X’s as a result of X and fourth centromeres going to the same pole. There should 
be no similar decreased recovery of detachments due to aneuploidy or nondisjunc- 
tion in the case of detachments involving the Y chromosome. Inequality of rate of 
recovery of complementary exchange products based on chromosome length might 
be expected,* ® although the expected higher frequency of recovery of shorter ex- 
change products was not found in the case of the tandem acrocentric compound X.’ 

The sc’- Y chromosome is shaped approximately like a normal Y, with a small 
piece of the X including the normal allele of yellow (y*) near or at the tip of the long 
arm, presumably. A total of sixty-eight detachments obtained from females with 
this marked Y were tested for segregation of y+ and for fertility of males when either 
the short arm (sc’'-Y*%) or the long arm (sc- Y“) of the Y is added. The stocks 
tested include some recovered from mass matings of treated females, as well as those 
reported in Table 2. Thirty had y+ linked with the X, including nine which were 
fertile with sc’!- Y* and therefore were X-Y". Thirty-eight showed the phenotype 
yellow?; of these, nineteen were found by fertility tests to be X-Y*, since males 
were fertile with the addition of sc- Y but not with sc”!-Y°. Four were found to 
have fertility factors of both arms of the Y linked with the X and are designated 
X-Y"Y°*; these are fertile with either arm of the Y added and are fertile with no 
free Y. Nine cases were found to have X and chromosome 4 linked, while six 
have not been identified as to type of exchange. Thus approximately 80 per cent 
of the observed detachments arise from exchanges with the Y. Considering the 
possible differences in probability of recovery of exchanges of the different types, 
there seems to be a slight increase over the no-Y females in total number of ex- 
changes and a marked reduction in the frequency of exchanges between X and 4. 

With a 2-armed Y chromosome, sc: Y“, lacking fertility factors of the short arm 
but having the bb+ locus, there is no apparent change either in detachment rate or 
in the rate of X-4 translocations as compared with the sc’- Y. Ten of seventy-seven 
cases showed X-4 linkage. On the other hand, the l-armed short arm of the Y, 
sc’'.Y°, allowed a marked increase in detachment rate, with apparently no change 
in the frequency of X-4 translocations, eight out of ninety tested showing X-4 link- 
age. The increase in exchange frequency between X and Y° thus is not at the ex- 
pense of exchanges between X and 4. 

Use of the closed long arm of the Y, Y““, should allow the recovery only of de- 
tachments which did not involve exchanges with the Y chromosome. It was found 
that the detachment frequency did not differ significantly from that in attached-X 
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females lacking a Y. ‘The exchanges were found to be predominantly X-4 trans- 
locations, twelve out of fifteen tested showing X-4 linkage. 

Discussion.—Detachment of attached X’s regardless of presence or absence of a 
Y chromosome is a 2-hit process,® * * varying nearly as the square of the dose at 
low dosages. The occurrence of a high frequency of X-4 translocations in oécytes 
in which a Y chromosome is not present indicates close association of these chromo- 
somes during meiosis, perhaps due to homology. Exchanges of this type cannot be 
detected by ordinary genetic methods in free-X females, since the usual results 
would be exchange of centromeres and of heterochromatin. It would be necessary 
to mark the centromeres either of X or of 4 to detect this type of exchange. 

If the reduction in frequency of exchanges between X and 4 when the Y chromo- 
some is present is to be explained on the basis of random competitive rejoining of 
broken ends of chromosomes, the detachments which are the result of exchanges 
between X and Y should largely be in addition to, rather than in the place of, those 
due to exchanges between X and 4 (as well as being in addition to long deficiencies 
which have been reported recently*). A great reduction in X-4 exchange frequency 
might be achieved through very high frequency of breakage of the Y, resulting in a 
very great increase in the frequency of detachment. Failure to find this type of 
result does not completely invalidate the hypothesis, since one might expect an ap- 
preciable number of exchanges between Y and 4 to be occurring at the expense of 
X-4 exchanges. Information on Y-4 translocations in this system must be obtained 
before any decision can be reached concerning the alteration in X-4 exchange fre- 
quency. 

An alternative hypothesis is offered that the presence of a Y chromosome in some 
way “protects” chromosome 4. This “protection” might be due to pairing rela- 
tionships operating in such a way that X and 4 are more widely separated when a 
Y is present, hence less likely to be involved in mutual exchanges. 

The failure of the closed long arm of the Y, Y““, to bring about a reduction in 
exchanges between X and 4 may be of significance. On the basis of the competi- 
tive rejoining hypothesis, there must be no breaks in the Y“° which can combine 
with breaks in the X. Since all such exchanges presumably would be eliminated 
as dicentrics, the rate of recovery of detachments should be very low. On the 
other hand, if there is ‘protection’ of chromosome 4 by the Y, this effect must be 
inoperative in Y“°, owing to lack of some element which is present in the other Y 
chromosomes employed or perhaps to alteration of the pairing properties of Y“° as 
aresult of its ring shape. This ring chromosome differs from the other Y chromo- 
somes employed in that the bobbed locus is not present; all others have the normal 
allele of bobbed. The questions raised in this connection may be resolved partially 
through the use of a true l-armed Y™. 

The frequencies of detachments recovered using sc3+ Y as compared with se’!+ Y* 
differ by a factor of approximately 2. Multiplying the absolute detachment fre- 
quency by the frequency of X-4 linkage among the detachments shows, however, 
that there is not a reduction in X-4 exchanges with the doubling of X-Y exchanges; 
the value of 0.14 per cent is obtained for se’!- Y°, 0.12 per cent for sc- Y, and 0.10 
per cent for se-Y” for the dosage of 2,000 r. The lower detachment frequency 
found with use of a complete Y as compared to an acrocentric short arm of the Y 
must mean that the presence of a long arm of the Y lowers the frequency of ex- 
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changes between X and the short arm of the Y. This might likewise be interpreted 
as “protection,” with the possibility in this case that the ‘“‘protected” short arm 
might have a high frequency of exchange with the long arm. Investigations of ex- 
changes with free or with attached arms of the Y chromosome should be of value in 
answering some of the questions raised. 

Summary.—In attached-X females lacking a Y chromosome, most of the induced 
detachments involve exchanges with chromosome 4+. The presence of a Y chromo- 
some increases the frequency of induced detachment, 80 per cent or more of the 
cases resulting from exchanges between X and Y. With only Y® present, there is 
a much greater increase in detachment frequency. It is suggested that the long 
arm of the Y “protects” the short arm against induced exchanges with the X and 
possibly vice versa. With sc’- Y, se- Y“, or se’’+ Y*, there is a marked reduction in 
the frequency of exchanges between attached X’s and chromosome 4. This re- 
duction is not obtained with Y“°, possibly due to its shape or to its lack of the 
proximal part of Y*. 

* Work performed under Contract No. W-7405-Eng-26 for the Atomic Energy Commission. 
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GLOBAL VERSIONS OF THE CENTRAL LIMIT THEOREM* 
By RaLeH PALMER AGNEW 
MATHEMATICS DEPARTMENT, CORNELL UNIVERSITY 
Communicated by S. Bochner, July 16, 1954 


1. Jntroduction.—Our terminology and notation agree with those in the book! of 
Cramér, to which we refer for an exposition of the facts set forth in this introduc- 


tion. Let &, &, ... be independent random variables having the same d.f. (dis- 
tribution function) F(x). Thus, for each k = 1, 2,... 

Pr{&S2j = F(z), -2<2<-@, (1.1) 
where F(x) is a real monotone increasing function for which F(— ©) = 0 and 
F(+o) = 1. Let g(t), defined by 

lh) = _ @ a dF (zx), mm. rh <0), (1.2) 


denote the c.f. (characteristic function) of F(x). Wesuppose that 
So. x dF(r) = 0, So. x? dF(z) = 1, (1.3) 


so that F(x) has mean 0 and standard deviation 1. 
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ted The d.f. F,,(z) and the c.f. %,(¢) of the sum & + & + ...+ &, are such that 
rm on(x) = [¢(x)]”", and hence 
yee [e()]" = fr ~ e dF, (2). (1.4) 
> in 
The d.f. of the combination 
sin a Lee 
no- (1.5) 
~*” 
the 
p ig is then F,,(n/ x), and we denote this by F,,(x). Its c.f. is [g(n — 744) |", that is, 
m se “ mie S 
te [o(n-"1)]" = S72. e# dF, (2). (1.51) 
) in The hypotheses (1.3) imply that formulas (1.3) hold when F(x) is replaced by F,,(x). 
re- A special case of the central limit, theorem asserts that, for each individual x in the 
the interval —© <r< om, 
lim F(x) = (zx), (1.6) 
where (x) is the Gaussian d.f. defined by 
®(x) = (29) fe" du. (1.61) 
It is the purpose of this paper to investigate the constants (,'”’, C2'”’, . . . defined 
by 

C,? = f2. |F,(z) — (x)? dz. (1.7) 
If itis known that, for some p > 0, lim C,,‘”’ = 0, that is, 

lim f°. |F,(x) — (x)|%dx = 0, (1.8) 
the result is a global version of the central limit theorem which complements the 
local version (1.6) in which values of x are considered one at a time. After having 
proved that C,,°? —~ 0 by use of a different method, the author discovered that the 
global and local versions of the central limit theorem are much more closely related 
than he had suspected. As we shall show in section 3, it is possible to pass from 

one to the other by application of theorems on convergence of sequences of d.f. 
of We shall show that if p > !/2, then C,‘” + Oasn—> @. 
ma Among the questions which remain unanswered, one of the most interesting is the 
lis- question whether the convergence of C,,‘”’ to0 is monotone. It is easy to conjecture 
that the answer to the last question is affirmative when F(z) is a binomial or uniform 
1) distribution satisfying (1.3) and perhaps for every distribution satisfying (1.3). 
2. Existence of the Constants C,'” in (1.7).—We first suppose only that F(z) is a 
ind d.f. for which the first moment 
So. v dF(x) 2.1) 
2) a ean cree = oe 
exists. Existence of this integral implies existence (as finite numbers) of the in- 
tegrals in the left members and hence also the integrals in the right members of the 
3) equalit ies 
Se. (—2) dF(x) = f°... F(x) dz, So? x dF(x) = fo” [1 — F(x)] dx. (2.11) 
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This shows existence of the integral 
So. | F(x) — H(a)| dz, 2.2) 


where H(z) is the unitary function for which H(z) = 0 when x < 0 and H(z) = | 
when x 2 0. Since [F(x) — H(x)] is measurable and | F(x) — H(z)| < 1, it fol- 
lows that [F(a) — H(x)]isin class L, foreach p 2 1. If F(x) and G(x) are two df. 
having first moments and p 2 1, then both [F(x) — H(x)] and [G(x) — H(zx)] are 
in L,, and hence their difference [F(x) — G(x)] is in L,. The assumption that 
F(x) possesses a first moment implies that the functions /’,,(2) defined in section | 
also possess first moments. Since (x) possesses a first moment, we conclude 
that, if F(x) possesses a first moment, then the constants C,,‘” in (1.7) are all finite 
when p 2 1. 
We now use the second of hypotheses (1.3). If « < 0, then 


0< 2x°F(x) = f*. x2 dFly) < ft. ydFly) S$ fe. ydFly) = 1, (23) 
and hence 0 S F(x) S x~*. If x > 0, then 
0s 21 — F(x) = SP dF) SSL ydFiy) s Sf2. ydF(y) = 1, (24) 


and hence 1 — x~? S F(x) S 1. If G(x) is another d.f. having the properties in 
(1.3), then 0 S G(x) S x-* when x < 0, and 1 — x~-? S G(x) S 1 whenz > 0. 
Therefore, 







































|F(z) — G(a)| Sa, ot #0. (2.5) 


Let h(x) denote the positive function for which h(x) = 1 when | x| < land h(x) = 
a-2when |z| = 1. We then have 


| F(x) — G(x)| < h(z), —-o <r< om, (2.6) 


and hence, for each p > '/2, 
oe Te . 4p e 
\F(x) — G(a)|\? dx Ss (h(x) dx = cs 7 (2.7) 
—o —o 2p — 


Since conditions (1.3) are satisfied when F(z) is replaced by F(x) and by ®(x), we 
can replace F(x) by F,(x) and G(x) by (zx) in (2.6) and (2.7). Therefore, the 
constants C,,‘” in (1.7) are finite when p > !/». 

As we shall see, the main usefulness of formulas (2.6) and (2.7) lies in the fact that 
the dominating integrable function [h(x) ]? is the same for all pairs of d.f. F(a) and 
G(x) satisfying (1.3). 

3. The Condition lim C,'” = 0.—Let p > '/2. Because of facts set forth in the 
introduction, we can prove that lim C,,'”’ = 0 by proving the following more general 
theorem. 

THEoREM 3.1. Let p> '/s. If Go(x), Gi(x), Go(x), . . . ts a sequence of d.f. (dis- 
tribution functions) for which 


S-.. 2 dG,(z) = 0, S-...2? dG,(z) = 1, n=0,1,2,..., (8.11) 
and if, for each individual x in the interval —-” <x < @, 


lim G,(x) = Gy(x), (3.12) 


n—-> © 
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then 






























4 lim f°. | G(x) — G)(x)|? dx = 0. (3.13) 
é The hypotheses imply, as in section 2, that 

Lf. lGn(z) — Go(x)|? < [h(2) , (3.14) 
are 

sal where the function [A(x) ]’ is independent of n and is integrable over —~ <z2< ™, 
~ Hence (3.13) follows from (3.12) and the Lebesgue criterion of dominated con- 


uli vergence for taking Jimits under integral signs. This proves Theorem 3.1. 


‘ite The next theorem, a sort of converse of the previous one in which we do not need 
the hypotheses (1.3), shows that the local central limit theorem with conclusion 
(1.6) is a direct consequence of the fact that C,,°? > 0 as n > © and the elementary 
fact that &(2) is continuous and increasing over —© < 4 < o. 
3) THEOREM 3.2. Let G(x), Gi(x), Go(x), .. . be d.f. for which Go(x) is continuous and 
increasing over —-~ <x< om. If, for some p > 0, 
4) lim Sx |Gn(x) — Go(x)|? dx = 0, (3.21) 
sin then 
zi! ; 
lim G,(2) = Go(x) (3.22) 
9) uniformly over —-27 <4< om. 
pe Suppose the theorem is false. Then, for some ¢« > 0, there is an increasing se- 
quence of values of n such that for each n in this sequence there is a point 2, for 
which 
.6) | . 
|G,(t,) — Go(re)| 2 «. (3.3) 
Choose a constant b such that 
.7) 1—¢2S Q(z) S51, 228, (3.4) 
Pe and consider those values of n for which z, > b. Since 0 S G,(z,) S 1, (3.3) and 
‘he (3.4) imply that G,,(x,) — Go(z,) < —«, and hence 
G(x) S Giz.) —-e Sl —e, (3.41) 
hat - 
ae so that 
G, (2) Ss ] = ¢& b < x < Ln (3.42) 
the : 
ral From (3.4) and (3.42) we see that 
€ 
lis- (i(z) =~ G(x) >. 9? fee In 
Hence 
11) 
° Y Y Dp € f ‘ ‘ 
|Gp(a) — Go(x)|? dx 2 (x, — b) (5 (3.43) 
12) : : me ‘ ? 
If the sequence x, is not bounded above, we obtain a contradiction of (3.21). If 
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the sequence x, is not bounded below, a similar argument furnishes a contradiction 
of (3.21). 

It remains for us to consider the case in which the sequence x, is bounded, say 
—A <2, S A. We may replace the increasing sequence of values of n for which 
the points x, exist by a subsequence over which lim z, exists, say lim x, = a. The 
alternative case being analogous, we consider only the case in which 













Gr(Xn) ae Go(Xn) Pe € 


(3.5) 
for an increasing sequence of values of n. Choose @ such that 0 < @< «and Go(a) + 
6< 1. Then, for each sufficiently great n in the sequence, there is a unique point 
Yn for which y, > x, and 















Go(Yn) al Go(Xn) + 0. (3.51) 












When zx, S x S yy, We have 


G,(x) — Go(x) 2 Galan) — Go(x) 2 Golan) + 6 — Go(x) = Golyn) — Go(a) 2 0. 
3.52) 










Therefore, 


S20 |Gr(a) — Go(x)|? dx = Sf" |Go(yn) — Go(x)|? dz. (3.53) 





The fact that G(x) is continuous and increasing implies existence of b such that 
b> a, lim y, = b, and 


lim f2" |Go(yn) — Go(x)|? dx = f | Go(b) — Go(x)|? dx (3.54) 







as n becomes infinite over the sequence we are considering. The last integral in 
(3.54) is a positive number, say B, and it follows from (3.53) and (3.54) that 


Se |G(x) — Go(x)|? dx > B/2 (3.55) 



















for an infinite set of values of n. This contradicts (3.21) and completes the proof 
of Theorem 3.2. 
* This research was supported by the United States Air Force under Contract No. AF18(600)- 


685, monitored by the Office of Scientific Research. 
' Harald Cramér, Mathematical Methods of Statistics (Princeton, N.J., 1946). 


ON THE QUOTIENT OF AN ANALYTIC MANIFOLD BY A GROUP OF 
ANALYTIC HOMEOMORPHISMS* 
By W. L. Batiyt 
DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY 


Communicated by S. Bochner, June i4, 1954 


[ntroduction.—Let D be a complex analytic manifold, not generally compact. 
Let G be a properly discontinuous group of analytic self-homeomorphisms of 2. 
We assume first of all that the factor space D/G is compact, i.e., that there exists a com- 
pact subset © of D such that GE = D. Our second assumption requires some pre- 


liminary definitions. Let B be a fiber bundle over D, 7 the canonical projection of 
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B onto D, and G* a group of homeomorphic bundle maps of B onto itself. For 
q* « G* let 1g* be the induced map of the base space D onto itself. Suppose that we 
are given a homomorphism h of G into G* such that (1) rh(g) = g, g « G, and (2) 
oh h(g) (b) = b if gr(b) = (b),b¢€ B. Then we say that (B, h) is a G-bundle over D 


on 


‘he Now suppose that F is an analytic complex line bundle over D, and let G* be the 

group of analytic self-homeomorphisms of F which are bundle maps. Let h be a 
J homomorphism of @ into G* as above. Then (F, h) is called a G-complex line 
5) bundle. Let (F, h)* be the associated principal bundle. Note that if 2) « D, g «G, 
Rv and U is a sufficiently small neighborhood of z, then, with reference to a fixed pair 
int of co-ordinate neighborhoods containing U and gl, respectively, h(g): 2 —~'(U) > 


r'(gU) is given by h(g) (z, £) = (gz, f(z)°¢), where f is analytic and 0in U. We 
- speak loosely of f as the last co-ordinate of h(g). We obtain naturally from (F, h)* a 
1) principal G-positive-real-line bundle by mapping the nonzero complex numbers 

(* into the positive reals R+ by a: ¢— |¢|? and applying @ to the co-ordinate 

functions of F and to the last co-ordinates of the bundle maps h(g) in an obvious 
. way. Denote this principal G-bundle with fiber R+ by (| F|?, |h|*). It is easy to 
V2) show that (| F 2) always has a C® cross-section invariant under the group 
h|2(@). Suppose that (|F|2, |h|?) has a C® cross-section § invariant under 
|h| 2(@) such that the (1, 1) form 





9 














93) 
y = 100 log § 

rat 

is positive definite everywhere on D (note that y is well defined and invariant under 
; (7). ~~ we say that (F,, h) is positive and write (Ff, h) > 0. Our second assump- 
m4) tion on D and G is that there exists a positive G-complex line bundle over D. This is 
in all we assume. Our main result is ihat the factor space D/G is an algebraic variety 

imbedded in a projective space. More precisely, there exists a regular mapping ® of 
Pr D into a projective space P™ of appropriate high dimension M such that: 
59) 1. can be factored by way of D/G 
oof D/ 
)0)- cceeiitediamrmntge GU a: Fh 





where V is a homeomorphism. 

2. is locally biregular outside the set of fired points of nonidentity elements of G. 

3. &(D) zs a locally irreducible algebraic variety in P™. 
OF 4. The set of singular points of &(D) is of complex dimension < (n — 2), where D 
has complex dimension n. 

Notice explicitly that we do not assume G without elements having fixed points 
inD. Otherwise, our result would be contained in a result of Kodaira.! 

When D is a bounded domain in complex Euclidean space, the second assumption 
is always fulfilled. In fact, we may let F be the product bundle D X C, let 


ct. h(g)(z, ¢) = (gz, jo(z)~'+¢), where j,(z) is the determinant of the Jacobian of g at z, 
>. and let § be Bergmann’s kernel function. In this case, the mapping @ is achieved by 
m- taking as the homogeneous co-ordinates of ®(z); ho(z), . . . , har(z), where ho, . . . , har 
re are a linear basis for the automorphic forms of a certain suitable high weight ;? 


moreover, &(D) is precisely the variety given by the homogeneous algebraic equa- 
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tions among the homogeneous co-ordinates (ao, . . , Xs) satisfied by ho, . . , hag on D. 
This case is included in results obtained by H. Cartan,’ using an approach different 
from ours. A further case arises when D is a nonsingular algebraic variety and Ga 
finite group of analytic self-homeomorphisms of D. 

A Sketch of the Proof —The proof proceeds essentially along the lines of Kodaira’s 
proof of his result.!_ A difficulty lies in the fact that Kodaira’s proof depends on the 
theory of harmonic integrals for nonsingular, compact, complex analytic manifolds, 
whereas our quotient.space D/G is not, in general, nonsingular. Without attempt- 
ing to develop a theory of harmonic integrals on manifolds with singularities, we 
indicate how the principal features of the theory of harmonic integrals can be imi- 
tated in our situation. 

By a G-open set in D we mean an open set © in D such that gO = © for g « G. 
We have a metric on D invariant under G. Hence a small geodesic ball U’ about 
z e D is carried onto itself by the group G,, leaving z fixed, and, since G is properly 
discontinuous, this means that G, is finite. We may assume U so small that 
gU n U is empty unless g e G,. We call such a small geodesic ball G-small. It is 
clear that G-open sets of the form GU, where U is G-small, generate the ring of all 
G-open sets. 

Let S be a sheaf over D. Suppose that with each g ¢ G there is associated a com- 
patible endomorphism ¢g, of S, mapping, in particular, the stalk over gz isomorphi- 
cally onto that over z, such that 4, = ¢),¢,.. We wish by means of such endo- 
morphisms to give an identification of the stalks over points congruent under G. 
For this identification to be well defined, we must have ¢g, the identity on the stalks 
over points left fixed by g. But this may not always happen. Hence we define 
the subsheaf Sg of S by 


Sc = V {se8S,|¢,(s) = 8 for gz = 2}. 
zeD 


(Naturally, Sg also depends on the g,’s.) It is easily verified that Sg is a subsheaf 
of S. By a G-section of Sg we shall mean a section of Sg over a G-open set, invariant 
under the g,’s. Let M be the module of differential forms in a small neighborhood of 
zeD. LetgeG. Then with g there is associated the dual mapping g* which, in par- 
ticular, takes M onto the module of differential forms in a small neighborhood of g~!z. 
The mapping g* then gives an endomorphism of the sheaf S of germs of differential 
forms. We call Sg corresponding to these ¢g,’s the sheaf of germs of G-invariant 
forms, or simply G-forms. By aG-form we shall mean a G-section of Sg. Suppose, 
more generally, that (F, h) is a G-complex line bundle. Let My be the module of dif- 
ferential forms with coefficients in F in a small neighborhood V of z « D. Let 
xeMy. ThenyeFy @ M, where M is the module of differential forms over V and 
Fy the module of sections of F. Define (gp, »| Mr) = (hg-))|Fy) @ (g*| M). 
Then gy, , defines an endomorphism ¢, of the sheaf of germs of differential forms 
with coefficients in F. The resulting Sg is called the sheaf of germs of G-forms with 
coefficients in (F, h). By a G-form with coefficients in (/, h) we mean a G-section 
of this Sq. 

Let D be a divisor on D. Suppose that for each z « D there exists a neighbor- 
hood U of z such that the local equations of D can be chosen invariant under in 
@u. Then D is called a G-divisor. A G-divisor defines a G-complex line bundle in 
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a canonical way. Moreover, by the mth Jacobian bundle we shall mean the G- 
bundle (Ff, h), where F is —m X Jacobian bundle of D and the last co-ordinate of 
h(g) is j,(2)~™ with respect to a fixed pair of co-ordinate neighborhoods. It will be 
seen later that the mth Jacobian bundle is well defined for m divisible by an integer 
mo Which depends only on D and G. 

The following are easily proved facts at our disposal: 

1. The operators d and 6 carry G-forms into G-forms. If (F, h) is a G-complex 
line bundle, 0 is already well defined and d can be defined® with respect to a C” 
cross-section of (|| ?, |h|?) invariant under |h|2(@); then these so-defined opera- 
tors carry G-forms with coefficients in (F, h) into G-forms with coefficients in (F, h). 

2. If Wis the positive definite form of type (1, 1) of the Introduction, dy = 0, so 
our structure is Kahler. 

3. An integral fo = Joya ¢ of a G-form ¢ of rank 2n over D/G may be de- 
fined in an obvious way. Let 7 be a 2n — 1 G-form of class C”. Then we have 
JSo/a dn = 0. The inner product (¢, 7) of G-forms ¢, 7 of same degree is defined 
by (v, ) = Soya ¢a*h. We then have the adjointness formula for G-forms, 





(dy, n) = (¢, 5n), 


and the adjointness formula for G-forms with coefficients in a G-complex line 
bundle,® 
(Og, n) = (¢, dn). 


4. The following local properties hold: 

(i) The Poincaré lemma with respect to G-forms. 

(ii) Weyl’s lemma.’? Here we may take as Laplacian A = 6d + dé for ordinary 
G-forms and O = od + do for G-forms with coefficients in a G-complex line bundle. 

5. The sheaf of germs of G-forms ¢ C” (over G-open sets) is fine. We simply 
show the existence of a G-invariant partition of unity subordinate to any finite 
covering of G-open sets. 

6. The following orthogonal decomposition formula holds: 


Ck UF, h) = 0C#*~(F,h) @0C#* + (F, h) © HE “(F, h), 


where CZ; ‘(F, h) is the module of G-forms of type (p, q) with coefficients in (F, h) and 
H@"(F, h) is the module of G-forms C”, g, of type (p, q) with coefficients in (F, h) 
such that (0d + d0)¢ = 0. 

7. From 3 and 4(i) we get 


H4(D, 02 (F, h)) & ZB *(F, h)/oCz 4 ~ '(F, h), 


where HZ(D, 22,(F, h)) is the direct limit of the cohomology modules of the nerves 
of the coverings of D by G-open sets whose coefficients are G-sections of O2(F, h), 
the G-sheaf of germs of holomorphic G-forms of degree p with coefficients in (F, h), 
and where 

Ze 7 (F,h) = {ee CRF, h)| de = Of. 


8. H2(D, Q2(F, h)) S HR (FP, h). 
9. Let m; = order of G,. Let mo = ym (m;). mo is a finite integer because of 
the compactness of D/G. Then the Nmoth Jacobian bundle is well defined (N is 
an arbitrary integer), since (j,(¢))”" = 1 if gg = ¢. 
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10. Suppose (F, h) is a G-complex line bundle, and let mK be the moth Jacobian. 
bundle. If m(F, h) — mK > 05 


He(D, NG(F,h)) =0, q2>1. 










Using 1-10, we can imitate Kodaira’s proof! and obtain the mapping ® satisfying 
1, 2, and 4 of our main result. There is a theorem of Osgood* * which implies that 
a mapping of a neighborhood V of 0 in C” into C” given by n analytic functions in 
V such that the inverse image of a point is discrete is an open mapping. Using 
this, we show that (D) is a locally irreducible analytic manifold. Using Chow’s 
theorem," '! we obtain part 3 of our main result. 













* Refer to the literature indicated for the notation used. 
+ The author was a National Science Foundation predoctoral fellow during the preparation of 
this paper. He wishes to express his warmest gratitude to Professor K. Kodaira for many helpful 
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ON EXTREMAL QUASI-CONFORMAL MAPPINGS. 1 
By Murray GERSTENHABER* AND H. E. Raucu 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF PENNSYLVANIA 


Communicated by E. Hille, June 26, 1954 


1. Introduction: Quadratic Differentials—Riemann, in his fundamental me- 
moir, “Theorie der Abel’schen Funktionen,” observed that two algebraic, i.e., 
closed, Riemann surfaces of the same genus p (= !/2 connectivity) could not neces- 
sarily be mapped on one another conformally but that, remarkably enough, the 
totality of “conformal classes” of these surfaces depended in an imprecise way on 
3p — 3 continuous parameters, which he called moduli. Among the various at- 
tempts to make the dependence more precise, current attention centers on an ab- 
stract and speculative attempt by Teichmueller! based on the remark—a simple 
consequence of the Riemann-Roch theorem—that the number of complex-linearly 
independent analytic quadratic differentials is the mysterious number 3p — 3. 
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ian. A quadratic differential on a Riemann surface S is a covariant tensor a(z) such 
that a(z) dz? is formally invariant under change of local uniformizing parameter. 
It is analytic if the a(z) are so. 
Teichmueller claimed that, even if two surfaces, S; and So, of the same genus are 
- not conformal, one can “get from S; to S,” via an analytic quadratic differential and 
a a special quasi-conformal mapping. To be more precise one needs the concept of 
es dilatation of a homeomorph ism f? Iff takes the form w = w(z) in local parameters 
‘ zabout p e S; and w about f(p) ¢ So, then dw/dz at p depends on the approach to p 
= unless f is conformal at p. For fixed p the ratio of the supremum of | dw/dz| over 
Me all directions to the corresponding infimum is D,/, the dilation at p. If supy.s, Dp f 
= Df < ~, then f is called quasi-conformal and Df is its dilation. Clearly Df > 1 
and Df = 1 if and only if f is conformal. Now Teichmueller claimed: 
mn. of A. In every homotopy class of homeomorphisms of S, on S2 there exists an essen- 
ipful tially unique quasi-conformal f, called extremal, for which D, f, where defined, is a con- 
stant D = Df. 
316, 


At its exceptional points this f is conformal; to any other p e S; corresponds the 
re unique line element +dz, for which | dw /dz\ assumes its supremum. Clearly, since 
every infinitesimal circle about p goes over into an ellipse about f(p), one goes from 
z to w by stretching z along every parallel to +dz, by D, the ratio of major to mor 
axis. Teichmueller’s next assertion is: 

B. The line elements +dz, defined by an extremal quasi-conformal f satisfy » = 





ee a(z) dz,? > 0, where n is an analytic quadratic differential on S,.* 
Finally Teichmueller claimed, and Ahlfors‘ affirms, that A and B can be deduced 
rt 1, by solving the obvious minimum problem Df = min among all f in a given homotopy 
E class. But experience shows that derivatives are harder to minimize than integrals 
er of derivatives, and so the aim of this note and the following is to formulate a problem 
in the calculus of variations whose solutions satisfy A and B. Questions of existence 
are then reduced to the study of double integral problems as treated by Morrey. 
Schiffman,® et al., but these are not broached here. 
2. Metric Considerations and the Dirichlet Integral—Now one may always intro- 
duce a conformal metric 7 = \ dwd#on S». Referring it to the co-ordinates of S; 
by means of f, one obtains 
f*n = \(w(z)) dw(z) d®(z) = a(z) dz? + a(z) dz* + b(z) dz dz = 
E dx? + 2F dx dy + Gdy?, (I) 
where a = [(E — G)/2]-— iF and b = (E + G)/2, as is easily verified (w(z) is not 
analytic, in general). For future convenience the terms in the middle expression are 
mul called of type (2, 0), (1, 1), and (0, 2) respectively, and it should be noted that the 
_ decomposition is independent of the local parameter z on S,; in particular, a(z) dz? 
— is a quadratic differential; and f is conformal if and only if a(z) is identically zero. 
the But even if f is not conformal but merely quasi-conformal—extremal or not—the 
i line-element argument of section 1 is still valid; and, on rewriting relation (1) as 
ab- f* = b dz dz + 2Re {a dz*t, (2) 
rple 
arly one sees that B with the word “analytic” deleted is true for an arbitrary quasi-con- 
formal f. 
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Furthermore, it is well known in the theory of minimal surfaces and particularly 
Douglas’ justly acclaimed solution of the Plateau problem® that the condition that 
a(z) be analytic is a necessary condition for extremizing the Douglas-Dirichlet functional 


l 
I( f*n) = ff (EB + G) dx dy = = ff bdz A dz, (3) 
“ S, = 8, 


where A denotes exterior product. To formalize the argument one needs the fol- 
lowing: 

Definition: A homeomorphism f of S; on S for which, given n on So, the term a(z) 
dz? in the decomposition of {*n is an analytic quadratic differential is called harmonic 
with respect ton. Then one has: 

THEOREM 1. (Given n on So, let f be a homeomorphism of S, on Se such that for 
any homeomorphism g homotopic to f, I(f*n) < I(g*n). Then f is harmonic relative 
to n.? 

A new proof of this known theorem is given in section 3, the proof consisting 
essentially of a calculation of the first variation of J(/*) in complex formalism (z 
and 2). This is needed in the remainder of the investigation and is otherwise very 
interesting per se.°® 

Finally, the property of extremal quasi-conformality of a homeomorphism f is 
independent of any metric 7. In the next note, therefore, the following extremal 
problem will be proposed: For any 7 on S, find a harmonic homeomorphism f 
minimizing I(f*n); then maximize I(f*n) over all 7. It is asserted, subject to the 
verification of certain continuity theorems, that a solution f of this problem is the 
extremal quasi-conformal mapping of A which, thanks to Theorem 1, will satisfy B. 

3. Proof of Theorem 1.—Any variation of f may be obtained by preceding it by 
a homeomorphism of S; onto itself. Let S; be covered by a finite number of com- 
plex co-ordinate neighborhoods N; with corresponding local complex parameters 
z; A homeomorphism H of S, onto itself will be called “small” if, given any point 
p of S;, there exists a co-ordinate neighborhood N among the N; such that both 
p and H(p) are in N. If z is the local complex parameter for NV, then H may be 
considered locally as a function of z and is representable by z ~ H(z). It is con- 
venient to set H(z) = z + A(z), and A(z) is then the “variation” of z. The homeo- 
morphism obtainea by following H by f will be denoted by Hf. It produces an- 
other form (Hf)*n on S;. In the neighborhood N, f*n may be considered as a func- 
tion /*n(z), and, if it is so denoted, then (Hf)*n (z) is f*n(z + h(z)). The theorem 
is proved by comparing the Dirichlet integrals of f*y(z) and f*n(z + h(z)). 

From relations (2) and (3) one sees that the form f*n(z + h), is 


2Re(a(z + h) d(z + h)*) + b(z + h) d(z + h) d(z + h), 


and to find its Dirichlet integral it is necessary to compute the coefficient of dz dz in 
this. If A, like all homeomorphisms, is assumed to be continuously differentiable 
at all except a finite number of points, with those exceptions, d(z + h) = dz + 
h,dz + h,dz = (1+ h,) dz + h;dz.°_ Therefore, 


f*n(z + h) = 2Rej\a(z + h)[(1 + h,)? dz? + (hz)? dz?) + 
Qa(z + h) (1 + h,) h; dz dz} + b(z + h) d(z + h) d(z + h). 
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The last term is purposely left unexpanded, for, if f is a differentiable function of 
z, the Dirichlet differential of df df is just '/.df A df. The Dirichlet differential of 
f*n(z + h) is therefore 


tRe[a(z + h) (1 + hz)hg] dz A dz + (i/2)b(z + h) d(z +h) A d(z + h). 


The hypothesis on the homeomorphism /f implies that if H is any small deformation 
of S, onto itself, then /((Hf)*n(z)) — I(f*n(z)) = I(f*n(z + h)) — I(f*n(z)) = 0. 
As (7/2) SSs, b(z + h) diz + h) A d(z +h) is identically equal to (7/2) SSs, 
b(z) dz A dz = I(f*n(z)), it follows that 


I((Hf)*n(z)) — I(f*n(z)) = 27 ffs, Rela(z + h) (1 + h,)hg] dz A dz. = (4) 


This is the exact variation of the Dirichlet integral J and will be denoted by 4/. 

It will now be shown that from 6/ > 0 for all H, it follows that a is analytic in a 
neighborhood JN, i.e., that the component of f*n of type (2, 0) is an analytic quad- 
ratic differential. This will prove that f is harmonic relative to 7. 

Let p be a point of N, where f*n is continuously differentiable. By considering a 
possibly smaller neighborhood of p, it may be assumed that f*y is in fact continu- 
ously differentiable throughout the neighborhood. This smaller neighborhood 
will henceforth be denoted by NV. If f*n is expressed here in the form 2Re[a(z) dz?] 
+ b(z) dz dz, then a(z) has continuous first partial derivatives with respect to both 
zand 2 throughout N. 

Let U be a smaller neighborhood of p topologically a disk and with analytic 
boundary, such that N > U. It is then possible to define in N a function h(z) 
having the properties that (i) A(z) is continuously differentiable throughout N; (ii) 
for any complex number a with | a| < 1 the transformation z ~ z + ah(z) defines 
a homeomorphism of NV onto itself which maps U into U and is the identity outside 
U. (In particular, h must vanish outside l’.) The homeomorphism of N onto it- 
self defined by z > z + a h(z) may be extended to all of S; by letting it be the iden- 
tity outside U and will be denoted by H“(z). 

For each H the integral expressing 6/ need be extended only over U’, and, by the 
hypothesis on f, one has that 


2i [fu Rela(z + ah) (1 + ah,)ah;] dz Adz > 0. 


As a is continuously differentiable in U, the expression on the left is 27 
SSu Refa(z)ah;] dz A dz plus terms of lower order with respect to a as a — 0. 
The dominant term must be non-negative. As a may have any argument, it follows 
that {ft a(z)h; dz A dz = 0, and this holds for all functions h satisfying condi- 
tions i and ii. By Stokes’s theorem, Jfp ah;dz A dz = —foy ah dz — Sfu ha; 
dz \ dz. But h vanishes on 0U; therefore, [fy ha; dz A dz = 0. The condi- 
tions on h are not stringent, and it follows that a; = 0 in U, that is, a is analytic in 
U. Hence a(z) dz? is an analytic quadratic differential in a neighborhood of any 
point where f* is continuously differentiable. At the finite number of exceptional 
points a is, by the assumptions on f and n, bounded. One may then extend a(z) 
dz’ to be an analytic quadratic differential on all of S;. . The component of f* of 
type (2, 0) is therefore an analytic quadratic differential on S,, whence f is by defini- 
tion harmonic with respect to n, and the theorem is proved. 
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* A Frank B. Jewett postdoctoral fellow at the time this work was begun. 

‘ Abhandl. preuss. Akad. Wiss., Jahrgang 1939, No. 22, 1940. 

2 Henceforth all homeomorphisms, when expressed in terms of local parameters, will be assumed 
at least C!, with the possible exception of a finite number of points. 

3 The exceptional points are thus the 2p — 2 zeros of n, where f is conformal. 

* J. d’analyse math., 3, Part 1, 1, 1953. 

5 Cf. Ann. Math., 48, 274, 1947. 

6 Cf. Am. J. Math., 61, 545, 1939. 

7 For the purposes of proof » will be assumed at least C! and f, C*, except at the exceptional 
points where the first partials will be assumed bounded. 

8 The variational condition in the Plateau problem can be remarkably elucidated with these 
ideas. 

® hz and hz denote the formal derivatives 0h/dz and Oh/03. 


REGRESSION ANALYSIS OF TIME SERIES WITH STATIONARY 
RESIDUALS* 
By UtF GRENANDER AND Murray ROSENBLATT 
UNIVERSITY OF STOCKHOLM AND UNIVERSITY OF CHICAGO 
Communicated by Marshall H. Stone, July 15, 1954 
Let y, = x, + m, be a complex-valued stochastic process with mean value Hy, = 
m, and residual x, weakly stationary, that is, the covariances Ex,x, = r; — , exist and 
depend only on the time difference 1 — rv. The parameter set 7, ¢t ¢ 7’, is the set 
of integers. The regression m, is assumed to be of the form 


p 
m= 2) ag, (1) 
k=1 
" 
¢1' ) 
where the regression vectors 9 = | : k = 1,..., pare assumed to be known. 
yg ) 
(k) 
YN 


The spectral distribution function! of the residual x, is assumed to be absolutely 
continuous, and the spectral density f(\) positive and continuous. Assume that a 
time series (a finite realization of the process) y;, ..., yy is observed and that one 
wishes to estimate the regression coefficients from these observations. 

We shall consider linear unbiased estimates of the regression coefficients. Let 


mM, Cy 
= , c= 
my Ry 
Relation (1) can then be written 
m= 4, 
where @ = (9, ..., g®). We could estimate ¢ by the vector c,* that mini- 


mizes the quadratic form? 


(y — &c)'(y — ®c). 














n. 
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We call the minimizing vector c,* = (&’6)~'’y the least-squares estimate. This 
vector is unbiased, that is, 


Ec,* = c, 
and the covariance matrix 
E(e,* — c) (e,* — c)' = (&’6)—'6’RO(6'6)—!, (2) 
where R is the covariance matrix R = {r;-4;j,k = 1,... ,N}. The covariance 


matrix 2 is nonsingular, since f(A) is positive. However, the best linear unbiased 
estimate of cis the Markov estimate, 


Co* = (®’R-'46)—'6’R-ly, 


in the sense that the covariance matrix of any other linear unbiased estimate of c 
is greater than or equal to that of c)* 


E(co* — c) (co* — c)’ = (#’R-'@)—". (3) 


The aim of this paper is to examine the behavior of these two estimates as N > 
and, in particular, to find out when the least-squares estimate is as good as the 
Markov estimate asymptotically. This is of especial interest in two contexts. 
If we do not know the covariance structure of x,, we are led to the least-squares esti- 
mate, since the Markov estimate requires knowledge of the spectrum of z,. On 
the other hand, we may know the covariance matrix R, but it may be too costly 
or tedious to carry out the inversion and compute R~!, in which case we are again 
led to the least-squares estimate. 

The results are obtained by what we might call a generalized harmonic analysis 
of the regression vectors. In order to carry this out, we have to assume that the 
regression vectors satisfy certain conditions. Let 


N 
= F le? k= 1,...,p 
j=1 

We assume that dy“ — © as N > o (if x, is normal, such a condition is required 
if we are to be able to estimate c consistently) and that 

Dyin 

: 
by 


as N —~ o for any fixed h. We also assume that the limits 


N 





(9) k 
itn Gt” 
lim “ = M,%» 
N—> by Gy 
exist (let ¢,” = 0 if <0). Consider the sequence of matrices 


Mya {Mer Reh ph ee... =H. 


rhis is a positive definite sequence of matrices in the sense that for any p vector z 
the sequence 


2’My,h=...,-1,0)1,..., 
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is a positive definite sequence. The sequence of matrices can then be seen to be 
the Fourier-Stieltjes coefficients 
h 7 
M, = fz, ¢" dM(a) 

of a nondecreasing matrix-valued function M(A), MQ) > M(d2), if Ar > Av. It is 
assumed that My) = M = M(mr) — M(—7) is nonsingular. The conditions im- 
posed on the regression vectors are easily seen to be satisfied by the usual sort of 
regression, such as a polynomial or trigonometric regression. However, one should 
note that they are not satisfied by an exponential regression. 

Let Dy be the diagonal matrix 


We consider the first-order terms of the covariance matrices (2) and (3) as VN > o, 
that is, 


DyE(e,* — c) (e,* — ¢)'Dy = Dy('®) | DyDy —'®’R€Dy~'Dy(®'8) Dy 
and , 
DyE(eo* — c) (eo* — c)'Dy = Dy('R-'@)—'Dy. 
THEOREM 1. We have 


lim DyK(e,* — c) (ep* — c)'Dy = Mf, f(—d) dM(\) M-! 29 (4) 


No 
and 
: ’ - ™ 1 ~1 
a. Dy E(eo* — c) (co* — c)’Dy = ( fen am) 27. (5) 
Clearly, 


M-J*_ f(—r) dM(a) M-' > (J po aM) 


Both limiting matrices are nonsingular. This theorem is obtained by approximating 
f(A) above and below by finite positive trigonometric polynomials or inverses of such 
trigonometric polynomials. 

It is of particular interest to find out when the two limiting matrices (4) and (5) 
are equal, that is, when the least-squares estimate is asymptotically as good as the 
Markov estimate. 

We call the set of points \ on which the matrix-valued function (A) increases 
the spectrum S of the regression. 

THEOREM 2. There is a unique minimal decomposition of the spectrum S into q 
nonoverlapping sets Ei, ..., 2a,q < p, such that the increments 


M; = M(E,) =Sz, €MQ) 


satisfy 


M iM —1M, —= 5;;M,. 
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This decomposition is minimal in the sense that, given any set E;' of a decomposition 
N,',..., Eq’ satisfying these conditions, there is a set BR, ¢ E;’. 
THEOREM 3. The equality 


a a? 
M~' f=, f(—d) dMQ)M~ = (J amy) (6) 
=? f( — ) 
holds if and only if f(A) = ¢: > 0 is constant on each of the sets Ei,i1 = 1,..., q. 
TuroreM 4. Equality (6) holds for all positive continuous f(r) if and only if the 
elements E;,i = 1,..., q of the regression spectrum are all points. 


It is of especial interest to note that equality (6) holds for all positive continuous 


f(d) in the case of a polynomial regression or a trigonometric regression or a more 


general mixed regression where the regression vectors are of the form 


( ith 
eo” =e, y= O;1,...,0 ] 
§ + 1+ the 
ar ae t’e' ; 


y=0,1,...,p (7) 


(m+... Im -1 t+ m—1 + vr) itd \ 
g. T Pm—1 Lyk) Sage t’e m | yp —{) hi h ory) 


Thus, if the regression vectors are of the form (7), the least-squares estimate is 
asymptotically as good as the Markov estimate. 

Perhaps the case of greatest interest is that in which the process x, and the re- 
gression vectors ¢ are real. Then f(A) = f(—A), and M(A) = M(—d). It is 
convenient to introduce N(A) = M(A+) — M(—A—-),0 <A <27,N = N(m) = 
M. The limiting matrices (4) and (5) can then be written as 


2aN-' fo" f(A) dN(A)N- 


r 1 ; : —1 
2 (J. f) aN 0») ‘ 


respectively. The analogues of Theorems 2, 3, and 4 can be seen to hold with 
N(A) in place of M(A). It can be seen that the least-squares estimate is asymptoti- 
cally as good as the Markov estimate if the regression vectors are of the form 


and 


er = 008 tA, v= 0, be oo ey 1p 
| a eo ae = — j , 
1p. Pn—i te TE oe 2 oom Mn, yv=0,1,. 1Dm (8) 
gr” =U sinih, v =0,1,..., 2m 
- be " = rage ae 
Ps ee 2Pm—1 m 1 + ») = ?’ sin tdms y=QO0 L nei 


Obviously the sine terms do not appear if \; = 0. 

It is of considerable interest to find out whether these asymptotic results are 
effectively valid for small or moderate sample size. A few numerical computations 
that have been carried out seem to indicate that they may be. It would be of espe- 
cial interest to carry such computations further. 

The results cited in this paper are discussed and proved in some detail in a forth- 
coming book by the authors on statistical analysis of stationary stochastic proc- 
esses. 

* Based in part on research supported by the Office of Naval Research at the Statistical Re- 
search Center, University of Chicago. 
1S$ee J. L. Doob, Stochastic Processes (New York: John Wiley & Sons, 1953). 
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2 A’ denotes the conjugated transpose of the matrix A. 

® Discussion of related problems is to be found in the two following papers: U. Grenander, “On 
the Estimation of Regression Coefficients in the Case of an Autocorrelated Residual,’ Ann. Math, 
Stat., 25, 252-272, 1954; U. Grenander and M. Rosenblatt, ‘‘An Extension of a Theorem of G, 
Szego and Its Application to the Study of Stochastic Processes,’’ Trans. Am. Math. Soc., 76, 112- 
116, 1954. 


ON POSITIVE GREEN’S FUNCTIONS 
By G. A. Hunt 
CORNELL UNIVERSITY, ITHACA, NEW YORK 
Communicated by S. Bochner, July 23, 1954 

Consider on Euclidean space E* an elliptic differential operator Au = >0(0 da;): 
(a;,(Ou/Ox;)), with the a, differentiable and symmetric in the indices. It is 
desirable to have a criterion in terms of the a;; for deciding when A has a positive 
Green’s function vanishing at infinity—that is, a function K(a, y), tending for 
each fixed y to zero as x tends to infinity, such that 

A Si K(x, yfy) dy = —f(x) 
for every smooth function f which vanishes outside a compact set. 

I give here a partial solution. The method does not yield a great deal, but it 
has the good feature that it applies to similar problems on certain homogeneous 
Riemannian spaces. 

Let |a| = (a? +... + 2,2)’, and for a function f on E* define 

r) = fof (ex) do, a a 
where OQ is the rotation group on E*, do is the Haar measure on O attributing total 


measure 1 to O, and z is any point reas \z| =r. Itis clear that 
far(|a| asf, JOf(x) " 
d 
pe \ — u \ oa; f as 


whenever the positive measure yw has rotational symmetry. 

We suppose given a positive function A(r) such that A (|x| )doa,(a)rAaAy 2 DM 
for xe E* and all \,, and writeI = f° r'~*A(r) dr. 

THeorEM 1. If I < ©, there is a solution ¢ of the problem 

Au(z) = Ofor|x| >1, u(x) >1as|a| \ 1, u(r) > Oasr—> o. 

For each p > 1 ate is certainly a solution ¢g, of the problem: Au = 0 for 1< 
| x| < p, u(x) > las || \ i u(x) > 0 as 2 7 p. Let us continue ¢, to all of 


E* by making it zero for | a| p and one for | |x| <1. Ifl<p<ao,wehave0 $ 
%, = v- & 1 by the sas of the maximum. Consequently, g(x) = lim ¢,(2) 
ro @ 


exists for all x, lies between zero and one, and satisfies Ag = 0 on the region |x| > | 
(Harnack’s theorem). Moreover, g(x) > | as | 2| — 1, because g(x) S g(x) <1. 
It is conceivable, however, that yg may be identically one; this, essentially, is what 
we must disprove. 

Lemma. Jf B(s) > 0 and b = f,* B(s) ds < &, the problem 
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f Ba C3) (=) ds = minimum, Y(r) = a, v(s) > 0ass—> ~ 


has a solution, and the minimum is a/b. 
Making the change of variables s’ = f° B(u) du, ¥’(s’) = (s), the problem be- 


comes 

b ‘\ 2 

dy’\? 

| (+) ds’ = minimum, v'(0) = a, y'(b) = 0, 
0 as 


whose solution is ~’(s’) = a — as’/b. 
Continuing the proof of the theorem, I note that ¢, also solves the problem 


Ou Ou ee ju(z) > las |z| \ 1, 
2a ay 5 as = minimum, 


u(x) —> Oas |x| 7 », 


and that the minima of the various problems for p 2 2 are bounded above by the 


minimum VV for p = 2. Let r> 2 be fixed, and suppose p> 2. Then 


M2 > au oer dx =f i i "> Po dx = 
“ v; 
el >? 


rh S>d 
i d¢,\" 4. fa (03,(|2|)\? 
A-1(| >| » ?) dx = |) 4 p . as 
by A ‘(| 2}) i) are our, ( da 


r ae M4 Pe 
: »—1 §4@,(8))? 
) A-1(s)s* 1 OPAL ls, 
“4 (8) ‘a f as 


where w is the (k — 1)-volume of the sphere | 2-| = 1. The last integral is not less 
than the minimum in the problem 


| ; 1 (ot (2) ds = minimum V(r) = ¢, (r) Y(s) > 0ass—> @ 
r ; ocd ds ‘ Pr | - , _ ) 


So, taking the lemma into account and rearranging, 


M(°? 3 
?,(7) <3 | A(s)s'~* ds’ 
210 J f 


Letting p— ~, one obtains the same inequality for g(r). This proves the theorem 
and even provides a majorization of 2(r). 

CoroLtary. There rs 8 < 1 such that g(x) < 8 for |x| = 2. 

So far the argument is quite rigorous. I only sketch the proof for the Green’s 
function, however, because I require the following statement, whose proof is not 
published. 

If the a;; and 0a;;/0x,; satisfy the preceding hypotheses and also global Lipschitz 
conditions of order one, there is a stationary Markoff process X(¢, w) (f = 0, ina 
probability space Q) such that (1) for each w the function X(t, w) of ¢ is a con- 
tinuous curve in E*; (2) the transition probabilities have the form 


Pr | X(t, w) € F| X(s,w) = a} = fp p(i— ss, 2, y) dy, s<l; 
(3) p(t, x, y) is continuous in ¢, x, y for d > 0, and Op/dt = A,p. 
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the equation 0u/dt = Au and that the Green’s function for A is 


(x,y) = fr” pit, 2 


provided that the integral e de Sv dy f° ptt, x, y) dt is finite for every bounded 
set F. Also, the g(x) constructed above is the (conditional) probability that a path 
starting from x will sometime hit the sphere | y| = 1, and the expectation of the 
total time 7'(x, w) a moving point X(t, w) starting from x will spend in a set F is 


&| T(x, w)| = Sr dy f° pit, x, y) dt. 


I shall now majorize & | T(x, w)} when F is a bounded set—say, the ball | y| < 1 
First, note that the expectation of the time it takes a moving point, starting from 
z on the sphere |z| = 1, to reach the sphere | w| = 2, is bounded by a constant y 
independent of z. The argument now proceeds in a way which can be made rigor- 
ous. Consider a path X(t, w) starting from 2; if it hits F, which it does with prob- 
ability g(x), let T:(x, w) be the first ¢ for which |X (é, w) | < 1 and 7,(2, w) the first 
time ¢ (which is finite with probability one) after 7,(x, w) for which | X(é, w)| = 2. 
If X(t, w) is in F for some t > T(x, w) (an event of conditional probability less then 
B), let T(x, w) be the first ¢ > T2(x, w) for whic )| = 1 and 7;,(2, w) the first 
t > T;3(x, w) for which | X(¢, w)| = 2. Taking the corollary into account, one sees 
that with probability one this process ends in a finite number of steps NV (a, w). It 








is clear that 


N(2,w) 


T(x,w) S YS (Ty(x w) — Tx~1(2, w)): 


1 


so, making use of the Markoff property, 


Sr dy So pit, x, y) dt = &{ T(x, w)} S (a) > B" y = d¢(z), 
n=0 
with 6 independent of x. 
Consequently, 


Sv dy Jo do fo” pit, ox, y) dt S 52(| x!) 


and from this one can deduce that fo K(x, y) do > 0as | a| +o. Thus we have 
sketched a proof of the following theorem 

THEOREM 2. If the aj; satisfy the sail hypotheses, then A has a positive Green’s 
function K(x, y) = Jo” pit, x, y) dt and fo K (ox, y) do > 0as | | > oo, 


| 
x) > -, 











It is not hard to show, granted this much, that p is the fundamental solution of 























A SIMPLE DEFINITION OF ANALYTIC FUNCTIONS 
AND GENERAL MULTIFUNCTIONS 


By Kart MENGER 


ILLINOIS INSTITUTE OF TECHNOLOGY 


Communicated by Marston Morse, March 17, 1954 


Weierstrass defined a general (multivalued) analytic function as a class of prop- 
erly related power series; Riemann introduced it as a single-valued function the 
domain of which is what is now known as a “Riemann surface.’’ In Weierstrass’ 
theory and its modern variants, single-valued functions, which are classes (of a 
certain kind) of pairs of numbers, are not special cases of multivalued functions. 
The latter are classes the elements of which are functions and not pairs of numbers. 
Riemann’s approach presupposes the definition of a Riemann surface. For a given 
function the surface is not, in general, uniquely determined. In the present paper, 
Weierstrass’ definition will be modified by a procedure which lends itself to sig- 
nificant generalizations. ! 

Analytic Functions and Multifunctions——An analytic function is a class f of 
ordered pairs of complex numbers such that (1) no two pairs belonging to f contain 
equal first, and unequal second, members; (2) the domain of f—briefly, ‘“dom f’’— 
is a connected open set in the Argand plane; (3) / is differentiable at each point of 
dom f. If the pair (z, w) belongs to f, then w is denoted by f(z). If F is any class 
of ordered pairs of numbers, then f & F will symbolize the fact that F includes 
all pairs belonging to f. In this case, F is called an extension of f, and f a restriction 
of F. 

One simple concept, that of a motion, is all that is needed for the introduction 
of richer classes of pairs of numbers under the name of ‘analytic multifunctions’ — 
I avoid the term ‘multivalued function,’’ since “function,’’ as defined above, has 
the connotation of single-valuedness. A motion in F is an ordered pair (¢, 7) 
[t, t:] of complex-valued continuous functions, defined for every ¢ of the real interval 
ty <4 < t, and such that every pair (¢(¢), n(¢)) belongs to F. 

I now define an analytic multifunction as a class F of ordered pairs of complex 
numbers such that 

I. Each pair in F belongs to an analytic function & F; 

II. For any two analytic functions & F, fo, and f,, there is a motion (¢, 7) [0, 1] 
such that 

a) For some e > 0, the segmental motions (¢, 7) [0, «] and [1 — e, 1] belong to 


foand f,, respectively. 


b) If 0 < t < 1, then, for some « > 0, the motion (f, 7) [t — «, t + e] belongs to 
some analytic function & F. 

In I, it may happen that ¢(0) = ¢(1) and n(0) ¥ n(1), or even that ¢(0) = ¢(1) 
and 7(0) = n(1), without fo and f; being identical near z. = ¢(0). In this case, due 
to a well-known identity theorem concerning analytic functions, fo(z) # fi(z) for 
every z in a punctured circle about 20. By the method of Poincaré-Volterra, it can 
be shown that for every 2 the set of all pairs in F having the first member 2 is de- 
numerable. Moreover, one can define when two motions are homolopic in F and 
can introduce the homotopy group of a multifunction. If f is an analytic function 
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¢ F, there exists at least one maximal extension of f in F, that is, an f,; which is 
saturated with regard to the property of being an analytic function 2 f and & F. 
A multifunction F is a function if and only if every function & F has exactly one 
maximal extension in F, namely, F itself. 

Let f be a restriction of F and z be a point in the boundary of dom f. The point 
2) belongs to what may be called the F-boundary of dom f, if there are no two ana- 
lytic functions fo and g such that (1) fo Sfand fo Sg SF; (2) 2 belongs to dom 
g and to the boundary of dom fo. This reduction of the boundary to the F’-boundary 
reflects the classical fact that the only significant points on a branch cut are the 
branch points. The point 2 is an F-branch point of f if F contains three analytic 
restrictions, 91, g2 Sf, and g3, such that (1) z belongs to the boundaries of dom 4; 
(¢ = 1, 2,3); (2) the sum of the three domains contains a punctured circle about 
z; (3) g: U ge and go u gs are analytic functions, whereas g; U go U gz is not. 
Here u symbolizes the union of sets of pairs of numbers. The previously men- 
tioned identity theorem concerning analytic functions implies that F-branch points 
are isolated on the F-boundary of f. 

General Multifunctions.—Where does analyticity enter into the proposed defi- 
nition of a multifunction F? Only in one place, namely, in the proof that the F- 
branch points of any restriction of F are isolated. Consequently, general continuous 
multifunctions can be introduced just like analytic multifunctions. The class of 
continuous functions (with open connected domains) becomes the foundation of a 
superstructure of multifunctions if, in conditions I and II, the word “analytic’’ is 
replaced by “continuous.’’ Even the isolated character of the F-branch points 
of any continuous function SF can be guaranteed by the following special as- 
sumption: if two continuous functions S&F, f,; and fo, assume equal values at a 
point zo, common to their domains, and at all points of a sequence converging 
toward z, then f; and fz are identical in some circle about z. An example of such 
a continuous multifunction, which is real-valued and therefore nonanalytic, is the 
class of all pairs (z, arg z) forz # 0. Without the special assumption, a continuous 
multifunction may have branch lines; e.g., the class of all pairs (7 + yz, ny) for all 
real x and y and all integers n has the z-axis as a branch line. 

Clearly, the procedure just described can be extended in various ways: 

1. From the Argand plane to any more than one-dimensional Cartesian space 
S (and even to more general spaces including, e.g., the Riemann sphere). 

2. From the analytic or continuous functions to any given class § of functions 
whose domains are connected open subsets of S. 

3. From isolated branch points to branch sets of a dimension not exceeding a 
prescribed number k < n. 

4. From functions on S to mappings of S into another space S’. 

If the last generalization be here disregarded, the procedure yields multifunctions 
in S which locally belong to the given class ¥ and whose branch sets are at most k-dimen- 
sional. Special cases of such multifunctions have been studied (although not quite 
in this form) in the theory of magnetic plates, in Sommerfeld’s theory of black 
bodies,? and, of course, in differential geometry. 

Polyfunctions.—F¥or given S, §, and k, by the sum, F, + F», of two multifunctions, 
I mean the class of all pairs (z, w) such that there exist two numbers, w; and tw», 
satisfying the following conditions: 
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(z,wi) € FP), (2, We) € Fe, and w, + we = w. 


The domain of F,; + F» is the intersection of dom Ff; and dom F2, and therefore 
open but, even if /; and F» are single-valued, not necessarily connected. If f; and 
f, are two analytic functions which cannot be analytically continued beyond their 
domains, and if the intersection of these domains is not connected, then f; + fo is 
the union of a finite or denumerable number of analytic functions. Similarly, 
multiplication (defined in analogy to addition) is not closed in the realm of general 
multifunctions nor, in particular, for the classical analytic multifunctions or even 
the analytic functions. 

One-place substitution can be defined only for functions whose domains are 
classes of numbers.* By FF». (the result of substituting 2 into F;), I mean the 
class of all pairs (z, z,) for which a number 2 exists such that (z, 22) « Fz and (2, 2) € 
F,. If F, is the exponential function and F, the logarithmic multifunction, then 
FF» is the class of all pairs (z, e'°* *), that is, the identity function ¢, consisting? of all 
pairs (z,z). The class FF, of all pairs (z, log e*) is the union of denumerably many 
functions, « u («+ 2r7) u (« + 477) uU.... This union of different functions is 
also an object of the classical theory. Moreover, the latter includes functions with 
multiplicities, such as the exponential function on the Riemann surface of the loga- 
rithmic multifunction (which is. with denumerable multiplicity), the square on the 
surface of the square root (which is « with the multiplicity 2), ete. This is a com- 
plication of which the theory here outlined obviously is free. 


By an analytic poly-m-function (‘‘m’’ for “multi’’), I mean any element of the 
smallest trioperational algebra® generated by the class of all analytic multifunc- 
tions, that is, the smallest set of classes of ordered pairs of complex numbers which 
‘an be obtained from analytic multifunctions by the operations of addition, multi- 
plication, and substitution. This set includes, in particular, the vacuous class of 
pairs of numbers. 


'! Part of the work on this paper was done as part of Project DA-11-022-ORD-1494, sponsored 
by the Office of Ordnance Research. The ideas here outlined will be elaborated in the author’s 
forthcoming book, A New Approach to Mathematical Analysis and Its A pplication in Physical Science. 

2 Cf. Sommerfeld’s chap. xx, in Frank-Mises, Die Differentialgleichungen der Mechanik und 
Physik (2d ed.; Braunschweig 1935), Vol. 2. 

‘If the domain of a generalized “function” does not consist of numerically characterized 
elements, in other words, if with each element of an abstract or nonarithmetical class a number is 
paired (e.g., with each act of observation of a certain type, the numerical result of that act), then 
this pairing does not lend itself to any substitution of a numerically valued function, even though 
such pairings can be added and multiplied just like analytic functions. Since substitution is the 
functional operation par excellence, I regard the fact that pairings of numbers with abstract or 
as a symptom of the overemphasis on the 
variable quan- 


” 


nonnumerical elements have been called “functions 
algebraic and vectorial operations in current analysis. I have proposed the name 
tity” for such pairings, which include the pressure of a gas and the volume as a sphere, as well as 
functions, such as the logarithm of a positive number and the product of two numbers. Cf. ‘‘The 
Ideas of Variable and Function,” these PRocEEDINGs, 39, 956-961, 1953. These ideas are elab- 
orated in chap. vii of the second edition of the author’s book Calculus: A Modern Approach 
(Chicago, 1953) and in the forthcoming paper “Variables in Mathematics and in Natural Seience,” 
Brit. J. Phil. Sei. (1954). 

‘Cf. the publications quoted in n. 3. Strictly speaking, FF. is the restriction of « to the class 
of numbers +0. 

° Cf. “Tri-operational Algebra,’’ Rept. Math. Colloquium, 2d ser., 5-6, 3-10, 1945, and ‘‘General 
Algebra of Analysis,” ibid., 7, 46-60, 1946. 
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ON THE SYMMETRY PROPERTIES OF POWERS OF REPRESENTATIONS 
OF THE TWO-DIMENSIONAL UNIMODULAR UNITARY GROUP 


By Francis D. MurNAGHAN 
INSTITUTO TECNOLOGICO DE AERONAUTICA, SAO JOSE DOS CAMPOS, BRASIL 


Communicated July 2, 1954 


The 2-dimensional unitary unimodular group has one irreducible representation 
of each dimension, and we denote by |m}, where m is any nonnegative integer, the 
irreducible representation of dimension m + 1; thus {0} is the identity representa- , 
tion, {1} is the self-representation, {2} is the 3-dimensional representation, namely, 
the 3-dimensional rotation group, and so on. It is customary to give prominence to 
the 3-dimensional rotation group and to speak of the representations |m} of the 2- 


dimensional unitary unimodular group as spin-representations (7), where 7 = m/2, ( 


of the 3-dimensional rotation group, so that (1/2) denotes the 2-dimensional unitary 
unimodular group, (1) denotes the 3-dimensional rotation group, and, generally 


(j) denotes the (27 + 1)-dimensional spin-representation of the 3-dimensional rota- r 
tion group, j being either a nonnegative integer or half an odd positive integer. In P 


the question we propose to discuss here, the fact that the underlying group is 2- 
dimensional, and not of dimension greater than two, is an essential simplification, 
and so we shall use the notation |m}, rather than the more familiar notation (j), I 
but anyone who finds it more convenient to use the latter notation need only set 
{m} = (j), where j = m/2. It is well known that |m}2, the square of {m}, may be 


is the completely symmetric part, and {m} ® 1‘ is the completely alternating part, 
822 


written as the sum of a symmetric part and an alternating part; we associate the ( 
symmetric part with the partition 2, and the alternating part with the partition 1’, F 
of 2, and we denote these two constituents of {m}2 by {m} @ 2 and tm} @ 1, 
respectively. The dimensions of |m} @ 2 and {m} @ 1? are those of the repre- 
sentations {2} and {12} of the (m + 1)-dimensional linear group, namely, */s(m + 
1)(m + 2) and '/sm(m + 1), respectively; each of them is a representation, in ( 
general reducible, of the 2-dimensional unitary unimodular group, and we have the a 
¢ 
{m} @ 2 = {2m} + {2m — 4} +... 
Ci 
(the last term being {0} or {2} according as m is even or odd); T 
2 ‘ ‘ ‘ ¢ > fi 
{m} @ 12 = {m-1} @ 2 = {2m — 2} + {2m —6} 4+... U" 
(the last term being {2} or {0} according as m is even or odd). In the same way 
may be written as a linear combination of three representations each of H 
which is associated with one of the three partitions 3, 21, and 1° of 3 (each of these _ 
constituent representations appearing the number of times indicated by the dimen- (0 
sion of the corresponding irreducible representation of the symmetric group on tt 
three symbols). Thus, on denoting the constituent representations in question by rs 
{m\ @ 3,{m} @ 21, and {m} @ 1°, we have {m}* = ({m} @ 3) + 2({m} @ 21) + ° 
({m} @ 1°). {m} @ 3 is the symmetrized Kronecker cube of {m}, and {m} @ : 
is the completely alternating part of {m}%. Similarly, {m}* = ({m} @ 4) + , 
3({m} @ 31) + 2({m} @ 2?) + 3({m} @ 212) + ({m} @ 14), where, again, {m} @ 4 ‘i 
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of {m}4, and so on. We propose in the present note to furnish the analysis of 
| m} @ (A), where (A) = (Aj,A2,...) is any partition of an integer k > 2, and m is an 
integer > 1, for those values of m and k for which mk < 20, it being understood that 
10} @ (A) = {m} @ 0 = {0}, {m} @ 1 = {m}, and {1} @ (A) = {a}. 

Our calculations are based on the following two facts which seem to have escaped 
attention. First, we have the recurrence formula 


{ml} Ok = (jm —2} @k) + (im —- 1} Q@k—1)(jm+k—-1t - 
pm + k — 3) + ({m} @ k — 2), 


where a} 1b} is furnished by the familiar rule la} {pt = fa + b} ob finite Bi 28 as 
—- {|a a bt. For example, 


{3} @ 3 = {3} + ({4} + {o})({5} — {3}) + {3} = fo} + fot + {3}. 
Second, 
{m} @1* ={m—k+1} Ok. 


This is a generalization of the well-known result {m} @ 12 = {m — 1} @ 2. For 
example, 


{3} @ 18 = {1} @3 = {3}. 
In order to obtain {m} @ (k — 1, 1), we use the relation 
{m} @ (k — 1,1) = ({m} @k — 1){m} — {m} @k 


(which follows from the relation |k — 1, 1} = {k — 1}{1} — {k}). Similarly, in 
order to obtain |m} @® (k — 2, 2), we use the relation 


{ml} @® (k — 2,2) = ({ m} @k—- 2)({m!} @ 2) — ({m} @k) - 

im} ® (k — 1, 1) 
(which follows from the relation {k — 2,2} = {k — 2}{2} — {k} — {k — 1, 1}), 
and so on. In the case where k = 3 we have two alternative formulas for the 
calculation of {m} @ 21, corresponding to the two relations {2}{1} = {3} + }21}; 
(14f1} = {21} + {13}. Of these two formulas, the second is the simpler be- 
cause |m} @ 12 is simpler than {m} @ 2 and {m} @ 1% is simpler than {m} @ 3. 
Thus, to calculate {3} @ 21, we set 


{3} @ 21 = ({3} @ 12){3} — ({3} @ 1% = ({4} + {o}){3} — {3} = 
t7} + {5} + {3} + {1}. 


However, the availability of the two formulas furnishes a convenient check on our 
calculations; a second check is furnished by the fact that the dimension of {m} @ 
(\) is that of the representation {\} of the (m + 1)-dimensional linear group. If 
the partition (A) = (Ay, ..., Ap) of & contains p nonzero parts, this is (m + dy) !(m + 
he — 1)!... (m+ dr, — p + 1)! times the dimension of the representation (A) of 
the symmetric group on k symbols divided by m(m — 1)? ... (m — p + 2)?! 
((m — p + 1)!}? times k! (the divisor being m!k! if p = 1). In particular, the 
dimension of |m} @ kis (m+ k)! + milk. 

In listing our results (Table 1), we do not give {m} @ (A) if the number p of 
nonzero parts of (A) = (Ax, ..., Ap) is greater than m + 1, since {m} @ (d) vanishes 
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when this is the case. This is a generalization of the well-known result that 
i} = {1! @ (A) vanishes if (A) contains more than two nonzero parts. Thus, 
since the representation |), \o} of the 2-dimensional unimodular unitary group 
equals {\: — de}, the analysis of {1}* is of the form {k} + ajk — 2} + afk — 
4 + ..., where c, = k — 1 is the dimension of the representation (k — 1, 1) of the 
symmetric group on k symbols, c. = k(k — 3)/2 is the dimension of the representa- 
tion (k — 2, 2) of this symmetric group, and so on. We do not list }m} ® (A) when 
the number p of nonzero parts of (A) is m + 1, for this is the same as |m} @ (Q’) 
where (A’) = (Ai — Amst, A2 — Amoi, ---, Am — Amy) 18 the partition, with less than 
m + 1 nonzero parts, of k — (m + 1)\n41. For example, {2} @ 1° = {2} @0 = 
fol. (2! @ 321 = {2} @ 21; {3} @ 31% = {3} @ 2 = {6} + {2}; and so on. 
We conclude with the remark that the relation {m} @ 1° = {m —k + 1! @ k, 
together with the relation {m} @ 1* = {m} @ 1”~**! (which is an immediate con- 
sequence of the fact that {m} is of dimension m + 1), imply the relation {m} @ 
k = {kt @ m, which is known as Hermite’s Law of Reciprocity. We have given 
previously! another proof, based on the relation ({m — 1} @ k){k — 1} = ({m} @ 
k — 1)}m — 1}, of this law. 


I have to thank Professor E. P. Wigner for calling to my attention the importance 
in spectroscopy of the problem here discussed. 


1F. D. Murnaghan, these PRocEEDINGs, 37, 439-441, 1951; Anazs. Acad. Brasil. cienc., 23, 347— 
368, 1951. 


A VON STERNECK ARITHMETICAL FUNCTION AND RESTRICTED 
PARTITIONS WITH RESPECT TO A MODULUS 


By C. A. Nicot anp H. 8. VANDIVER 
UNIVERSITY OF TEXAS 
Communicated June 28, 1954 


Introduction.—V on Sterneck! in 1902 employed the number 


y(n) n 

Pree g(n/(k, n))" (k, mn)’ 
where y(n) is the Euler indicator (or totient) and y(n) is the Mébius number. Also, 
(k, n) is the greatest common divisor of k and n; k and n integers with k > 0,n > 0 
and where (0, n) = n. (Further, unless otherwise stated, all small italic letters 
used will denote nonnegative rational integers.) He encountered (1) in investi- 
gating the number of different ways an integer s may be expressed as the sum of ¢ 
integers if all the summands involved are reduced to their least residues modulo 
m>1. He proved that 


> &(k, d) = i if (k, n) = n, es 


din 0 otherwise. 


He derived? also the relation, if d|n, 


Dx, d) = 


nifd = 1, 
0 otherwise, 


(1b) 
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and noted the multiplicative property 








b(k, m)d(k, n) = &(k, mn); 





(m,n) = 1. 





[Note that (1) reduces to g(n) when (k, n) = n and to u(n) when (k, n) = 1.] In 
the opinion of the authors, (1) will be very much employed by number theorists in 
the future and will eventually become classic in the literature of number theory; 
hence it will be convenient to give it a name. We shall call it the ‘Von Sterneck 







number.”"* 
In 1920 Ramanujan‘ considered the number 


CLh= > a’. a = erin (2) 


re a aed 











(Throughout this paper, a will have this meaning.) Such exponential expressions 


are called Ramanujan sums. He showed that 


Clk) = >> pwln/d)d. (3) 


d\k, d\n 








He also employed C,,(/) in certain infinite series. 
In 1936 Hélder® obtained the following fundamental result: 


Ci,(k) = ®(k, n). 











(4) 


He gave two proofs of (4) and indicated another. He also obtained results con- 
cerning infinite series involving C,,(k). 

Nicol® obtained a number of relations involving (1). In particular he related 
@(k, n) to the theory of the restricted partitions of the natural numbers which have 








the two generating functions 





ID (Ch: ote: 2’). (5) 


s=1 














Apostol and Anderson’ defined a number which is a generalization of (1) and 
treated infinite series involving it. It is possible that some of the theorems of the 
present paper concerning (1) may be extended to the generalized number, but we 
have not yet examined this possibility. 

Gagliardo® gave another proof of Hélder’s theorem. 

In the present paper we obtain a number of properties of (1) which involve the 
&’s only, or the ®’s and roots of unity, the proofs all depending on the use of (4), 
excepting Theorem VIII. We then follow with an application of the same methods 
to a problem involving the partition of integers with respect to a modulus (Theorem 
VII). 

This article is introductory to another article which we hope to publish in which 
arithmetical and analytic methods are applied to obtaining properties of the co- 
efficients of (5) involving, usually, the Von Sterneck numbers. 

Properties of the &-F unction.—We shall have need in several places in this article 
of a result which we shall state as Lemma 1. 

Lemma 1. Jf din, c = ¢g(n) and 


7, ee (6) 


are the positive integers less than n and prime to n, then the number of the integers im 
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the set (6) which are congruent to t modulo (n/d), with (t, n/d) = 1, is (d/b)¢e(b), where 
bis the greatest divisor of d such that (b, n/d) = 1. 
Consider the set 


In - 
: in eg, (6a) 
d 
ry ; 
eck where j = 0,1,...,d¢@ — 1. The integers (6a) are prime to n/d but not necessarily 
prime to n. We shall select those which also have the latter property. Suppose 
that b > 1, where b is defined in the lemma. Thed integers 0, 1,..., d — | may be 
(2) written as 
m + hb, (7) 
ons : 
wherem =0,1,..., b—landh=0,1,..., (d/b) — 1. If we use these values 
of j in (6a) and reduce the resulting expressions modulo b, we obtain the set of 
(3) residues 
n _ 
t+-m (7a) 
d 
(4) corresponding to each value of h. Since (b, n/d) = 1, this set modulo b reduces to a 
cyclic permutation of 0, 1,...,6 — 1. By definition, d = d,b, where each prime 
on- factor of d; dividesn/d. To obtain the integers in (6a) which are prime tod and con- 
sequently prime to n, it is sufficient to select those which are prime to 6, since they 
ted all prime to n/d and hence to d;. This subset of (6a) reduces modulo 6 to the 
ave g(b) integers in the set 0, 1, ...,6 — 1 which are prime to b. Now there are d/b 
values of b corresponding to each set (7a). Hence the number of integers prime 
ton is (d/b) ¢(b), as stated in the lemma. If b = 1, each prime factor of d occurs 
(9) in n/d, so that each element in (6) is prime to n. There are (d/1)¢(1) or d such 
numbers, and the lemma also follows for this case. 
ind As one application of this lemma, we shall give another® proof of Hélder’s theorem. 
the Suppose that a, 7, and n are positive integers, and let 
we ? Oni/y 
S= DI a”, where a = e""", 
(r, n) = 1 
sii Also, let (a, n) = d and kd = a, so that (k, n/d) = 1. Denote the positive integers 
4) less than n and relatively prime to n by 
ods 7; ee where c = g(n), 
em ie pl ‘ : 
and the positive integers less than n/d and relatively prime to n/d by 
ich n 
‘iia tte. oe where v = ol" ; 
( 
‘cle Then the exponents appearing in S are 
kar, kdre, + 4x8" kdr,. (8) 
Consider the set 
n 
t+J-, 
d 
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where 2 = 1,2,...,vandj7 = 0,1,...,d— 1. We then employ Lemma | and 
note that the number of integers prime to n for each ¢ in (9) is the same. Alto- 
gether, there are g(n) positive integers less than n and prime to n. Consequently, 
(8) may be replaced by kdt, kdto, . . . , kdt, repeated g(n)/g(n/d) times. We shall 
assume the well-known relation” 









(9a) 





(r, n/d) = d 






Since (k, n/d) = 1, we have 


Also, we note, 








> af? = _ (a*)"” a ( 11 ) 


(r,n) = 1 (r,n) = 1 









Then, from (10) and (11) and noting the remark just above (9a), we have 


g(n/d) (r, n/d) = 1 g(n/d) d 






which proves (4). 

We shall now prove (la). (Von Sterneck’s proof was more complicated.) Sup- 
pose d divides n and n > 0. Consider the set of positive integers less than 
n + 1 whose greatest common divisor with n is d. These may be written as da, 
day, ..., da)», wherev = g(n/d) and a, dz, ... , a, are the g(n/d) integers less than 
n/d and prime to n/d. Then we have, by (4). 


> &(k, d) a 2 =. og Ok 


d\n d|n (r,d) = 1 















since a“ is a dth root of unity. In the last relation, as d ranges over each divisor 
of n, we obtain the kth powers of all nth roots of unity, so that we may write 


> a I *) =lta°*t+a*t+...¢ a", 
€ 


d\n 
and this sum is zero unless a* = 1, in which case it is n, which gives (1a). 
If p is a prime such that p = 1 + tn, then the decomposition of (p) into prime 
ideal factors shows that one of them is (g' — 1, p), where g is a primitive root of p. 
Hence a = g‘(mod p), and, from (4), we have 
~ gg” = O(k, n) (mod p), 


(r,n) = 1 
or 


yg” = &(k, n) (mod p). (12) 


(r,n) = 1 





This proves a result due to Moller."! 















ind 
ito- 
tly, 
all 


Ya) 


10) 


sup- 
han 
da, 
han 


‘1s0r 


rime 


D. 
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Suppose & > 0, and consider the number 


n 


Ld Hk, s). 


s=1ad|/s 


From relation (1a) we note that the inner sum is zero unless ¢k, s) = s when it is 
s. Thus this number is the sum of the divisors of k less than or equal to n, which 
we shall denote by o(k, n). 

However, we also have 


~ ¥ ee = & "| P(I:, 8), 
s= 1 ; a 


where |n/s] is the largest integer in n/s. Therefore, we may state. 
THEOREM I. We have the relation (Paper N, [9]; note also the two equations 


following tt) 


> Bi sh ence. ds * @ 


s=1 


where k > 0 and o(k, n) is the sum of the divisors‘of k less than or equal to n. 
Consider the number 


(n/d) — 1 


>~ a“ &(dt, n). 


t=0 
Using the Hélder theorem, we have 


(n/d) — 1 (n/d) — 1 


=. «* @(dt, n) = > 2 oom z Pr a 


t=0 t=0 (r,n) = 1 
which may be written as 


> a Me ai’) ob ete Ae ee am a!” d@)—1)(7—s)) 


(r,n) = 1 


(r—s) 


The sum of the terms in the parenthesis is zero unless a = 1,in which case it 
isn/d. Since ais a primitive nth root of unity, a’ = | only if r=s mod (n/d) 
and the number of nonzero terms in the sum taken over the positive integers less 
than n and prime to n reduces to the number of positive integers less than n and 
prime to it which are congruent to s modulo (n/d), multiplied by n/d. If (s,n/d) # 
1, then nor exists. If (s, n/d) = 1, then, by Lemma 1, the number of integers less 
than n and prime to n which are congruent to s modulo (n/d) is (d/b)¢(b), where b 
is the largest divisor of d such that (b, n/d) = 1. This gives 
TuroreM II. Jf n> 0 and d\n, then 
“— J (n/b) e(b) tf (s, n/d) = 1, 


—tds 3, 
t=0 ” as laa ca 10 otherwise, 


r—s) 


(14) 


where a = e**"/" and b is the largest divisor of d such that (b, n/d) = 1. 
Suppose that p is a prime and g is a primitive root modulo p. Then we shall 


show that, modulo p, 
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p—-1 oe ~ 
ice taal a ep ed, (15) 
a —e oe (0 otherwise. (16) 


In relation (14), if we choose d = 1, thenb = 1. Also select n = p — 1, where p is 
an odd prime with (s,p — 1) = 1. Then (14) becomes 
e 


p-l 


L a H(t, p-lhb=p-1. (17) 


t= 
If p is a prime ideal divisor of (p), then we may write a = g (mod p) for g a primi- 
tive root of p. Then (17) becomes 


Pp 


1 
g°®(t, p — 1) = —1 (mod p) 
1 


iM! 


t 
or 


=  —1 (mod p), where (s, p — 1) = 1, 


which is relation (15). In a similar way, relation (14), for (s, p — 1) ¥ 1, becomes 
relation (16). Relations (15) and (16) may be combined to give the following re- 
sults: 

The only integral roots x of the congruence® 


p-—l 
> x'&(t, p — 1) = —1 (mod p) (18) 
t= 1 

are the g (p — 1) incongruent primitive roots of p, and the only integral roots y of 


the congruence 
a3 
Dd y'&(t, p — 1) =0 (mod p) (19) 
t=1 


are the integers in the set 0, 1,...,  — 1 which are not primitive roots of p. 
THeorEM III. Jfn> 1, then 


Sethe sic ane. (20) 


= 1] 


Proof: From (4) we have 
dX sh(s,n)= Ds Ye a= YP (a + 2a" +... + na). (20a) 
1 


s=1 s=1(r,n) = 1 (r,n) = 


As is well known, however, 


; . l 
a’ +207 +...+ na” = n( : i aa ). 
ll sna 


Consider now the function 


F.@)= WW @— a’). 


(r,n) = 1 
Differentiating and letting x = 1, we have 


FQ) sia 
F,(1) (r,n) = 1(a’ = 1) 
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D) Holder!” has shown that F’,,’(1) is equal to (1/2) ¢(n) when n > 1 and not a power of 
6) a prime, and F’,’(1) is equal to p*(p — 1) whenn = p*. 

Also, as is well known, 


is 
F,(1) = J 1 if n contains two prime factors, 
\p otherwise. 
7) ; 
Therefore, we have, from (20a) and (20b), 
\i- n 
no(n 
dX s&(s,n) = — o ) + ng(n), 
_ 2 
which establishes the theorem. 
If n > 0 and 2) n, consider the expression 
: S=1+a"+a"+...+a"™. 21) 
Since a’ is a (n/2)th root of unity, it follows that the above may be written, using 
as (4), as 
a 
P o(1 %) (22) 
Uy ¢ P 22 
d|n/2 2d 
Now consider 
8) 
d " . 
ot d|n d 
This may be written as S — aS, using (21), since 
9) a ta* +... ta" = aS. (24) 
We also have, from (21), that S = 0; except when k = n/2, when it is n/2. Then 
S — a‘S reduces to 0, except when k = n/2, when it reduces to n, and we have® 
Theorem IV. 
0) THEOREM IV. If n is even and n and k are positive, then 
- n lnifk = n/2, 
> (—1)% a k, = 4 bv 
djn d 0 otherwise. 
) Let w = &"/""/ and d|n andd <n. Then, by (4), 
n n n n/d e 
dY Fs, = Zz y @ = dz a’, 
s=]1 d (r.n/d) =1 s=1 a’ =] 
after we set s = k(n/d) + 8’, withO < k <dand0 < s’< n/d. Hence the second 
n 
factor on the right is zero. When d = n, the sum >> ®(s, 1) = n._ This gives 
s=1 
(1b). Von Sterneck’s proof was more complicated. 
Consider the cyclotomic function 
G(x) = II (x — a’) (25) 
(r,n) = 1 
Ib) and the problem of finding a convenient way of determining its coefficients when 
Written in polynomial form. First, the coefficients are, in absolute value, the ele- 
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mentary symmetric functions which are formed by the distinct g(n) primitive roots 
of unity. In view of (4), the sums of like powers of these roots are all known in 
terms of the Von Sterneck numbers. 

Suppose G(x) = a + pao '+...+pands;= >> a” = (i, n), where i = 


(r,n) = 1 
1,2,...,¢; ¢ = g(n). Consider the Newton formulas as applied to G(x). 

















8 + 18-1 +... + Dp-1h + ip, = Ofortz = 1,...,€¢. 


Let us assume that p, (a < c) has the property that it can be expressed as a poly- 
nomial in the Von Sterneck numbers with rational coefficients. Using the relation 


..- op 1) pa+1 = 0, 







Sati + DiSa + . 













it follows immediately by induction that p,, has the same property. Since this is 

obvious for p:, it then follows that each of the coefficients in the cyclotomic polynomial 
25) may be expressed as a polynomial in the Von Sterneck numbers with rational 
coefficients. 

We shall now consider two inversion formulas. 

THeoreM V. /f g,(s) = >> f(d)®(s, n/d), where f(n) is an arithmetic function, 


din 

















then 


n 


l 
S(n) = PD Yn(8). 


Ns =1 


Proof: rom the hypothesis we have 


Lae = 2-2 j(d)e(s,") =VI@LX #(s.i) 
l s= 1 ’ 


s=1din d\n 


From relation (1b) the only term in the right-hand member which is not zero occurs 
when d = n, and thus the theorem follows. 


TuroreM VI. If g,(d) = >> f(k) ®(k, d), where f(k) is an arithmetic function, then ; 
k t 


=] 


! 
fin) = 1 gad). : 


The proof of this theorem is the same as in Theorem V, except that relation (1a) is 
used. 


Restricted Partitions Modulo n.—Consider the expression 
n- 1 
= ae (t, n), 
t=0 


where m is a positive integer and a is a primitive nth root of unity. Then, from 
(4), we have 


n— 1 


n-—1l1 
=, a "@(t, n) = z. a" > ra = ai” = 
t=0 


t=0 (m,n) = 1 (tam, n) = 1 


— .. 2 Cie... Ser 


(m,n) = 1 (rm, n) = 1 
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where y = (4%) + 22 +...+ 2%, — 8). Thesum1 + a” +...+a”~”” is zero 
Y 


unless y = 0 (mod n) or a, +... + &m=s (mod n), in which case it isn. Then, if 
P,,(m, 8) is the number of solutions of the congruence x; + ... + 2m =s (mod n), 
in which the order of the terms is taken account of and (x;,n) = 1,7 = 1,2,..., m, 


we have 


n—-1 


> a “b(t, n) = nP,(m, s). (26) 


t=0 


Suppose that d is a positive integer which divides n, and let Sq be the set of positive 
integers less than n + 1 whose greatest common divisor with n isd. Then we may 
write 


n— 1 


y a “a(t, n) = OD a “O(a, n), 
= 0 


t din Sd 


where the inner sum is taken over the set Sz. However, for a positive integer D 
in the set Sg, the number (D, n) remains constant, and we may write 


n=— tI 


> a “b"(t,n) = > od, n) > a 
d\n Sd 


t= 0 


However, from (4), the inner sum may be replaced by ®(s, n/d). Therefore we 
have Theorem VIII. 
THEOREM VII. 


n 
>» &"(d, n)® (s, ) = nP,,(m, s). (27) 

d\n d 
where P,,(m, s) 1s the number of solutions of the congruence 2 + % +... + lyn = 
s (mod n), two solutions (21, V2, ... , Xm) and (ay', x2’, . . . , Xm’) being regarded as the 
same if and only if x; = 2;' (mod n), where (a;, n) = 1 fori = 1, 2,..., mand 


m>0O,n> 1. 

For small values of m, at least, the right-hand member of (27) may be evaluated 
independently of the use of the ¢-number. For m = | in relation (26) we have a 
special case of Theorem II, when d = 1 in the latter statement. Form = 2ands = 
0 we determine explicitly the value of P,,(m, s) as g(n). This® gives (26) for 
m=2,s=0. Ifwetakem = 3,s = 0, then we have to find the x,’s with (x;, n) = 
1,7 = 1, 2, 3, such that x; + x2 + x3 = 0 (mod n), which is equivalent to x,’ + 
v2’ + 1=0 (mod n), (a1’,n) = 1 and (22’,n) = 1, so we have the problem of finding 
the number of sets of consecutive positive integers < n and prime to n. 

Schemmel!’ stated that n II (1 — (2/p)) is the number of pairs of consecutive 

p 


positive integers less than n which are each prime to n. This is equivalent to the 
number of positive integers less than n which are relatively prime to n and less than 
nand which remain prime to n when unity is added to each of them. In this con- 
nection we have Theorem VIII. 


lHrorEM VIII. Jf nis an integer greater than one, then 
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2 
> &(d, n) = nil (1 = ), 
d\n Pp Pp 
where the product is taken over all prime divisors of n. 


Proof: From relation (1), we have 


= u(n/(d, n)) 
» &(d, n) = y(n) > o(n/(d, n)) 


Since d\n, then (d, n) = d, and thus 





ud) 
din din g(d) 


> &(d, n) = g(n) > 


Since u(d) and ¢g(d) are multiplicative, then the number >> &(d, n) is also multi- 
din 


plicative. Hence it is sufficient to consider n = p* where p is a prime and Gis a posi- 
tive integer. Then 


l 2 
&(d, p’) = of ry(i _- }= (1 = ), 
Pm (d, Pp) p—l I Pp 


which proves the theorem. 
Von Sterneck!* found an expression involving (1) for the number of representa- 


tions of s as the residue modulo n of a sum of m elements chosen from 0, 1, ..., 
n — 1, repetitions allowed but with order of the summands disregarded. 
Whiteman determined the number of sets of integers 4, Ue, ..., Uss 1 Su; < 
m, — 1 for which a"... a,"* = 1, where a; = erm = 1,2,..., 8. If mis the 
L.C.M. of these m,’s, this is the same as finding the number of representations 
V1, Ho, . . y Ly IN MX, + Mae +... + M2, =0(mod m). This is equivalent to the 


problem! previously discussed, except that our representations are restricted to 
certain types of repetitions. 


1 Sitzber. Akad. Wiss. Wien, Math.-Naturw. Klasse, 111 (Abt. Ila), 1567-1601, 1902. Cf. 
also Bachmann, Niedere Zahlentheorie (Leipzig: B. G. Teubner, 1909), 2, 230-232. 

2 Bachmann, op. cit., p. 232. 

3 We hope that other mathematicians will check our conclusion that Von Sterneck was the 
first investigator to isolate this number and obtain a few of its properties. 
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tent these sums had been considered before; ef. Hardy and Wright, The Theory of Numbers (Oxford, 
1938), p. 62. 

5 Prace Mat.-Fiz., 43, 13-23, 1936. 

6 These PROCEEDINGS, 39, 963-968, 1953. This paper will be referred to as ‘Paper N.’’ The 
present article is the result of joint work of Nicol and Vandiver, except as follows: Theorem I is 
due to Nicol (Paper N, [9]), who found a proof less direct than the one given here. He discovered 
Theorem IV and found a different proof. He obtained (26) for m = 2, s = 0 and obtained an- 
other type of proof for that case. The same remark applies to (18). All the results ascribed to 
Nicol were stated by him without proof in the reference given at the beginning of this footnote. 

7D. R. Anderson and T. M. Apostol, Duke Math. J., 20, 211-216, 1953. 
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9 Hélder gave two proofs of (4) and indicated another, and in addition we have another proof, 
as noted, by Gagliardo. If the reader desires to examine a proof of (4), our proof may be a con- 
venience, as the articles of Hélder and Gagliardo were published in journals which are not readily 
accessible, at least in the United States. 
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« Cf, . Landau, Handbuch der Lehre von der Verterlung der Primzahlen, 2, 572-573, 1909. 


1 Am. Math. Monthly, 59, 228, 1953. 
12 Cp. cit., p. 15. 


13 J. f. Math., 70, 191-192, 1869. This and other of Schemmel’s results were proved by Bach- 


mann, Niedere Zahlentheorie (Leipzig: B.G. Teubner, 1902), 1, 91-94. 
14 Sitzber. Akad. Wiss. Wien, Math.-naturw. Klasse, 114 (Abt. IIa), 711-730, 1905. 
4 Trans. Am. Math. Soc., '74, 87, relation (4.11), 1953. 


SPHERICAL SPINORS IN A EUCLIDEAN 4-SPACE 


By A. Pats 
ti INSTITUTE FOR ADVANCED STUDY, PRINCETON, NEW JERSEY 
- Communicated by J. Robert Oppenheimer, June 29, 1954 
1. In a foregoing paper! some qualitative properties of baryons (nucleons and 
hyperons) in interaction with bosons were discussed in terms of a model in which 
these particles are characterized by quantum numbers corresponding to certain 
representations of the group of real 4-dimensional rotations with det. = +1. 
the purpose of this note to append these considerations by explicitly constructing 
ta- the spinors which, in this model, describe the baryon states. This is a simple 
oat : matter, and it is hoped that it may elucidate the group theoretical argument given 
previously. 
< : It was important for the previous discussion! of charge-independence to note the 
the i connection between 3- and 4-dimensional rotations. This connection has long been 
ns ' known. In physical terms it is exhibited by the classical result? that the motion of 
the a (3-dimensional) spherical top can be represented by the motion” of a 4-dimen- 
to sional spherical pendulum, i.e., a point moving over a 3-sphere in 4-space. 
quantum-mechanically spherical top eigenfunctions can be written, in terms of 
Ct appropriately chosen variables, as eigenfunctions of a particle with intrinsic spin 
zero constrained to a 3-sphere, i.e., as 4-dimensional spherical harmonics.* 
4 all harmonics come into play, however; see the end of sec. 2.) Our present problem 
the { consists in constructing the eigenfunctions (spherical spinors) for a spin-!/2 particle 
| subject to the same constraint. 
oa In Section 2 the 4-dimensional spherical harmonics are expressed in a form some- 
what different than usual,‘ so as to bring in evidence formal connections with the 
Legendre polynomials.’ Recursion relations are obtained which closely resemble 
The those holding for associated Legendre functions and by means of which it is easy to 
1? find the spherical spinors (sec. 3). Some features of the spinors in Euclidean 4-space 
pi have been studied by Schrédinger* in a little-known paper on the Dirac equation in 
rye a closed spherical universe. The procedure followed here is a different and rather 
‘ simpler one. 
Thus the main interest of this note is one of method. The results may perhaps 
- be of use for field theoretical calculations, such as, for example, the relativistic two- 
eee: body problem. 
dily 2. In a Euclidean 4-space (variables x, ..., xs) we introduce the ‘angular mo- 
mentum”’ operator 
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ra) ) 
Ky = -1 (x — 2% , ih = 1,...,4, 
OX, Ox; 


KyKsy + Ko3Ky4 + Ky Ky = 0. 
















which satisfies 








The spherical harmonics are solutions of? 

















[1/> Kix? — A|Y = 0. 
It is well known that \ = j(j + 2),7 = 0,1,2,.... With Ag = '/2 (Kg, + 


K.a), a, B, y = 1, 2, 3 eyel., it follows from relation (1) that K+? = K~-? = K?, so 
that (2) can be written as 


[K? —k(k+1)] Y = 0, k = j/2 = 0, '/2,1, .... (3) 
Put 
x = rsin 6/2 cos ¢, 2 = rsin 6/2 sin ¢, 0<d0<7 
t3 = rcos 0/2 cosy, x4 = rcos0/2sinn, O< @ < 27; OX < 2m. (A) 


We have 
day day daz day = r* dr dQ, dQ = '/4 sin 6 dé dd dn. (5) 


Equation (3) becomes 


o? a) 1 oO? I 0? 
otg 8 eiaces iors ° —~— + k(k+1)|Y =0. 
E an 06 ” 2(1 — cos 6) 0¢? = 2(1 + cos 6) On? aE | 
Put cos 6 = x, and write the harmonics normalized over the unit sphere as 


1 : 
, ON. ~1/2¢ . (m + n)o + i(m — n) 
} aunt g, n) = 2N kmn femn(®) 5 eo - ie die (6) 
ray 


where Nim,’ is a normalization constant (see Eq. [16] below). Then 


]2 ] 2 2 ») ; if 
i i gg ee ee 1) fama =o @ 
dx? dx 1 — 2? 
Here —1 < x < landm,n = —k, —k + 1..., + k, corresponding to the 


(2k + 1)2-fold degeneracy of the eigenvalue. m, n are (half-)integer when k is 
(half-)integer. For m,n > 0 the (unnormalized) solutions (finite at 2 = +1) canbe 
written as follows: 


a) k integer: 


a Lee l at x hiatal: ENE 8) 
| a ( 1) 1 ry DpDrlk (2); ( 
xv 


b) k half-integer: 














S. 


30 


5) 


he 


be 
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: 1 =a 1/o(m + n) d 
Tran ( ) l+2z (1 + 2) A. 
(P,-1/,(X) ok Py41)(%)); 
_ 
Dn = (1+2)"——» (Do = 1). (9) 
dx 
As D,,D, = D,Dm, one has 
Sem ag —_ (10) 


a symmetry property inherent in equation (7). Thus the fim, are a kind of “twice- 
associated” Legendre polynomials; note that {.mo = Pm. One has the identities 


1/ Hime —'/or 
(1 — 2?) as + (mz + n)(1 — 2*)” “Sima = — Se, mt 1, 0 (11) 
dx 


\/o df, mn 


1 — 2? 
( sis dix 


— (mx + n)(1 — 2?) “femn = (k + m)\(k — m + Df, m1, ny 
(12) 


ti, m+1, n i 2(ma + n)(1 or 22)~ "fy. m, n + 


(k -+- m)(k —m + I) fi, ml1,n = 0. (13) 


Aecording to equation (7), one has 


So) fen(Drrale) de= 0, kk, (14) 
while 
Nemn = Sai) Pemn(&) dx (15) 
is given by 


Z k m)! (k n 
News = Een ete (16) 
2k + 1 (k —m)!(k — n)! 
valid for integer as well as half-integer k. 
We have now to consider the cases: m or n, or both m,n < 0. From the dif- 
ferential equation (7) it follows that fy, —m, n(x) and fx, m, _ n(x) are proportional to 


fk, m, n(—x) and that fy, _m, _n(x) is proportional to fx, m, x(x). We choose the 


proportionality factors so that the following conditions are true for all m and n: 
(1) the symmetry relation (10); (2) the identities (11)—(13); (3) the expression (16) 
for the norm as defined by equation (15). This is achieved by putting® 


k - 

fe, —m, n(X) = —1)" a ry ~ m;, a(— 
k- 

Ses m, ~a(2) = (-— 1) ae ra > Sima (— x) 


m+n (k — mk — n)! 


(k + mk + n)!° Sic, m, n(&)- 


he. —m, —n (x) == (+1 ) 
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The spherical harmonics Ym. satisfy 


ms | | | ?, n) Vermin’, ?,; n) dQ = Oxx’Omm'Onn’, 
Y; eri (i) ?, n) =(— re, mM, ght d, n), Vimn(9, ?, —n) = Yinm(9, q, n), 


Ys RS OO) « Ca) Sy |) {OO (17) 
We tabulate the first few fimn: 
fo, oo = 1, fi lgxtte Via gn (1 = a)’, 
fi, 0,0 = @q, fi, 1o= —(1 oa x?) hi, ti=l— 2; 


= ite =e o.0—'/s 1/; 
fs A VS Z : (1 = r) ; (1 + 32), far, thi —3.2 (1 + eyed ~ae x), 
. Poe ee 3/o 

for, 3/2, 3/2 = 32 A (1 — x) ‘ 

We recall finally that equation (3) is essentially the eigenvalue problem for the 
spherical top. However, a one-valuedness condition on the eigenfunctions (they 
should be unchanged if 6 > @ + 27) admits integer / only.° In our present case 
there is no such restriction; all k-values are allowed. 

3. Now we come to the spinors in 4-space. Let ¥,(a, ...,%4),@ = 1, ...,4 be 
a general spinor field in this space. We wish to develop y, with respect to its angu- 
lar dependence, in a complete set of orthonormal spinors. For this purpose we in- 
troduce the “total angular momentum”’ operator I,,: 

Tin = Kip + Sir, 
Six = —(t/4) [6:, Bel, BiB + BrBi = ix. 


Here S,, is the spin -!/2 tensor expressed in terms of the (4 X 4) Dirac matrices 
8; We find the desired complete set if we have all solutions of 


Hy = ny, H = V/oS in K ix. (18) 
Now one shows that 


H(H = 2) = V/oK ix”, 


—eor eam > => 


(the analogue of L S (L S + 1) = L?in3 dimensions), and hence, as also [H, Kim] 
= 0, 


n(n + 2) = 4k(k + 1), so n = 2k or —2k —2. 
Defining J.~ = 1/2(0g, + Jaa), we find 


T *? = '/.K x? + /4(1 + Bs) (CH + 3/2), Bs= — BiB 2838s. 
I +2 and / ~2 commute with H as well as with 8;. Putting J *2 = i*(i* + 1), one 
has 
f= 1, 9 = 2k: t= k+l/, i- =k; 
feel I n = —2k — 2: it = k—'/y, i =e 
6B = —l, n = 2k it = k, mg 1/9; 
6B=—-1l, n= —2k—2: it=k, i-=k—%. (19) 
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Relations (19) comprise the semispinor representations used in the previous paper 
+ *4 wal 1 


which were there characterized by (2+, 7~)*, |¢+ — a-| = '/2. To solve equation 
(18), we work conveniently in a special representation: 8; = pti, 7 = 1, 2, 3, 
8, = ps, hence B; = p;, in which equation (18) becomes 


('/2( + ps) rk + 1/o(1 — ps) T K ly = ny. (20) 


ft 4 
*e™ \o - 1)’ 


where / is a two-by-two unit matrix. Thus equation (20) decomposes as follows: 


Here p; is a diagonal matrix: 


(7 K *)pt = nyt, (21*) 
(1 Ko = We, (217) 


where the superscript + or — on y denotes that we have the eigenfunctions belong- 
ingtoB; = +1. wt and y~ have the following general structure: 


VY 0 

a Yo oe 0 
libs 07 * V3 
0 WV 


instead of which we shall write, in obvious notation, 


pee A ae V1 
. ies) - 2% 7 


Introducing the angular variables of section 2, one finds that K ~(6,¢,7) = AK *(@, 
¢, —n). Hence it suffices to solve equation (21+) to have all solutions. Equation 
(21+) can be written as 


i ‘ re) l / 
i [nisin (@ + n) — 72¢08 (+01. + saa } 71 cos (@ + n) + rosin (od + n)} % 


a) fa) - : ) 
vi a + '/,) cotg O(71 cos (6 + 9) + 72 8iIn (@ + n)) —T3y X 
(2 a >) |v mer 
a6 | On iin: 
Put 
y+ = U-1 oe i a ; 
¥ X; / exp > (@ + ») 


) 


x = aot es + er ee aad (x = COS 0). 
g(x) 
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840 
Then 
1/ l =i/, , 
k — x)” - + {(m — '/2)x + n}(1 — 2%) ‘ f= — ['Yeln + 1) + m)g, 
AL 


oe _1/, , 
q ~ 2)” : — {(m + !/)a + n}(1 — 2?) ‘i g = ['/o(n + 1) — mf. 
ax 


These equations are easily solved with the help of identities (11) and (12). Specify 
v+ fully as Y+i4. i, m,n» The final normalized results are 


L (Vk + m+ 2 Ve, min)! 
V2k +1 Vk — m+ 1/4 reve ji 


+ es 
v k+%/a,k, mn —~ 


yt = | . _— Vk - m + "/y Yx, m—'/o,n P 
k—'/2,.k, mn ~~ —— , ° 
V 2h + l Vk os m 1/4 ) k, m+'/o. n 


Examples: 


1 /1\t i nr | 
Do nome= —Z(g)s Mwe-mere()i @ 
T “ Tv “ 


2) yr ] — cos o/2e"™ , y 
“ O tf, Othe. ™ gi cz ‘ . ; ’ 0, 1/2, 0, —1/2 = 
Q2r\ V3 sin 6/2-e* in se 


1 (he sin =. (23) 


2a \cos 6/2+e"” 


The pair of constant spinors given by equations (22) corresponds to the nucleon 
doublet denoted by (!/2, 0)+ previously; the pair (23) was called (0, '/2)* in the 
foregoing! and corresponds to the (A°, A+)-doublet; etc. Note that /stWtis, i, m,n 
= my" s+, 1-. m. ny La ts, t-. ms nw = TW t+ t-. ms ny 80 that m, n are just the quan- 
tum numbers which were called! 73+, 73~, respectively. 

The solutions of equation (21~) can now be written down immediately, using 
relations (17) and (19): 


y as ee ] Vik + n we VW, Fy, m, n—'/2 
k, k+1/o, m,n — ae ape Ge Ui <r ’ 
k+ V/ 2k + 1 Vk iy + 1/, Ke m, n+ 1/2 


l ee PE TE ey 
V 2k + 1 Vk +n i VW, Yr, m. 2 +/s 


v k,k—'emn = 


Finally, one may ask how to set up a spinor calculus for the case of a Euclidean 
4-space, along lines similar as in the space-time description with its indefinite met- 
ric. We shall briefly indicate how this is done: 

a) In space-time one proceeds like this:!! call & (a = 1, 2) an “undotted co- 
variant spinor” if it transforms as £’, = > sA ats, det Aaal|| = 1, and call 94(8 = 
1, 2) a “dotted covariant spinor’’ if it transforms as &* does. (Note: £&, ng are 
here constant [i.e., not explicitly space-time-dependent] spinors.) One introduces 
contravariant spinors of either kind and then gives a construction recipe for tensors 
of various ranks in terms of products of spinors. Formally the spinor calculus in 
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space-time rests on the homomorphism of the proper Lorentz group with the 2- 
dimensional complex linear group of transformations A,g with det. = +1 (a 6- 
parameter group). This property is characteristic for the indefinite metric of the 
space-time manifold; it is not valid in the Euclidean case. 

b) In the latter case one can formulate the problem as follows: call & (a = 1, 2) 
an “undotted covariant spinor’ if it transforms as £,’ = )og Uagts, det ||wag|| = 1 
and Usga* = Ugg ', 1.e., u isa unitary unimodular matrix, depending on three param- 
eters. Call 3 a “dotted covariant spinor’ if it transforms like 9/3 = >on, 3, 
where v is another unitary unimodular matrix, whose three parameters range inde- 
pendently from those of wu over their domain of values. £,, ng are again constant 
spinors. One can now follow from here on the entire argument as, for example, 
given by Laporte and Uhlenbeck,!! with only this difference: £&* (n,*) does not 
transform like a covariant dotted (undotted) spinor but like a contravariant un- 
dotted (dotted) spinor. Examples of a &,: the spinors ~"1,,, 0, + 1/,, 0, and of a 
Na: W 0, 1/20, + 1/  —Thespinor calculus rests in this case on the property of the 4- 
dimensional real proper transformation group being essentially!? the direct product 
of two unitary unimodular groups, a property which in its turn does not hold for the 
Lorentz group. 

Erratum.—In the paper referred to in footnote 1 (p. 486, line 25), read 


—> —)> <p =p 


Ra tha KK -®, forks, oP ie — =O 


1A. Pais, these PRocEEDINGS, 40, 484, 1954. 

* F. Klein and A. Sommerfeld, Theorie des Kreisels (Leipzig: Teubner Verlag, 1898), Vol.2, 495. 

’ This was first noted by F. Hund, Z. Physik, 51, 1, 1928. The 4-dimensional spherical har- 
monics expressed in terms of the variables of eq. (4) below are also used by V. Bargmann, Z. 
Physik, 99, 576, 1936, in his discussion of the hydrogen atom. 

'Cf., e.g., Higher Transcendental Functions, Vol. 2, especially sec. 11.6 (Bateman Manuscript 
Project [New York: McGraw-Hill, Book Co., 1953]). I want to thank Professor A. Erdélyi for 
referring me to this work. 

5 We use the following definitions of the Legendre functions: 


m 


Pz) = a at (x? — 1F,, Py, m (2) = (— 1)" (1 — 2?)””2 = P,(2). 

°K. Schrédinger, Pontif. Acad. Sct. Commentationes, 2, 321, 1938. 

7 Summation is implied over tensor indices occurring twice. 

’|{m] = largest integer contained in m. With this convention, all f’s are real. 

* Closely related to this is the following theorem: If we rotate the two-dimensional sphere on 
which the normal spherical harmonics Yim are defined, then the Yim are linearly transformed. 
The transformation coefficients are the 4-dimensional spherical harmonics with integral k. Cf. 
Higher Transcendental Functions, Vol. 2, sec. 11.7; also H. Hénl, Z. Physik, 89,244, 1934. The har- 
monies with half-integral k occur in the transformation of spherical spinors in 3-space. Professor 
V. Bargmann has pointed out to me that with the help of these considerations one is led to ex- 
pressions for the spherical spinors in 4-space. I am indebted to him for an enlightening discussion 
on this subject. 


"y* and y~ separately are eigenfunctions because the expectation value of 8; is a good quantum 
number, If the latter is not true, one has to work with appropriate linear combinations of y* and 
¥~ with the same 7+, i>. 

0, Laporte and G. E. Uhlenbeck, Phys. Rev., 37, 1380, 1931. 

* Cf. B. L. van der Waerden, Gruppen von linearen T'ransformationen (Berlin: Julius Springer 
Verlag, 1935), esp. sec. 7. 











A REMARK ON THE EXTREMAL DISTANCE OF TWO BOUNDARY 
COMPONENTS 


By Kurt STREBEL 
INSTITUTE FOR ADVANCED STUDY 


Communicated by Marston Morse, June 4, 1954 


1. In a previous paper by the author! the following theorem was proved: 

THEOREM |. A domain G, of the complex z-plane with two distinguished boundary 
components 1, and T; of extremal distance \ < © can always be mapped conformally 
on a domain G,, with the following properties: 

a) The boundary components T°, and V1, (the images of T? and TY}, resp.) are circles 
\w| = ro and |w| = 7 > 1 with radial incisions at most on a set arg w. of linear 
measure zero. 

b) Every compact subset of the rest of the boundary is a set of radial slits.” 

c) X = (2r)— log (11/10). 

The mapping is uniquely determined up to a constant factor. 

The purpose of this paper is to characterize the domain G,, by an extremal prop- 
erty expressed by means of extremal length and to give a new proof of the unique- 
ness of the mapping based on this characterization. 

2. We need lemma 12 of Strebel! in a slightly generalized form: 

Let the domain G, be bounded by two analytic Jordan curves I) and I, a com- 
pact set of radial slits T, and a finite set of continua C of total angular measure 
Ag. The curve I> separates G, from a sufficiently small circle | 2| = ro, the curve 
I’; from a sufficiently large circle |z| = 7: > 7. By U(¢) we again denote the logar- 
ithmie length of the set of radial stretches arg z = ¢g lying between Ip and I’; and 
joining these two curves. Then we have, for the extremal distance of Ip an Ii 


with respect to G,, 
1 i 1 
> le, (1) 
r1— Jel(y) © 


where £ is the open set of directions which do not meet C. 

If the radial projection of the set I has linear measure zero, the estimate is im- 
mediately given by familiar arguments. In the general case we need an approxi- 
mation similar to the one used in the paper referred to above.! The do- 
main G,* has now to contain the continua C, which means that it is bounded by the 
=r, |z| = m1, and some arbitrarily chosen Jordan curves I’’ surrounding 


«& 








circles | z 
the set I. In the integration we have to exclude not only small angular intervals 
around the directions tangent to I'y or I, but also the angular intervals containing 
the continua C. The former can be chosen arbitrarily srhall, but not the latter. 
That is why, in the result, we now have only the integration over the complementary 
set E of these angular intervals.’ 

3. Definition: A domain G,, is called a radial-slit domain if it lies in a circular 
ring ro < | w| <r and separates these two circles and if the extremal distance h of 
| w| = ry and | w| = r, with respect to G, is equal to (2) ~! log (11/ro), 7.€., 18 equal to 
the extremal distance of the two circles with respect to the whole plane. 
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This definition is equivalent to properties a and b of Theorem 1. In fact, assum- 
ing G,, to be a radial-slit domain, we consider the intersection I” of the boundary 
of Gy with 75 + 1/n < w Sm — 1/nand the domain G" bounded by I” and the two 








circles | w| = ryand|w| =r. The extremal distance X, of | w| = ro and | w| =r 
with respect to G” obeys the inequality 


(2x)! log (r:/ro) = X > Ay > (2)! log (11/70) 


and is therefore equal to (27)~! log (r:/ro) for all n. As every compact subset of 
boundary components lying in ro < | w| <r, is contained in a certain I, it has to be 
a set of radial slits, and, because the area of I” has to be zero for all n, condition a 
follows. 

If, on the other hand, a and 6b are fulfilled, we can, for each fixed n, exclude from 
7", by finitely many Jordan curves of arbitrarily small total angular measure, 
those points of I” which belong to the incisions attached to the circles. The rest of 


Yn, 


I” being a compact set of radial slits, we can apply inequality (1) to G”: 


r —1 A r —1} r —1 
2a (toe ) o> he g (10g ) dg > (2r—- 0 (tog ) : 
Yo E Yo Yo 


Letting « — 0, we get 


An = (2r)—' log s 
To 

The extremal distance A, of the intersection of |w| = rm + 1/nand |w) = nm — 
1/n with G,,with respect to the intersection of 79 + 1/n < | w| <r, — 1/n with G, 
is, of course, <X, and > (27)! log (71 — 1/n)/(ro + 1/n) and tends to X as n in- 
creases (see Strebel,! or Wolontis,® p. 591). Therefore, X = (27)! log (r:/ro). 
This establishes the equivalence. 

4. Asis well known, the extremal length of a set of curves | y} in a domain G can 
be defined in the following way: One considers all nonnegative measurable func- 
tions p in the domain, including p = ©, for which the length L,(y) = S +p! dz| is 
defined and > 1 foreach y. The glb of the square of the norm|| p||? = /cJS  p? dx 
dy for all such p is equal to 1/A. The set of all admissible functions p is convex: 
With p; and po, also p = tp, + (1 — t)pe (0 < t < 1) belongs tothe set. If there are 
functions p with a finite norm, there exists a minimal sequence p, (|| Pn | — inf | ll), 
and, according to well-known theorems on L?-functions, this is a Cauchy sequence, 
saree |lon — P| = 0, and there exists a nonnegative L*-function p*, uniquely de- 
termined a.e. with ,“2, ||p, — p*|| = 0, and therefore || p*||? = 1/A. We call it the 
generalized extremal metric. In particular, if there exists an extremal metric in 
the strict sense, i.e., p* admissible and || o*||2 = 1/A < ~, it is uniquely determined 
ae. 

If jy}: ¢ | y} is a subclass of curves y, every p admissible for | y} is admissible 
for {y}i. If, therefore, \{ y}1 = \{ 7} < ©, {7}, has the same generalized extremal 
metric as } yt. In particular, if there exists such a subset with an extremal metric 
pin the strict sense, p is equal a.e. to the generalized extremal metric of | y} and is 
therefore independent of the subset { };. 

5. We consider a domain G, with two distinguished boundary components 
r; and I"! of finite extremal distance \ and all conformal mappings of G, on radial- 
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slit domains G,, where 2 and [! go into the boundary components which contain If 

the two circles, respectively (ro <1). The extremal distance X = (27)~! log (r;/ry) 4h 
of the two cireles | w! = ro and | w| = r, with respect to G,,is > A, because the curves 
yielding form a subset of the set of curves yielding \. If there exists a mapping 

with \ = X, it is uniquely determined up to a constant factor, for the function p = an 
(| w(z)| log (ri/ro))-! is admissible for X and || p||? = 2x (log (r1/ro))~! = 1/d. For 
two such mappings w(z) and w*(z), we have, therefore, | w*(z)| = const. | w(z)| 

a.e., Which yields w* = const. w. is 

We can now state Theorem | in the following way: of 

THEOREM 2. I/f G, is a plane domain with two distinguished boundary components dR 

? and VT! of extremal distance \ < ©, there exists within all conformal mappings of inc 
G, on radial-slit domains G,, ro < | w| <r, where the images T°, and T}, of T° and r 
contain the circles | w| = 7) and | w| = r, respectively, one which minimizes the ratio 

ri/ro. It is uniquely determined up to a constant factor. 

' Kurt Strebel, ‘Die extremale Distanz zweier Enden einer Riemann’schen Flaiche,”’ Ann. Acad. [es 

Sci. Fennicae, Vol. 179, 1955. A¢ 

2 T.e., contained in a circular ring Ry < |w|< R,, the extremal distance of the two circles with res 

respect to the complement of the set being (27)~! log (Ri/Ry). We drop the attribute ‘‘minimal” 

used in Strebel, op. cit., in order to shorten the expression. ee Ev 
3 We can express the result by saying that for the estimate (1) the set of radial slits is removable. 

‘ L.e., there exists a Jordan curve in Gy separating them. Tes 

5 Vidar Wolontis, ‘Properties of Conformal Invariants,” Am. J. Math., 74, 587-606, 1952. In 

red 

sceanaaes of 1 

res 

IRREVERSIBLE SYSTEMS, ENTROPY, AND RIEMANN SPACES l 

Che 

By ALEXANDER W. WUNDHEILER se 

ILLINOIS INSTITUTE OF TECHNOLOGY, CHICAGO, ILLINOIS For 

Communicated by P. W. Bridgman, July 13, 1954 0. 

Summary.—For a certain class of irreversible systems the following can be rigor- 1 

ously proved with the means of conventional thermodynamics. aie 

The “entropy change of the universe” associated with the production of an in- of ¢ 

finitesimal displacement d (whose terminals may be nonequilibrium states) has a giv 

nonnegative minimum R(d), the degradation, when the production is a best one. cha 

If, and only if, R(d) = 0, d is reversible. [R(d)]? = dR? is a quadratic form in the cha 

differentials of the state parameters. A geodesic in the state space with arc ele- sot 

ment dR is the path of a continuous best production (at minimum dissipation) of the I 

change from any of its points to another. col 

For a certain choice of the state parameters E,, (ge 

dR? = (gdE,)?. and 

dy a 

If only £, varies, dE; is the energy change, and this displacement can be best pro- rem 

duced using a single reservoir of temperature 1/¢(d;). Assuming that dH; > 0, bol; 

2q: = o(—d;) — o(d;). 
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If no heat reservoir is used, put ¢ = 0. The displacement d; is reversible when 
g: = 0, and then only. Put 


2f; = o(di) + o(—d,), 
and call 1/f; the temperature of the system for the displacement +d;. Then 
dS = >; fdk; 


is the entropy differential if it is integrable. In any case, if d’ is the displacement 
of the complement of the system in a best production of d, we have dS + d’S’ = 
dR. It is believed that gd; are Bridgman’s postulated “characteristic entropy 
increases’ and that this postulate has been proved. ' 


1. RESERVOIRS: TEMPERATURE AND ENTROPY 


1.1.—A change T in a set of reservoirs [X;] is determined by a set of numbers 
le], the energy increases of the X;. If = [e;], the [ke,] is the &-multiple kT of I. 
A change is cyclic if its complementary (i.e., the residual simultaneous change of the 
rest of the universe) can be zero. 

1.2.—Aaziom of Multiplication: If a Change of Two Reservoirs Is Cyclic, So Is 
Every Positive Multiple of the Change.—Hence either all the cyclic changes of two 
reservoirs occur in the same sense, or they are all reversible (equivalent reservoirs). 
In the first case, the reservoir with the positive energy gain is the colder one. This 
reduces the contents of Clausius’ form of the Second Law to a statement of existence 
of inequivalent reservoirs, and the Kelvin-Planck form to the assertion that a heat 
reservoir is always ‘‘colder’’ than a work reservoir. 

1.3.—Kelvin’s Axiom: For Any Three Reservoirs, a Nonzero Reversible Cyclic 
Change, and Every Multiple of It, Are Possible-—Let i = 1, 2,3, and let no two of the 
reservoirs X,; be equivalent. Let [e;] be a reversible cyclic change of this triplet. 
For any other reversible cyclic change of the triplet, the ratios of the e;'s must be the 
same: otherwise a combination of them would produce reversibly the change 
(0, e, —e], e ¥ 0, which is impossible. 

1.4.—F rigors.—Let ¢; + 0 be three numbers such that (ge) = die: + doe2 + des 
= 0 for all reversible cyclic changes of a reservoir triplet [X,]. Then, in view 
of ¢: + e2 + e; = 0, every cyclic change with (¢e) = 0 will be reversible. For a 
given choice of ¢,, the number (¢e) cannot have different signs for two possible 
changes of the triplet: otherwise, a linear combination of two such irreversible 
changes would produce a nonzero reversible cyclic change. Let the ¢; be chosen 
so that (ge) is nonnegative. 

If two numbers ¢» and ¢; > ¢o are assigned to two fixed reservoirs, X) and X, 
colder than Xo, a frigor ¢’ will be assigned to every reservoir X’ by the equation 
(ge) = 0, and a colder reservoir will have a larger frigor. If Xo is a work reservoir 
and ¢) = 0, 1/¢’ will be the Kelvin temperature if ¢’ > 0. If the “‘seale points,”’ 
¢ and ¢;, are changed, ¢’ undergoes a linear transformation. Hence, if two frigors 
remained unchanged, all of them do so. We shall, henceforth, use the frigor sym- 
bols to designate reservoirs. 

1.5.—Entropy for Reservoirs —If T! = [e;] is a change of a set of reservoirs [¢;], 
of) = }'¢,e; is the entropy change in T. If I is cyclic, then o(f) = 0 if I is re- 
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versible; otherwise, o(I) is positive. This relation is invariant under scale trans- 
formations ¢’ = ad + 8 (they transform a general entropy change o into o’ = 
ac + Be,e = >.e,). 

1.6.—Reduction of Reservoir Changes.—If oo, ¢: are two reservoirs, a reversible 
cyclic change of [¢o, ¢1, ¢:] will remove the gain e; from ¢; without an entropy 
change, and the residual change |e, e;], after this is done for every reservoir ¢,, 
is independent of the scale chosen: every change in reservoirs can be reversibly re- 
duced to a two-reservoir change. If e = doe; = 0, then e = eo = 0. Ife ¥ 0, and 
¢ is unrestricted, then a reversible cyclic change of [@o, ¢1, *], where 


Poo + pili 
* 
ro) = 
a 
will result in a residual change in one reservoir only. 
° 2. FINITE CHANGES OF SYSTEMS 
2.1.—Coentropy, Degradation, Pre-entropy.—A (possibly discontinuous) change 
ab of a system from state a to state b (they may be nonequilibrium states) is regular 
if both ab and ba are producible, so that the residual complementary changes IT and 
[’ occur in reservoirs only, [ and I’ being independent of each other. Since I and 
’ add to a cyclic change, the entropy change of the reservoirs o(f + I’) = o(f)+ 
a(I’) is nonnegative and has a nonnegative lower bound. Since fixing I’ does 
not restrict [, (1) will have a lower bound, cab, the coentropy change in ab, reached 
for the best production of ab. Put 


2Rab = aba + aab, 2Sab = aba — aab. 


The degradation Rab is nonnegative, symmetric in a, 6, and is zero when ab is re- 
versible, and only then. The pre-entropy Sab is skew-symmetric in ab. It is easy 
to show that when ab and a’b’ are simultaneous regular changes of 4 system and 
its complement, the swm of the pre-entropies is equal to the degradation: 


Sab + S’a'b’ = Rab. = R’a’b’: 


this relation shows the primary application of the pre-entropy (or entropy change 
when the entropy function exists). 
2.2.—I nequalities.—It is easy to prove that 


| Sabe| = [Sab + Sbe + Sca| < Rab and Rab + Rbe — Rea. 


If Sabe = 0, the state triplet (a, b, c) is entropic. This is always the case when ab 
is reversible. If all the changes are reversible (in a certain region), then Sab = 
Sb — Sa = Sob — Soa, where o is an arbitrary reference state; thus the Clausius 
entropy theorem is proved without restriction to equilibrium states, infinitely slow 
processes, and systems whose state space can be mapped on a plane (also without 
integration). An entropic region (where all state triplets are entropic) may exist 
without unrestricted reversibility. A simple argument shows that the set (a, , 
c, x) is entropic if (a, b, c) is so, and az is reversible. Jf every continuous one-dimen- 
sional set of states is a possible path of a best production, then Rab + Rbe — Rca = 0 
for every (a, b, c), and entropy exists. 
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3. INFINITESIMAL DISPLACEMENTS OF SYSTEMS 


3.1. Elementary Entropy and Degradation.—lf ab = d is an infinitesimal dis- 
placement, Sab = S(d) and Rab = R(d) are, respectively, odd and even in d, and 
both are homogeneous of first degree in d. Hence (barring singularities) S(d) = 
dS is a (possibly nonintegrable) linear, and |R(d)]? = dR? a quadratic form in the 
differentials of the state parameters. In the state space where dR is a Riemannian 
arc element, a geodesic is the path of the system in a best process (at minimum dis- 
sipation). 

3.2. The coefficients of the two forms can be determined for certain systems if 
the following conventional assumption is made: 

Axiom of the warmest reservoir: There exists a reservoir not colder than any other 
one. 

The ‘‘work reservoir’’ of conventional thermodynamics is one of them. If zero 
is assigned as the frigor of the ‘‘warmest”’ reservoir, all the other frigors will be non- 
negative, and the temperatures 1/@ will be nonnegative and unbounded. 

3.3.—Frigor and Ditropy.—A displacement d is ditropic if both d and —d can 
be best produced using only one (but possibly not the same) reservoir. Assuming 
the energy change dF positive, let ¢(d) and ¢(—d) be the frigors of the two reser- 
voirs, and 


2f = o(—d) + ¢(d), 2g = o(—d) — ¢(d); 


fand g are, respectively, the frigor and the ditropy in the directions +d. (If f # 


0, 1/f is the temperature of the system in the direction +d.) We have by 2.1, 
2d8 o(—d) — a(d) = o(—d)dk + od(d)dE = 2fdE, and 2R(d) = o(—d) + 
o(d) = ¢(—d)dE — o(d)dE = 2gdk. Hence, for a ditropic displacement d, 


I] 


dS = fdk, dR = gdE. 


Such a displacement is reversible if, and only if, the ditropy g = 0. Otherwise, the 
ditropy is positive. The frigor f, although it may be different for different direc- 
tions, still performs the same function as the reciprocal of the temperature in en- 
tropy relations for simple substances. 

3.4.—Ditropic Systems.—A displacement d is ergodic if both d and —d can be 
best produced using an unrestricted number of zero-frigor reservoirs. For an er- 
godic displacement, f = g = 0. A system is ditropic (in a region) if it has an en- 
tropy function S and if its every displacement is a sum of (i.e., the result of suc- 
cessively performed) ergodic displacements and a finite number (1) of ditropic 
displacements, d;. Let the co-ordinate lines in the state space include (at every 
point) N lines D;, every element of which is ditropic, and let the parameter along 
each D; be the energy HE: dH; = d,E. Since entropy exists, it follows from 3.3 
that 


és = ps if dE, 


which identifies the operational meaning of the coefficients of the entropy differen- 
tial. Now let dR? = ¥> g,;dE dE; (i,j = 1,2,3). Along D, only dE, 4 0; hence 


dR? =9,dE,. By 3.3, d:R = gdEj; hence gi = g?. By virtue of symmetry, 
d; + d;, i # j, produce the same degradation dR; hence g;; = 0 fori #7. Now 
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dR? = Yi(gdE;)’, 


and the coefficients of the square degradation differential have been identified. 

3.5.—Heat, Internal Energy, Thermal Parameters.—For ditropic systems these 
terms can be defined with precision. In a ditropic displacement dE = TdS, and 
for an arbitrary one 





dQ = > oe, ds; 
is the heat input of the system. Consider the resolutions of the type: 
dE =dQ+ Dd, p,dv, + dZ = dU + dZ, 


where S,, v,, and Z are independent state parameters. When the number of the v, 
cannot be reduced with preservation of this type, U is the internal energy, and S,, v,, 
and any function of them, are the thermal parameters. 

3.6.—Relations to Onsager’s Theory.—The derivation of the fundamental quad- 
ratic form in 3.4 not only removes the restriction to near-equilibrium states but also 
the phenomenological status of Onsager’s coefficients? which are now expressed in 
terms of the basic concepts of conventional thermodynamics. Moreover, a vexing, 
although inconsequential, difficulty in Onsager’s theory is removed by identifying 
the fundamental quadratic form as the square of the minimum entropy increase of 
the universe (i.e., the square of Bridgman’s “characteristic entropy increase”’') 
rather as the “internal entropy increase” of Onsager’s theory (d,S). The latter 
identification is dimensionally untenable, and analytically impossible if entropy is 
a function of state. 





1 P. W. Bridgman, The Thermodynamics of Electrical Phenomena in Metals (New York: Macmil- 
lan, 1934), pp. 57 ff.; The Nature of Thermodynamics (Cambridge: Harvard University Press, 
1941), pp. 136 ff.; Phys. Rev., 58, 845, 1940. 

2S. R. de Groot, Thermodynamics of Irreversible Processes (Amsterdam: North Holland Publ. 
Co., and New York: Interscience Publishers, 1952), pp. 4, 6, 228; Onsager, Phys. Rev., 37, 405, 
1931; 38, 2265, 1931. 


HIGH-TONE STIMULATION AND HEARING LOSS* 
By Irvine E. ALEXANDER AND THoMaAS F. O’BrRIEN, JR.f 
PRINCETON UNIVERSITY 
Communicated by E. G. Wever, July 27, 1954 


For more than fifty years the problem of tonal localization in the cochlea has oc- 
cupied investigators concerned with auditory phenomena. Its importance lies in 
the fact that the foundation for several theories of hearing rests upon the assump- 
tion that such localization exists, and much experimentation has been undertaken 
to describe the nature of the process. The details have proved to be difficult to fill 
in. 

Stimulation deafness studies have been especially fruitful sources of evidence re- 
garding cochlear action. Many techniques have been used. Some investigators 
employed the conditioned-response method; others used cochlear potentials. 
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Some were especially concerned with cochlear damage determined histologically, 
while many used a combination of these methods. The findings, in general, are 
somewhat favorable to localization of high tones in the basal turn of the cochlea. 
However, no such enthusiasm prevails when the localization of low tones in the 
apex is considered. 

A good part of the evidence for high-tone localization comes from studies in 
which intense pure tones are imposed on the ear and the locus of the resultant im- 
pairment to cochlear structures is noted. In a recent study! the senior author, 
utilizing the electrical potentials of the cochlea, observed a general rather than a 
specific loss in sensitivity to a tone of 5,000 cycles, despite the fact that this same 
tone produced cochlear destruction in the limited region 5-9 mm. up from the 
basal end. This result was in keeping with the findings of Wever and Smith? and 
of Smith.* 

Since there is no general agreement as to the relationship between cochlear po- 
tentials and a behavioral index of hearing, the present study was designed to explore 
the correspondence between an induced local lesion caused by a high tone (5,000 
cycles) and what the organism hears as measured by conditioning techniques. 

Apparatus and Procedure.—Guinea pigs were used as subjects. Four animals, 
with normal hearing as judged by the Preyer reflex, were conditioned to ten tones 
ranging from 200 to 10,000 cycles per second. The conditioned respiratory re- 
sponse method was the basic technique employed to determine threshold sensitivity. 
This method has been used by many investigators, among whom are Horton‘ and 
Kemp.> For a detailed description of the conditioning procedure used here the 
reader is referred to O’Brien.® 

Each animal, wrapped in a pneumatic sleeve, was suspended from the ceiling of a 
soundproof chamber 48 inches deep, 55°/, inches high, and 96 inches long. The 
inner surface of the chamber was lined with acoustic tile. Changes in air pressure 
produced by breathing movements were transmitted to a tambour which had at- 
tached to it a paper flag, which in turn scanned a photoelectric-cell amplifier. 
The amplifier output then activated a Brush electromagnetic penmotor, type BL 
902 A. An electromagnetic event-marker recorded the onset, duration, and ces- 
sation of the stimulus tone. The records were made on paper tape moving at a 
speed of 0.64 inch per second. 

Sounds, in the form of pure tones, were provided by a General Radio oscillator, 
type 1304 A. Output was maintained constant by means of an a.c. output meter. 
Attenuators between oscillator and speaker permitted variation in tonal intensity 
by 1-decibel steps from levels well below to levels well above threshold. An elec- 
tronic switch, included in the circuit, prevented clicks at the onset and cessation of 
the stimulus. The speaker, University type PA 30, was fixed to the ceiling of the 
conditioning box. 

After conditioning thresholds were established and rechecked, each animal was 
anesthetized (dosage, 0.66 cc. Nembutal per 100 gm. body weight) and prepared for 
electrical recordings, which were taken from the round window. The operative 
procedure has been described fully elsewhere.’? A 1-microvolt level audiogram was 
obtained for tones from 200 to 10,000 cycles. What is measured here is the sound 
intensity needed to maintain a constant level of response as frequency is varied. 
As a next step, the ear was stimulated with a 5,000-cycle tone at an intensity of 1,000 
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dynes per square centimeter for a period of time long enough to cause a reduction of 
not less than 40, or more than 65, decibels in the 1-microvolt response to the 5,000- 
cycle tone. Past experience has shown that it is in this range, from 40 to 65- 
decibel loss, that some cochlear damage will be produced without sacrificing the 
entire structure.! After these readings were taken, the nonstimulated ear was de- 
stroyed by inserting a sharp probe through the tympanic membrane and cracking 
through the cochlear wall. 

A minimum period of five weeks was allowed for recovery, after which post- 
stimulation conditioning thresholds were determined. The animals were then 
sacrificed, and the temporal bones were removed and prepared for eventual micro- 
scopic study. The histological procedure has been reported previously.* One 
animal was discarded from the study because of a middle-ear infection which was 
evidently unrelated to its previous experimental treatment. 















TABLE 1 
IMMEDIATE Loss IN ELECTRICAL SENSITIVITY 


Duration of 





























Animal Original Sensitivity* Stimulation Loss 
No. (Dynes/Cm?) (Min.) (Decibels) 
994 0.01 10 65 
995 .014 5 62 
998 0.0032 10 53 
* Sound pressure in dynes per square centimeter necessary to produce a response of 1 microvolt prior to over- 
stimulation. 
TABLE 2 


HEARING Loss IN DECIBELS 


Frequency 


(eps) No. 994 No. 995 No. 998 Group Mean 
200 25 4 6 12 
300 44 28 36 41 
500 51 44 56 50 
700 50 38 48 45 

1,000 57 64 56 59 

2,000 48 60 64 57 
3,000 54 56 50 53 
5,000 56 70 58 61 
7,000 78 88 74 80 
10,000 66 79 80 75 


Results—The effect of overstimulating the ear with a 5,000-cycle tone is wide- 
spread and involves at least the range of tones studied here. The duration of 
stimulation required to produce this impairment and the actual loss as measured at 
this time are presented in Table 1. 

The loss in hearing five weeks or more after overstimulation is shown in Table 2. 
For every animal there was a reduction in sensitivity for all test tones. An examina- 
tion of the group averages shows that the lowest tones were least affected, the tones 
between 500 and 5,000 cycles were more seriously and more uniformly affected, and 
the highest tones of 7,000 and 10,000 cycles were most severely affected. Figure | 
portrays this trend graphically. 

Interindividual effects of the 5,000-cycle stimulation were reasonably uniform. 
For any tone tested, not more than 21 decibels of hearing loss separated the most 
sensitive from the least sensitive ear, and for some tones as little as 6 decibels of 
loss would describe the variation. The individual loss curves are presented in 
Figure 2. 
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The histological results are in keeping with previous findings with regard to the 
locus of injury." * In all cases the portion of cochlear structures damaged by this 
overstimulation was between 4 and 10 mm. from the basal end. This region in- 
cludes the upper basal coil of the cochlea and the lower half of the second coil. 
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Fic. 1.—Mean loss in threshold sensitivity for all animals. Each animal is used as its own 
control. 
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_ Fic. 2.—Individual curves representing loss in threshold sensitivity after stimulation with a 
5,000-cyele tone. 


In one case, No. 995, some mild impairment could be traced to the very beginning 
of the basal turn. In no case, however, was any structural impairment noted in the 
upper third of the cochlea, from approximately 11 mm. to the apex. 
Discussion.—The implications of this study will now be obvious. It has indeed 
been demonstrated many times that overstimulation by a high tone will cause 
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structural changes in the basal end of the cochlea. ‘The inference that has usually 
been drawn from such findings is that a high tone affects only a specific place in the 
cochlea and, further, that it is this specific action that forms the basis of pitch per- 
ception, or, more generally, there is said to be a progressive representation of pitch 
in the cochlea. From our results we would argue that these inferences are not sup- 
ported. There is no question that traumatic effects are localized in some degree, 
but this localization may be tangential to the basic problem of how we hear. 

The argument uses as its main support the finding that restricted cochlear lesions, 
in an area reserved by some for tones roughly between 1,400 and 5,000 cycles,* 
affect all tones between 200 and 10,000 cycles and that there is little difference in the 
magnitude of the loss in the wide range of 500-5,000 cycles. In addition to the 
widespread extent of functional impairment caused by this single tone, we have here 
the strange circumstance that the greatest hearing loss obtains for tones higher in’ 
frequency than the injury-producing tone. It seems unlikely that an organ which 
has been assigned such minute and highly differentiated powers of discrimination 
could act in so gross a fashion. What we are saying is essentially this: over- 
stimulation with high tones will produce lesions in the basal end of the cochlea, and 
there is a rough progression of injury locus toward the apex as the frequency of the 
stimulus is decreased. However, this limited pattern is not reflected in either the 
electrophysiological activity of the cochlea or in the behavioral responses of the 
organism. Hearing loss, measured in either of thése ways, remains reasonably 
uniform. 

In order to account for these results, it would seem more fruitful to invoke an ex- 
planatory principle which would look to structures beyond the cochlea to account 
for pitch discrimination. These results are interpreted as unfavorable to an ex- 
planation of pitch perception on the basis of a simple localization of action in the 
cochlea. It is clear that the primary action of the cochlea can supply only gross 
information as to the frequency of a stimulus. 

Summary.—Specific and localized cochlear lesions were produced in the ears of 
guinea pigs by use of an intense high-frequency tone (5,000 cycles). Loss in sen- 
sitivity for tones between 200 and 10,000 cycles was measured immediately after the 
injury and after a period of five weeks or more. It was found that hearing loss was 
uniform and severe for tones between 500 and 5,000 cycles, least severe for the low- 
est tones, and most severe for tones higher in frequency than the injury-producing 
tone. The results are discussed with regard to cochlear localization and pitch per- 
ception. 

* This research was supported in part by the Office of Naval Research under Contract N6-onr- 
270, T.O. 3, Project NR140-322, and also by Higgins Funds allotted to Princeton University. 

+t Now at Yale University Medical School, New Haven, Connecticut. 
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AQUANTITATIVE THEORY OF VISUAL EXCITATION FOR THE SINGLE 
PHOTORECEPTOR* 


By Conrap G. MUELLER 


COLUMBIA UNIVERSITY 


Communicated by C. H. Graham, July 2, 1954 


Introduction.—The paradigm for experiments! using the single fiber-sense cell 
preparation in the study of vision involves three elements: (1) a source of radiant 
energy, the light stimulus; (2) a sense cell containing an absorbing material that 
converts the radiant energy into a form available to a final element; (3) the nerve 
fiber. Any comprehensive account of what happens in the visual system must in- 
volve at least these three components. The large body of facts and theory avail- 
able with respect to each of these components provides a guide in the formulation 
of visual theory. 

Physical experiment and theory provide a description of the discrete or quantum 
nature of light and the probabilistic nature of radiation and absorption. The rate 
at which quanta are emitted and the time distribution of these events can be 
handled by treating quantum emission as a Poisson variable, be it in one dimension 
(e.g., time) or in two dimensions (e.g., two dimensions of space for discussion of area 
effects). Thus the description of the first step in the visual process is uniquely 
specified for us by the physics of the experimental situation.? 

We are in a less favorable position for understanding what happens after quantum 
absorption in the case of the visual sense cell, but the general data of physical 
chemistry and the recent researches of Wald, Morton, Dartnall, and others provide 
a foundation of fact and theory concerning the photochemistry of visual substances. 
Similar comments apply to the study of activity in single nerve fibers. 

The purpose of the present account is to show that the facts and theory now 
available in the areas cited permit a theory of the activity in the single sense cell 
and nerve fiber in simple visual systems and that this theory accounts for much of 
the experimental data obtained in recording from single optic nerve-sense cell 
preparations. 

General Theory.—The formulation involves three assumptions: (1) Every quan- 
tum absorbed changes one molecule of the absorbing substance (in the present no- 
menclature, from state A to A* or from substance A to A*). The description of 
this stage follows directly from the physics of the stimulus as outlined above. 
(2) Every “changed” or ‘‘activated”’ molecule resulting from this quantum absorp- 
tion will subsequently go to a third state (or substance) A’ and in so doing will 
provide a stimulus to the nerve fiber. The rate at which the second change in state 
occurs and the latter’s time distribution are not uniquely forced on us, because of 
our incomplete knowledge of the several molecular states or substances involved. 
Nevertheless, this stage can be simply formulated, as a first approximation, in terms 
common to many reactions following photon bombardment. The assumption made 
is that the probability that an activated molecule in state A* will go to state A’ 
is the same for all intervals dt and that the probability in the case of any molecule is 
independent of the state of any other molecule. (This is the kind of treatment 
given, for example, to radioactive-decay problems in physics.) It follows, of course, 
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that the average rate for an aggregate of molecules is proportional to the number 
of molecules present. (3) Subject only to the restriction that the second change of 
state occur outside the refractory period of the nerve, each such second change of 
state (A* — A’) initiates a nerve impulse. 

For the purposes of the present formulation, we shall employ the notation of a 
change in state, such as is found in discussions of energy bombardment of simple 
atoms. It is realized that the reaction with which we are dealing is of a different 
order of complexity and that it is not possible to specify the exact energy exchanges 
involved following quantum bombardment. As a guide to the discussion, the 
“states” in our notation may be thought of, for example, as the rhodopsin, lumi- 
rhodopsin, and meta-rhodopsin of the visual cycle as outlined by Wald.* (The 
reader is referred to discussions by Dartnall‘ and Ball, Collins, Morton, and Stubbs® 
as to the energetics involved.) 

Briefly restated, then, the hypothesis assumes that every quantum absorbed 
changes one molecule of substance A to an activated state A*. This molecule sub- 
sequently changes to a third state A’. Each of these second changes of state in- 
itiates a nerve impulse unless the change of state occurs in the refractory period of 
the fiber. 

Specific Soluttons.—The problem for any theory of single-fiber activity in a sen- 
sory system is one of describing the distribution of nerve impulses in time. Thus 
we require a statement of the time of occurrence of each impulse or, alternatively, 
the succession of time intervals between impulses (or their reciprocal, frequency of 
discharge) as a function of time and of such experimental variables as the intensity 
of the source, the duration of exposure, etc. 

First, consider the case of test stimuli of short duration. In this case we assume 
that all the quanta that are going to be absorbed are absorbed within the ex- 
posure period and that this period is small compared to the time constants of the re- 
maining system. The theory will subsequently tell us exactly what we mean by 
“short.” For the present we will leave the term qualitatively specified as in the 
range between 0.001 and 0.01 second. Since we consider the first step, the change 
in state from A to A*, to be coincident with a quantum absorption, what we see 
in the nerve response to stimuli of short duration are the characteristics of the 
second step in the theory. 

For events distributed randomly in time in the sense stated in assumption 2 above, 
the time intervals between events are distributed in such a way that 


=e", (1) 
P,=re", - (2) 


where Ps, is the probability of an interval of length greater than ¢, P, is the prob- 
ability of an event in the interval between ¢ and t + dt, and r is the average rate of 
occurrence of the events. This argument is usually proved for the case where r is 
considered constant over the interval of interest, but it can be generalized to include 
cases where r is a function of time or number of previous events. Since r in our 
case is proportional to the number, N*, of molecules in state A*, equations (1) and 
(2) can be rewritten in terms of N*. For example, equation (1) becomes 


= be 
Ps. =e kN t 


(3) 
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- Equation (2) can be solved for the mean interval in the usual manner, and the 
af average length of interval between events is 
of 
] l 
. sles r kN*° (4) 
le In order to extend the results from intervals between changes of state A* — A’ 
ut to intervals between impulses requires a consideration of the refractory period of the 
” nerve fiber—a period following the initiation of a nerve impulse during which the 
1" probability of obtaining a response to a stimulus is considered to be zero. For the 
i- case involving a refractory period it can be shown that the above equations hold if 
~ we change the time dimension to (t — 7), where 7 is the duration of the refractory 
s° period. Thus equation (3) becomes 
. Ps, = geen (5) 
D- and equation (4) becomes 
n- 
of t=rt+ a (6) 
kN* ° 
ail Consider a situation in which the probability is small of having two changes in 
US state (A* — A’) within an interval r. One example involves the case where the 
ys response to a brief test stimulus contains only two or three impulses, i.e., where 
of kN* is small. In such a case it should be possible to use the number of impulses 
Ly as an indicator of the number of molecules changing state. Since, as indicated 
above, the average rate is directly proportional to the number of molecules in the 
ae activated state, a plot of successive intervals (¢ — 7) should yield a linear function 
ae in the reciprocal of the number of activated molecules at the beginning of each 
al interval and, if our approximation holds, in the reciprocal of the number of im- 
ny pulses remaining in the discharge. The intercept of this linear function should be 
he zero. That this theoretical expectation is in agreement with experimental results 
- is shown in Figure 1. Figure 1 shows a plot of the time interval between impulses 
ee (minus a constant, the refractory period) plotted against the reciprocal of the 
he number of impulses remaining to be given. This figure shows a plot of some data 
of Graham and Hartline. Here the data points represent the average of three 
e, 


measurements on a single fiber. Each of the three responses yielded five impulses. 
Figure 2 shows some unpublished results of Mueller in which the stimulus yielded 
1) seven impulses. 

For situations in which the number of quanta absorbed in a brief test flash is 


2) sufficiently large so that N* is approximately a constant for the first few impulses, 
ma the initial frequency should be described by the time-interval relations given in 
of equation (6). Since we assume that N* is approximately constant and equal to the 
ce number of quanta absorbed, N* becomes proportional to the number of quanta 
de presented to the sense cell, which, for constant exposure time, is proportional to the 
ur intensity of the stimulus. Thus, by substituting and solving for frequency (de- 
nd fined as 1/t), we have 

1 - 
3) fom Te /KI? “) 


where f; is the initial frequency of nerve discharge. 
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Fic. 1.—Time interval between impulses as a function of the reciprocal 
of the number of impulses remaining in a short burst of impulses. Data 


from Graham and Hartline, J. Gen. Physiol., 18, 917-931, 1935. 
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Fia. 2.—Same as Fig. 1. Data from Mueller 
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For large values of intensity this equation obviously has the limiting value of 
1/7, and the function is a monotonically increasing function for increasing intensity. 
A simple way of testing this prediction is to return to equation (6), perform the ap- 
propriate substitution, and take the logarithm of both sides of the equation. This 
yields 

log (t — r) = —log K — log /. (8) 


Thus we should observe a straight line with a negative slope of unity if we plot log 
(t — r) against the logarithm of intensity. This relation is expected to hold for ex- 
posures that are short compared with the time constants of the system. For longer 
exposure times the prediction is complicated by the fact that the density function 
has two components: one an increasing and the other a decreasing function of time. 
We shall treat this problem in more detail in a moment. Suffice it to say that the 
initial frequency in the long-exposure case is a function in VI rather than in J. 
Thus, for long exposure, equation (8) becomes 


log (( — r) = K — '/2 log J. (9) 


Figure 3 shows data of Hartline and Graham’ for long exposures, covering a range 
of 1-10,000 in intensity. The solid line is the theoretical line with the required 
slope of —0.5. Figure 4 shows some data of Hartline* obtained for four different 
exposure times. The lowest theoretical line describes the data for the two longest 
exposure times, 1 second and 0.1 second. The upper set of data, containing only 
two points, represents the 1-millisecond exposure. Although more data are needed, 
these points are obviously consistent with the slope of unity shown by the upper, 
solid, theoretical line. The middle, dotted, curve is empirical and shows that the 
expected intermediate slope obtains for an intermediate exposure time of 0.01 
second. 

For small values of intensity the number of impulses given in response to a short 
test flash should correspond approximately to the number of quanta absorbed (and 
thus be directly proportional to the number presented). For higher intensity 
values, where the number of activated molecules is large, the rate at which mole- 
cules change state is large, and the nerve discharge should not follow the changes in 
state completely, since many of the molecules will be changing state (A* ~ A’) 
during the refractory period. The probability function for n impulses should cor- 
respond to the compound distribution involving not only the probability of getting 
m quanta absorbed (Poisson) but also the conditional probability that, given m 
activated molecules representing m absorptions, n of these will change state outside 
the refractory period. Thus the theory predicts the probability of responses in- 
volving any number of impulses. As yet no published data are available for testing 
this solution.® 

One aspect of the probability problem just discussed is encountered in analyzing 
the experimental situation in which we measure the average number of impulses 
given in response to short test flashes of differing intensity. In solving for the 
theoretical expectation of the number of impulses in the total discharge, we must be 
concerned with the distribution of time intervals, the parameter of which is chang- 
ing as a function of time. For low intensities the expected number of impulses 
deviates little from the number of absorbed quanta, as indicated above. For high 
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discharge as a function of stimulus intensity. 
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intensities there would be large losses at the beginning of the discharge and less 
fractional loss toward the end of the discharge. What we require is the expected 
number of changes of state that would follow the previous change of state by a time 
interval greater than the refractory period, i.e., we want the expected number of 
intervals between impulses greater than r. 

Consider the case in which N quanta are absorbed. The probability is unity 
that the first molecule to change state yields an impulse. The probability that the 
second change of state initiates an impulse is the probability that the interval be- 
tween the first and second is greater than 7, and we have seen that this probability is 
e"-)" Probabilities may be assigned similarly for an interval greater than 7 
between the nth and the (n + 1)th impulse. It can be shown that the average 
number of impulses expected when N quanta are absorbed is 


E = K(1 — e~*"®), 
where K is a constant involving the terms k and r. 


Figure 5 shows a plot of some data reported by Hartline and MacDonald." 
Here the number of impulses given in response to a short test flash is plotted as a 
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Fic. 5.—Number of impulses in response to a short test stimulus as a function of stimulus in- 
‘ensity. Data from Hartline and MacDonald, J. Cellular Comp. Physiol., 30, 225-254, 1947, for 
the dark-adapted condition. 


function of the intensity of the stimulus. The solid curve is the theoretical curve. 
At this point a question arises concerning the nature of the refractory period. 
Specifically, one would like to know whether the refractory period follows each 
‘mulation (change of state A* — A’) or follows only the occurrence of an impulse. 
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If we assume that a refractory period follows each stimulation (whether an impulse 
is propagated or not), we obtain as an expression for frequency of nerve discharge 


f = Kle~*”, (10) 


where KJ refers to the rate at which the molecules change state. 

Unfortunately, neither the data from nerve physiology using electrical stimuli nor 
the data resulting from optic-nerve recording provide us with an unequivocal answer 
to this question. It develops, however, that the functions in equations (7) and (10) 
are the same to a first approximation, as may be seen by expanding the right-hand 
side of equation (10) in series. Thus equation (10) becomes 

oe 

or 7" + KIr 
which is of the same form as equation (7). The difference between the two formu- 
lations shows up at the high intensities, i.e., large values of AJ. It may be noted 
that an important consequence of the second alternative (represented by eq. [10}) 
is the requirement that we get an ‘‘overloading”’ phenomenon, i.e., the function re- 
lating initial frequency and intensity of test flash passes through a maximum. 
At present there seems to be no direct evidence on this point. The empirical 
functions do not show the “overloading,” but there is some question of whether 
high enough intensities have been used to produce it. An increase by a factor of 
from 10 to 100 in the range of intensities shown, for example, in Figure 2 would per- 
mit a more conclusive test of these two alternative formulations. (It may also be 
noted that the assumption leading to eq. [10] suggests an account of some charac- 
teristics of the response, such as the “silent period,”’ to continuous stimuli of high 
intensity.) 

Now consider briefly the case in which the total time is large during which ab- 
sorptions are occurring. Under these circumstances the number of activated mole- 
cules, N* (and therefore the instantaneous probability of the nerve discharging), is 
not only spontaneously decreasing with time but is also increasing, owing to the 
stimulus being present. As indicated at the outset, this situation may be treated by 
the same stochastic model. 

As an example of an experimental result for which this type of analysis is ap- 
propriate, consider the expected latency of the first impulse for the case where in- 
tensity is constant but where exposure time cannot be considered instantaneous. 
Here we can separate the two components—the increase in instantaneous prob- 
ability density due to continuing stimuli and the decrease in probability density due 
to spontaneous change of state. By dealing with the interval to the first impulse, 
we are dealing only with the first component. In this case the probability of an in- 
terval greater than tis e-/°%°._ If the number of impulses in the total response is 
sufficiently large so that we are assured that the number of quanta absorbed during 
the exposure time is large, we can consider the steplike increase in the instantaneous 
density function to be continuous in time. For a constant intensity of illumination 
throughout the exposure time, the probability density is changing linearly with 
time. In this case we must evaluate the integral ,/o‘k/tdt, where t is the time of the 
first change of state, A* > A’. Evaluating the integral and taking the median 
value of t gives 
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It, = Bb, 
where 
B = —(log 0.5/k). 

Unfortunately, we do not know ¢ exactly, because it includes traversal time for the 
impulse to reach the recording electrodes. However, we know that t = L — T, 
where 7' is a constant and L is the measured latency. Therefore, the latency is re- 
lated to the intensity of the stimulus by the equation 

log (L — T) = 1/2 log B — '/. log I. (11) 


Thus the latency is independent of exposure time as long as it is greater than ¢ 
where ¢ is less than the latency. 

Figure 6 shows the results from a study by Hartline’ on intensity-time relations 
in the single fiber of Limulus. Here we have plotted the logarithm of (L — 7) 
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Fig. 6.—Latency of response as a function of stimulus intensity. See text for ex- 
planation. Data from Hartline, J. Cellular Comp. Physiol., 5, 229-247, 1934. 


against the logarithm of intensity. The lower set of points represents the 1- 
second exposure and has been moved downward on the axis of ordinates by five- 
tenths of a log unit, so as to separate the functions for inspection. The upper set 
of points represents the one-tenth-second exposure. The solid lines are theoretical. 

Obviously, a similar analysis can be made for the interval between the first and 
second impulse with minor modifications, such as taking account of the refractory 
period, eliminating considerations of conduction time, etc. Thus we are led to 
equation (9), already discussed in connection with measures of the initial frequency. 
In the case of latency as well as of initial frequency, the prediction is that the slope 
of the theoretical curve in a plot such as that shown in Figure 6 should be unity for 
very short exposure times. 
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A more detailed analysis of the latency problems shows that for stimuli of inten- 
sity / and for exposure times ¢, equal or less than (L. — 7’), the latency is given by 
the equation 


L= T+ 0.5t. + = ; (12) 
It, 

For exposures longer than L — 7 the latter term is substituted for ¢,. Equation 
(12) obviously reduces to expressions (8) and (9) under appropriate conditions. 
In addition, this equation relates to some data reported by Hartline where latency 
measures are reported for stimuli of constant energy but of varying duration. 
For this experimental case the last term on the right-hand side of equation (12) 
is a constant, and we expect a linear function of slope 0.5. Eight exposure times in 
Hartline’s experiment satisfy the restrictions of being less than L — 7’, and these 
points are shown in Figure 7 with the theoretical line of 0.5 slope. 
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Fic. 7.—Latency of response as a function of the duration 
of a constant-energy test stimulus. Data from Hartline, J/. 
Cellular Comp. Physiol., 5, 229-247, 1934. 


The present formulation emphasizes some of the “‘transient”’ features of the re- 
sponse of single sense cells to photic stimuli. It indicates the extent to which these 
features can be treated by a theory that does not require significant changes in the 
concentration of the absorbing material and/or complex chemical regenerative sys- 
tems in order to yield the quantitative characteristics of the nerve response. No 
attempt has been made to present a comprehensive review of all the data that can be 
treated by this type of formulation. Rather, the paper attempts to indicate some 
of the specific solutions that are possible and for which experimental data exist. 
For example, more extensive treatment can be given to the data of Hartline® on 
intensity-time functions, to the data of Riggs and Graham" on interpolated flashes, 
and to other relevant experimental results. 

Summary.—This article considers some consequences of a theory of visual ex- 
citation derivable from (1) the physical characteristics of the stimulus, (2) a general 
statistical property of chemical systems, and (3) the refractory period of nerve 
activity. Specifically, the experimental problems treated are (a) frequency of re- 
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sponding in short bursts of impulses, (b) initial frequency of responding as a function 
of the intensity of the stimulus, (c) probability of at least n impulses as a function of 
intensity, (d) average number of impulses in response to brief stimuli as a function 
of the intensity of the stimulus, and (e) latency of response as a function of the dura- 
tion and intensity of the stimulus. 


* This work was supported by a contract between the Office of Naval Research and Columbia 
University. 
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THE STRUCTURE OF THE SCHWANN CELL AND ITS RELATION TO 
THE AXON IN CERTAIN INVERTEBRATE NERVE FIBERS* 
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Communicated June 23, 1954 


So far as we are aware, it is a general rule that, in all types of peripheral nerve 
fibers capable of conducting propagated impulses, the axon (axoplasm) is separated 
from the exterior by Schwann cells. In heavily myelinated vertebrate fibers, 
Schwann cell constituents can be recognized chiefly by the nucleus which lies ex- 
terior to the myelin and by the cytoplasmic inclusions and limiting membranes 
visible at the node of Ranvier." 2 In myelinated invertebrate fibers, such as those 
of prawns and shrimps, the nucleus and cytoplasmic constituents lie on the inside 
surface of the myelin;* so far as we know, no electron microscope description of 
these fibers has as yet appeared. Unmyelinated vertebrate fibers, particularly 
the very thin ‘“C”’ fibers, were shown by Gasser”: ‘ to lie within the protoplasm of the 
Schwann plasmodium. The smaller fibers of squid and lobster nerves appear to 
resemble the vertebrate ‘‘C”’ fibers in that a number of fibers share common Schwann 
cells. In the squid giant and medium-sized fibers and in larger lobster fibers the 
Schwann cells lie around the axon, the nuclei indenting the axon. 
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The bioelectrically active surface membrane, whose properties have long been 
the object of study by electrophysiologists, would appear to be a constituent of 
the Schwann cell or to bound this cell. Accordingly, the detailed structure of ff 
this region of the fiber is of interest not only to the morphologist but also to the 
physiologist. Detailed description of our own studies of invertebrate nerve 
fibers, as viewed in thin sections with the electron microscope, will be given else- 
where. In the present brief report we wish to stress the membranous structures 
which occur in the Schwann cell and particularly the axon—Schwann cell interface, 
the special properties of which suggest that it may play an important role in the ff 
active processes of nerve. t 
































METHODS 


The squid giant fibers and the nerves of the walking legs of the lobster were ex- Ff 


cised in chilled sea water and immediately fixed in 1 per cent OsO, in buffered sea i 
water (artificial sea water, according to Hodgkin and Katz> or Harvey®). Imbed-  § 
ding was in normal butyl methacrylate (occasionally mixed with 5 per cent methy| 
methacrylate). Sectioning was done with the Minot microtome as modified by ff 
Geren and McCulloch.?, An RCA Model EMU-2 electron microscope equipped 
with an objective aperture was used. 
RESULTS 

t 


The Schwann cell may be extremely thin (except, of course, in the neighborhood 
of the nucleus). In the squid giant fiber it is of the order of 0.5 wu (0.2-1.0 u) thick. i 
In the lobster fibers the thickness may be similar, although in some cases the cell 
may be only 0.1 uw or less in thickness. The surface-to-volume ratio of these ff 
Schwann cells is obviously very high. In the squid giant fiber the volume of the 
Schwann cell is of the order of | per cent of that of the axon. However, this small 
protoplasmic mass and its limiting membrane may play a highly important role in 
nerve-fiber activity. 

Despite its thinness, the Schwann cell protoplasm contains occasional spheroidal 
mitochondria having structures similar to that to be described below. Other 
structures commonly observed in active cells, such as endoplasmic reticulum 
(either as diffuse lacework or as a tubular apparatus), Golgi material, etc., were not 
conspicuous, although no careful study has thus far been made of the Schwann cell 
protoplasm as such. 

1. Intracellular Dense Membranes.—Within the Schwann cytoplasm of the squid 
giant fiber there are several (3-6) layers, each about 150-250 A thick and showing a 
double-contoured structure with dense osmiophilic borders (Pl. I, Figs. 1-3). A 
similar layer can usually be seen at the axon—Schwann cell interface. The intra- 
cellular layers may appear to be rather convoluted, arising from the single dense 
membrane at the outer surface of the Schwann cell (PI. I, Fig. 3) and, by a tortuous 


Fias. 1-3.—Transverse sections through Schwann cell region of squid giant fibers. Fig. 1, 
note dense, osmiophilic layers (arrows). Fig. 2, continuity of intracytoplasmic membrane with 
limiting membrane at axon—Schwann cell interface. Fig. 3, continuity of intracytoplasmic mem- 
branes with dense layer (single) at outer surface of Schwann cell. Abbreviations: Az, axoplasm; 
B.M., basement membrane; C.7'., connective tissue; C.7'.C., connective-tissue cell; 1, mito- 
chondria; S.C., Schwann cell cytoplasm. 
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course, reaching the axon—Schwann cell boundary (PI. I, Fig. 2). Such cases may, 
in fact, represent the junction of two Schwann cells. Except for the outer bound- 
ary, beneath the basement membrane, the layers are double-contoured. In the 
lobster fibers usually only the double-contoured membrane at the axon—Schwann 
cell interface and the thin membrane at the outer boundary may be seen, although 
occasionally one other intracytoplasmic dense, double-contoured layer may be 
seen. 

It seems reasonable to identify these layers with lipid-protein, myelin-like layers. 
Their borders are osmiophilic, and they have a thickness equal, within the resolu- 
tion available, to those of the myelin sheath as demonstrated by X-ray diffraction 
methods* * but different from the layer thickness observed by electron micros- 
copy." They may be similar in nature to those recently observed within the 
cytoplasm of various types of cells. 

If squid nerve fibers are immersed in an artificial sea water devoid of Ca or Mg 
ions and are fixed in 1 per cent osmic acid in this medium or if, to the fixative, suf- 
ficient versene is added to deionize the divalent cations in sea water and in the nerve 
fiber, the dense-edged double membranes are observed to be spread apart. Thus 
divalent cations are implicated in the structural integrity of the intracytoplasmic 
dense-edge membranes, and the apparent width and shape of such membranes in 
sections is dependent on the local ionic environment. X-ray diffraction studies!! 
of nerve lipid double layers in dispersions of univalent and divalent cations at vary- 
ing ionic strengths have shown marked differences in the effects, on the distances 
between double layers, of sodium or potassium ions as compared with calcium ions. 
Since one can observe, in electron micrographs of fixed and sectioned preparations 
of purified lipids, spacings comparable to those obtained from X-ray diffraction 
data,'? it should soon be possible, by detailed studies of model systems, to clarify 
the molecular composition of the dense-edged lipid-protein double layers as seen 
in the electron micrographs (see also Finean!*). 

2. The Metatropic Reaction.—As a result of a study of the so-called ‘“‘metatropic”’ 
reaction of Géthlin,'4 Bear and Schmitt showed that the relationship between the 
sign and magnitude of birefringence manifested by the exterior region of inverte- 
brate nerve fibers and the refractive index of the immersion media indicated the 
presence of lipid and protein molecules oriented in the same manner as in the 
myelin sheath of vertebrate fibers; a difference was that the lipid was present in 
very much lower concentration. The oriented lipids were thought to exist in a 
“sheath’”’ surrounding the axon. These studies were extended to the squid giant 
fiber by Bear, Schmitt, and Young,'* who interpreted their results to indicate the 
existence of a myelin-containing layer analogous to the ‘‘metatropic sheath” of 
crustacean fibers, “lying immediately next to the axis cylinder, as belonging to 
cells presumably analogous to the Schwann sheath cells of vertebrate fibers. As in 
the fibers of prawns and shrimps and opposed to the relation in vertebrate axons, 
these cells lie inside the myelin-containing layer, being placed between it and the 
axis cylinder.’’ The thickness of the metatropic sheath in the squid giant fiber was 
estimated to be of the order of 3-4 u. In their illustrations (Pl. I, Figs. 1 and 2) 
Bear, Schmitt, and Young indicated that the ‘‘myelin-containing” layer surrounds 
the Schwann cell. 
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Examination of thin sections of invertebrate fibers with the electron microscope 
has greatly extended these observations and has permitted a more accurate descrip- 
tion of the structural basis of the metatropic reaction. As a result of eary studies 
of thin sections of squid giant nerve fibers in this laboratory,” '7 it was observed 
that the Schwann cell is bounded on the peripheral side by a layer 0.1-0.2 uw thick 
which could be readily differentiated structurally from the surrounding connective 
tissue. Because it was located in the region expected of the metatropic ‘‘sheath,” 
this rather amorphous, somewhat dense layer was tentatively known among our- 
selves as the ‘‘M”’ layer. 

Further study caused us to regard the rather amorphous ‘‘M”’ layer as a kind of 
basement membrane of the Schwann cell, having no obvious relation to oriented 
lipid-protein layers or to the metatropic reaction as observed in the squid giant 
fiber. At its outer surface, bounding the connective tissue, this layer has a smooth 
contour, while the inner boundary is markedly contorted (PI. I, Fig. 1) | Lobster 
fibers (of a size in which the metatropic reaction may readily be demonstrated) 
possess no such basement membrane or any other ‘“‘sheath’’-like structure belonging 
to the nerve fiber itself, peripheral to the Schwann cell. What, then, is responsible 
for the metatropic reaction in such fibers? 

Bear, Schmitt, and Young" estimated the thickness of the metatropic sheath to 
be about 3 yu. From electron micrographs of transverse sections of squid giant 
fibers we see that the Schwann cell is but about 0.5 yu in thickness and, outside the 
basement membrane of the Schwann cell, collagen filaments of the connective tis- 
sue are closely packed. Between layers of collagen filaments are also seen the 
thin processes of the elongate connective-tissne cells. If Bear, Schmitt, and 
Young were correct about the thickness of the metatropic sheath, the connective 
tissue would seem to represent the major fraction of the ‘‘sheath.” 

The collagen filaments obviously contribute to the birefringence of the sheath, 
the sign being negative with respect to the radial direction (due to negative form 
birefringence and the positive birefringence of the collagen filaments which are 
oriented tangentially and longitudinally). Dense, double-edged, osmiophilic in- 
tracytoplasmic layers (Pl. IT, Fig. 5, b) are also frequently observed running paral- 
lel with the elongate processes of the connective-tissue cells. Quite possibly 
these, like their equivalents in the Schwann eells, may contribute birefringence 
positive with respect to the radial direction and hence may play a role in the 
metatropic reaction observed in the usual preparation. 

While the above points bear in a practical sense on the metatropic reaction, the 
aspect which is of primary interest to the nerve physiologist is the evidence which 
the highly sensitive polarization optical method throws on the presence of oriented 
lipid and protein molecules in the nerve fiber itself. Because of the thinness of the 
Schwann cell (about equal to the wave length of light) and the very large diameter 


Fics. 4-5.—Fig. 4, transverse section through fibers of lobster leg nerves, showing peripheral 
distribution of mitochondria just under Schwann cell interface. Note outpocketing of Schwann 
cell cytoplasm marked by arrow (1). Enlargement of marked areas is shown in Figs. 7 and 8. 
Fig. 5, a and 0}, transverse section through fibers of lobster leg nerves showing mitochondria clus- 
tered under axon-Schwann interface. Note internal structure in mitochondria. Also note alter- 
nate layers of collagen filaments and connective-tissue cells (one with nucleus). Fiber diameter 
about 40 u. 
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of the giant fiber, it would be impossible to establish the precise locus of the oriented 
molecules by this method. 

However, certain observations help localize the oriented lipid molecules. Bear, 
Schmitt, and Young'® noticed that squid giant fibers which had been immersed in 
| M NaOH showed positive birefringence at the outer edge of the axon, presumably 
owing to radially oriented lipids. This relatively high concentration of alkali not 
only disoriented the connective-tissue constituents but possessed a refractive index 
sufficiently high to reduce or cancel any negative form birefringence. The infer- 
ence is that the optical positivity originates in the Schwann cell. This is supported 
by our own observation that careful treatment of giant fibers with a bacterial pro- 
teinase!® caused disintegration and disorientation of the connective tissue sufficient 
to reduce greatly or to abolish the negative birefringence of the connective tissue. 
Nevertheless, in glycerinated sea water such proteinase-treated fibers manifested 
positivity at their surface, presumably owing to the presence of oriented lipid 
molecules in the dense, osmiophilic layers in the cytoplasm of the Schwann cells 
and at the axon—Schwann cell interface. 

The question may be raised whether there is sufficient lipid material in the 
Schwann cell to produce a readily measurable positive birefringence and to make up 
the osmiophilic layers seen in the Schwann cell. Mitchison™ stated that, if all the 
lipid in the sheath of the squid giant fiber were concentrated in a continuous layer, 
that layer would have a thickness which could not be greater than about 15 A 
and might be as low as 5 A. Mitchison therefore claimed that the interpretation of 
the metatropiec reaction given by Bear et al.» '® must be incorrect and suggested 
that his own interpretation, offered to explain the birefringence of the red cell 
ghost, be applied also to the metatropic sheath of the nerve fiber. 

Mitchison based his calculations concerning the squid sheath on the lipid analyses 
of McColl and Rossiter.2”. Unfortunately, the values quoted in that paper should 
be multiplied by 1,000 because of a typographical error.?! Estimation of the 
amount of lipid present in the Schwann cell from the data of McColl and Rossiter 
depends upon the average thickness of connective tissue estimated to have re- 
mained on the fibers used for the analysis. However, on any reasonable assump- 
tion there would be more than sufficient lipid to make up half a dozen of the layers 
observed2? (assuming that they were composed entirely of lipid). 

It is probable that the layered structures depicted by the electron microscope to 
exist within the Schwann cell would constitute a Wiener mixed body of the type 
originally proposed to explain the metatropic reaction. If we assume that the 
dense-bordered layers paralleling the surface of the film seen in sections of the giant 
fibers have dimensions similar to those existing in life and that they are myelin-like, 
lipid-protein structures having a birefringence of about 0.011 (as found for myelin 
by Schmitt and Bear?*), it may be shown that retardations may be expected which 
are similar to those which are actually observed. 

3. The Axon Surface of the Schwann Cell and Its Relation to Mitochondria.— 
The mitochondria of the nerve fiber resemble those found in the perikaryon‘ and 
in other tissue cells. They are ovoid structures of widely varying size. Indeed, 
in the axon of the squid giant fiber some of the mitochondria (PI. III, Fig. 6) are as 
large as the sarcosomes of insect flight muscle (2-4 »). They also manifest various 
types of internal fine structure which may reflect real differences in architecture or 
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may depend on differences in response to the fixative. The most commonly ob- 
served structural types are those with internal honeycomb structure, those with 
internal lamellae, those consisting of concentric or helical membranes, and those 
with irregularly shaped inclusions (Pl. IT, Fig. 5, a and b, and PI. III, Figs. 6, 8, 
and 10). 

To call these particulates ‘‘mitochondria”’ seems justified on structural grounds; 
definitive proof would require studies of their chemical properties. Such studies 
are under way, employing extruded axoplasm of squid giant fibers. This material 
has unique experimental advantages because no harsh mechanical treatment, such 
as the use of homogenizers and other devices to disintegrate cells, is necessary. 

The differential distribution of the axonal mitochondria near the surface of the 
Schwann cell may have peculiar significance. In squid giant fibers the mito- 
chondria are found in the body of the axoplasm as well. However, in lobster fibers, 
especially the larger ones, and in the medium-sized squid fibers, the mitochondria 
are almost exclusively clustered just below the surface of the Schwann cell (see 
Pl. I, Fig. 4 and Fig. 5, a and b). 

The concept of mitochondria as the power plant of the cell?> has led morphologists 
to associate the presence of large numbers of mitochondria in tissue situations with 
the expenditure of large amounts of metabolic energy. The orientation of the 
sarcosomal mitochondria about the myofibrils is an example of such a case. If 
such a morphological diagnosis of function were applied to lobster and squid nerve 
fibers, one might suppose that the Schwann cell, spread out in a very thin layer 
around the axon, may be the site of high metabolic activity. The highly contorted 
axonal surface of the Schwann cell is in agreement with such a possibility. The 
high surface-to-volume relation in the Schwann cell in itself suggests high energy 
expenditure and the diffusion into the surface film of molecules presumably from 
Schwann cell protoplasm. 

The axon surface of the Schwann cell manifests outpocketings which are highly 
suggestive. Figures 4, 7, 8, 9, and 10 show structures evidently being produced by 
an active process of evagination of the Schwann cell membrane which closely re- 
semble the mitochondria lying immediately below the surface. To speculate that 
these mitochondria may, in fact, have been produced by the Schwann cell would 
require reliance on highly circumstantial evidence. The implications of this pos- 
sibility are sufficiently important to make its careful evaluation urgent. If the 
surface membrane of the Schwann cell actually can participate in the formation of 
mitochondria, it may be presumed to contain spatially organized enzymes such 
as are thought to occur in mitochondria. Such a situation would be of the greatest 
interest to nerve physiologists because of the possibility which it offers of providing 
a means of chemical coupling of energy-giving reactions with those concerned with 
the maintenance of polarization potentials and the propagation of action waves. 
Experiments devised to test the validity of this speculation are in progress. 


Fig. 6, section through squid giant fiber axon, showing mitochondria containing irregularly shaped 
inclusions. Note also the dense, osmiophilic granules (arrows) characteristic of squid axoplasm. 
Fig. 7, higher magnification of region of outpocketing of Schwann cell cytoplasm marked with 
arrow (1) in Fig. 4. Fig. 8, higher magnification of region marked with arrow (2) in Fig. 4. Note 
continuity of mitochondrial membrane with membrane at the axon—Schwann cell interface (arrow). 
Figs. 9 and 10.—Sections through medium-sized squid nerve fibers (ca. 20 » in diameter), showing 
what appear to be mitochondria connected with the Schwann cell cytoplasm. 
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SUMMARY 


Attention is called to the existence of osmiophilic layers in, and at the surface of, 
the Schwann cells of lobster and squid nerve fibers, and the suggestion is made that 
these are of lipid-protein (myelin-like) constitution and that they manifest the 
optical properties described as ‘‘metatropic.”” Attention is called also to the highly 
contorted nature of the axon border of the Schwann cell and to the fact that the 
mitochondria of the axon (particularly in lobster fibers) are preferentially located 
immediately below this interface. A speculation is offered about the possible re- 
lation between the Schwann cell and certain physiological processes in the nerve 
fiber. 
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ERRATA: Fluctuations in the Space Distribution of the Galaxies 


In the article of the foregoing title appearing in these PRocEEDINGs, 40, 541-549, 
1954, the following corrections should be made: 
1. Equation (6) should read: 


N(m’) = op f° [1 + Dr) ]r? dr, 
2. Equation (7) should read: 
A se 10° {m’ + 5 — Mo 

3. The caption under Figure 1, last two lines, should read: 

Pr = e~"”* (solid line); A = 10%" +9 0) 
4. Equation (17) should read: 

q = 10 (m’ + 5 — M1)/5 

5. Equation (19) should read: 


m Mz ’ 
, wp 

N et Se M f 2 . . 
N(m’) M.— M, { dM 3 dr [1 + D(r)] 


6. Page 545, ninth line from the bottom, equation at the end of the line, should 
read: 
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Erratum: A Nonenzymatic Illustration of “Citric Acid Type’? Asymmetry: The 
Meso-carbon Atom 


In the article of the foregoing title appearing in these PROCEEDINGS, 40, 499-508, 
1954, the following correction should be made: 

The words ‘‘reaction (4)’’ in the third line from the bottom of page 506 should 
read “reaction (3).”’ 


H. E. Carrer 
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